Fortran Coder

查看: 8498|回复: 4
打印 上一主题 下一主题

[数值问题] 照抄的程序不知道哪里不对劲,结果竟然是0

[复制链接]

11

帖子

5

主题

0

精华

入门

F 币
92 元
贡献
40 点
跳转到指定楼层
楼主
发表于 2014-7-24 15:02:18 | 显示全部楼层 回帖奖励 |倒序浏览 |阅读模式
[Fortran] 纯文本查看 复制代码
module newton_cotes
contains

subroutine solve(func,s,a,b,n)
implicit real*8(a-z)
external func
integer :: n,i
hstep = (b-a)/n
s=0

do i=1,n
c=a+i*hstep
d=a+i*hstep
call cotes(func,s1,c,d)
s=s+s1
end do
end subroutine

subroutine cotes(func,s1,c,d)
implicit real*8(a-z)
external func
h=(d-c)/500
call func(f1,c)
call func(f2,c+h)
call func(f3,c+2d0*h)
call func(f4,c+3d0*h)
call func(f5,c+4d0*h)
call func(f6,d)
s1=19d0*f1+75d0*f2+50d0*f3+50d0*f4+75d0*f5+19d0*f6
s1=s1*5d0*h/288d0
end subroutine cotes

subroutine fun1(f,x)
implicit real*8(a-z)
f=x**2+dsin(x)
end subroutine fun1
end module newton_cotes

program main
use newton_cotes
implicit real*8(a-z)
integer :: n
open(unit=11,file='result.txt')
write(11,101)
call solve(fun1,s,-2d0,2d0,8)
write(11,102) s
101 format(/,T5,'复合高阶Newton_Cotes法计算数值积分',/)
102 format(T5,'积分结果',F15.8)
end program main



分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

11

帖子

5

主题

0

精华

入门

F 币
92 元
贡献
40 点
沙发
 楼主| 发表于 2014-7-24 15:14:05 | 显示全部楼层
补充一下,这个程序是宋叶志版本fortran科学计算与工程的370数值积分部分的
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-5-14 10:18

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表