| 查看: 1442 | 回复: 4 | ||
fgwei2000铜虫 (正式写手)
|
[求助]
我编的Simpson积分法fortran程序给不出结果，大侠们看看哪里出了问题？
|
|
如题：我编的Simpson积分法fortran程序给不出结果，大侠们看看哪里出了问题？ 线性函数f=x在区间【0,1】上的积分，结果0.5.[ Last edited by fgwei2000 on 2011-10-29 at 23:07 ] |
回复此楼» 猜你喜欢
职称评审没过，求安慰
已经有35人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
回收溶剂求助
已经有6人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
申请26博士
已经有5人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 各位大侠，我们来一起讨论下菌种退化问题，怎么防止，书上说的不适用
已经有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