主题:两个简单积分的差别
大家好,我在用谭浩强的教科书学习积分,第八章P171矩形法的例子如下:(稍微有点简化,但是结果正确)
PROGRAM MAIN
REAL A,B,H
PARAMETER(N=100)
F(X)=EXP(X)
A=1.
B=0.
H=(A-B)/N
X0=A
F0=F(X0)
S=0
DO 101 I=1,N
SI=F0*H
S=S+SI
X=X+H
F0=F(X)
WRITE(*,*)I,S
101 CONTINUE
END
我在看这个例子之前自己先编了一个,但是结果与上面有很大区别,我查了半天也不知道什么地方出错了,请大家帮帮看看好吗》?下面是我自己编的程序:
program main
REAL A,B,H
INTEGER N
PARAMETER (N=100)
F(X)=EXP(X)
A=1.
B=0.
H=(A-B)/N
s=0
DO 100 I=1,N
X=float(I-1)*H+A
SI=F(X)*H
S=S+SI
PRINT *,I,S
100 CONTINUE
END
PROGRAM MAIN
REAL A,B,H
PARAMETER(N=100)
F(X)=EXP(X)
A=1.
B=0.
H=(A-B)/N
X0=A
F0=F(X0)
S=0
DO 101 I=1,N
SI=F0*H
S=S+SI
X=X+H
F0=F(X)
WRITE(*,*)I,S
101 CONTINUE
END
我在看这个例子之前自己先编了一个,但是结果与上面有很大区别,我查了半天也不知道什么地方出错了,请大家帮帮看看好吗》?下面是我自己编的程序:
program main
REAL A,B,H
INTEGER N
PARAMETER (N=100)
F(X)=EXP(X)
A=1.
B=0.
H=(A-B)/N
s=0
DO 100 I=1,N
X=float(I-1)*H+A
SI=F(X)*H
S=S+SI
PRINT *,I,S
100 CONTINUE
END