| 查看: 1471 | 回复: 4 | ||
fgwei2000铜虫 (正式写手)
|
[求助]
我编的Simpson积分法fortran程序给不出结果，大侠们看看哪里出了问题？
|
|
如题：我编的Simpson积分法fortran程序给不出结果，大侠们看看哪里出了问题？ 线性函数f=x在区间【0,1】上的积分，结果0.5.[ Last edited by fgwei2000 on 2011-10-29 at 23:07 ] |
回复此楼» 猜你喜欢
实验室接单子
已经有7人回复
-
假如你的研究生提出不合理要求
已经有11人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 各位大侠，我们来一起讨论下菌种退化问题，怎么防止，书上说的不适用
已经有12人回复
- ulipad 怎样让运行结果显示在消息栏中
已经有6人回复
- 在linux下编FORTRAN程序
已经有14人回复
- 求助Fortran编程
已经有8人回复
- 求助计算机高人，编写个Fortran程序
已经有13人回复
- 请教一个fortran小程序编译出错的问题，谢谢
已经有9人回复
- 如何用fortran编写高斯白噪声程序
已经有6人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】请问一下fortran的可视化编程
已经有7人回复
- 【分享】C/C++常用算法设计方法,fortran其他语言编程也可以参考下
已经有8人回复
dingxb
金虫 (正式写手)
迷途书虫
- 应助: 12 (小学生)
- 金币: 1674.5
- 帖子: 482
- 在线: 59.5小时
- 虫号: 34085
- 注册: 2004-01-06
- 性别: GG
- 专业: 原子和分子物理
2楼2011-11-01 14:54:30
dingxb
金虫 (正式写手)
迷途书虫
- 应助: 12 (小学生)
- 金币: 1674.5
- 帖子: 482
- 在线: 59.5小时
- 虫号: 34085
- 注册: 2004-01-06
- 性别: GG
- 专业: 原子和分子物理
【答案】应助回帖
fgwei2000(金币+2): 2011-11-01 22:10:32
|
建议参考下列代码 Program Simpsons implicit none real, Parameter :: a=0,b=1 integer :: n_points integer :: k,i real :: approx real :: del_x, height, area, t,x real :: n_points_1 Write(*,*)'Enter number of points' Read(*,*)n_points n_points_1 = real(n_points) del_x = abs(b-a)/(n_points_1-1.) approx = 0 t = a do k=1, n_points-1 x = t + del_x/2 area = del_x/6 * ( 4*functio_n ( x ) + functio_n(t) + functio_n(t+del_x) ) t = t + del_x approx = approx + area end do Write(*,*)approx Contains Real Function functio_n(x) Real, Intent(in) :: x functio_n=exp(-x) end function end program |
3楼2011-11-01 14:59:13
fgwei2000
铜虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1688.4
- 散金: 220
- 红花: 3
- 帖子: 349
- 在线: 438.1小时
- 虫号: 1265787
- 注册: 2011-04-14
- 专业: 核物理
4楼2011-11-01 22:09:10
fgwei2000
铜虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1688.4
- 散金: 220
- 红花: 3
- 帖子: 349
- 在线: 438.1小时
- 虫号: 1265787
- 注册: 2011-04-14
- 专业: 核物理
5楼2011-11-01 23:16:45