| 查看: 1791 | 回复: 4 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
fgwei2000铜虫 (正式写手)
|
[求助]
我编的Simpson积分法fortran程序给不出结果，大侠们看看哪里出了问题？
|
|
|
如题：我编的Simpson积分法fortran程序给不出结果，大侠们看看哪里出了问题？ 线性函数f=x在区间【0,1】上的积分，结果0.5.[ Last edited by fgwei2000 on 2011-10-29 at 23:07 ] |
回复此楼» 猜你喜欢
filecode，4个jtjc了
已经有13人回复
-
今天基金会出结果吗？20260819
已经有10人回复
-
时间戳变了，能看出什么问题？
已经有15人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有121人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有6人回复
-
今天放榜吗？
已经有15人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 各位大侠，我们来一起讨论下菌种退化问题，怎么防止，书上说的不适用
已经有12人回复
- ulipad 怎样让运行结果显示在消息栏中
已经有6人回复
- 在linux下编FORTRAN程序
已经有14人回复
- 求助Fortran编程
已经有8人回复
- 求助计算机高人，编写个Fortran程序
已经有13人回复
- 请教一个fortran小程序编译出错的问题，谢谢
已经有9人回复
- 如何用fortran编写高斯白噪声程序
已经有6人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】请问一下fortran的可视化编程
已经有7人回复
- 【分享】C/C++常用算法设计方法,fortran其他语言编程也可以参考下
已经有8人回复
fgwei2000
铜虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1688.4
- 散金: 220
- 红花: 3
- 帖子: 349
- 在线: 438.1小时
- 虫号: 1265787
- 注册: 2011-04-14
- 专业: 核物理
4楼2011-11-01 22:09:10
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
- 专业: 核物理
5楼2011-11-01 23:16:45
