//龙贝格积分
procedure QROMB(A, B,miu,xi:real; var SS:real);
CONST
EPS = 0.000001; JMAX = 20;
VAR
S,H:ARRAY[0..21] of real;
JMAXP,J,K,KM:integer; DSS:real;
begin
JMAXP:=JMAX + 1;//何用??????????
K:=5;
KM:=K - 1;//何用??????????????
H[1]:=1;
For J:=1 To JMAX do
begin
TRAPZD(A, B,miu,xi, S[J], J);
If J >= K Then
begin
POLINT(H, S, K, 0, SS, DSS);
If (AbS(DSS) < (EPS * Abs(SS))) Then Exit;
end;
S[J + 1]:=S[J];
H[J + 1]:=0.25 * H[J];
end;
ShowMessage('Too many steps.');//没有规定循环次数啊,这个干吗的????
end;
delphi常用数值算法集中133页的龙贝格积分是不是有问题啊?
- 写回答
- 好问题 0 提建议
- 追加酬金
- 关注问题
- 邀请回答
-
0条回答
悬赏问题
- ¥15 飞机曲面部件如机翼,壁板等具体的孔位模型
- ¥15 vs2019中数据导出问题
- ¥20 云服务Linux系统TCP-MSS值修改?
- ¥20 关于#单片机#的问题:项目:使用模拟iic与ov2640通讯环境:F407问题:读取的ID号总是0xff,自己调了调发现在读从机数据时,SDA线上并未有信号变化(语言-c语言)
- ¥20 怎么在stm32门禁成品上增加查询记录功能
- ¥15 Source insight编写代码后使用CCS5.2版本import之后,代码跳到注释行里面
- ¥50 NT4.0系统 STOP:0X0000007B
- ¥15 想问一下stata17中这段代码哪里有问题呀
- ¥15 flink cdc无法实时同步mysql数据
- ¥100 有人会搭建GPT-J-6B框架吗?有偿