有关矩形法求变上限积分的问题

红街咖啡

implementation
{$R *.dfm}
function f(x:Real):Real;
begin
f:=2*x;
end;
procedure TForm1.Button1Click(Sender: TObject);
var
i,x,y,s,b: Real;
n,j:Integer;
const
a=1;
H=99;
m=100;
begin
Button1.Caption:= '积分结果 ';
j:= StrToInt(edit1.Text); //分成j等分，edit1中先输入的值
i:=(b-a)/j;
y:=0;
for n:=1 to j do
begin
x:= a+(n-1)*i ;
y:=y+ f(x); //即f（a）+f（a+i）+f（a+2i）+···的值
end;
s:=y*(b-a)/j ; //总面积的值
for b:= a to m do
begin
If s=H then
break;
end;
If s=H then
b:=b;
end;

edit1.Text:='';
edit1.Text:= floattostr(b);
end;
end.

红街咖啡

谁能帮忙看下哪里出错了，或者有会矩形法求变上限积分的例子，给我借鉴借鉴，如果矩形法求变上限积分得以解决，必有重谢~~

709023996

没懂你意思，你把要求说一下，重谢不必互相学习

红姬

matlab么？

红街咖啡

红姬 发表于 2013-5-13 22:55

                               
登录/注册后可看大图

matlab么？

不是，Delphi7