照抄的程序不知道哪里不对劲，结果竟然是0

xiaorenwu

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	module newton_cotes

02

contains

03

04	subroutine solve(func,s,a,b,n)

05	implicit real*8(a-z)

06	external func

07	integer :: n,i

08	hstep = (b-a)/n

09

s=0

10

11

do i=1,n

12	c=a+i*hstep

13	d=a+i*hstep

14	call cotes(func,s1,c,d)

15

s=s+s1

16

end do

17	end subroutine

18

19	subroutine cotes(func,s1,c,d)

20	implicit real*8(a-z)

21	external func

22	h=(d-c)/500

23	call func(f1,c)

24	call func(f2,c+h)

25	call func(f3,c+2d0*h)

26	call func(f4,c+3d0*h)

27	call func(f5,c+4d0*h)

28	call func(f6,d)

29	s1=19d0*f1+75d0*f2+50d0*f3+50d0*f4+75d0*f5+19d0*f6

30	s1=s1*5d0*h/288d0

31	end subroutine cotes

32

33	subroutine fun1(f,x)

34	implicit real*8(a-z)

35	f=x**2+dsin(x)

36	end subroutine fun1

37	end module newton_cotes

38

39	program main

40	use newton_cotes

41	implicit real*8(a-z)

42	integer :: n

43	open(unit=11,file='result.txt')

44	write(11,101)

45	call solve(fun1,s,-2d0,2d0,8)

46	write(11,102) s

47	101 format(/,T5,'复合高阶Newton_Cotes法计算数值积分',/)

48	102 format(T5,'积分结果',F15.8)

49	end program main

xiaorenwu

补充一下，这个程序是宋叶志版本fortran科学计算与工程的370数值积分部分的

fcode

[Fortran] 纯文本查看 复制代码

查看代码

打印?

1	c=a+i*hstep

2	d=a+i*hstep

检查一下这两行是否抄错了？c 和 d 不应该一样

岸边的鱼

fcode 发表于 2014-7-24 15:47
检查一下这两行是否抄错了？c 和 d 不 ...

老大就是老大，我检查那本书了，的确是你说的问题，原书上是

[Fortran] 纯文本查看 复制代码

查看代码

打印?

1	c=a+(i-1)*hstep

2	d=a+i*hstep

楚香饭

楼主已失踪