| 查看: 2012 | 回复: 7 | ||||
| 当前主题已经存档。 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
heyansheng12345银虫 (小有名气)
|
[交流]
【求助】求助FORTRAN求积分程序
|
|||
| A等于0到X上的积分，用Fortran写出来的程序。 |
回复此楼» 收录本帖的淘帖专辑推荐
 Deform模拟 |
» 猜你喜欢
面上项目filecode邪修
已经有6人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有14人回复
-
2026国自然放榜时间
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
senlia
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 贵宾: 3.521
- 金币: 2178.1
- 红花: 1
- 帖子: 223
- 在线: 6.7小时
- 虫号: 849005
- 注册: 2009-09-16
- 性别: GG
- 专业: 计算数学与科学工程计算
★ ★
0112358(金币+2,VIP+0):辛苦了，欢迎常来 9-25 09:37
0112358(金币+2,VIP+0):辛苦了，欢迎常来 9-25 09:37
|
program f915 implicit none external f1,f2,f3 double precision a,b,simps write(*,*)'a,b' read(*,*)a,b write(*,*)simps(a,b,1d-7,f1) write(*,*)simps(a,b,1d-7,f2) write(*,*)simps(a,b,1d-7,f3) end program function simps(a,b,eps,f) implicit none double precision f double precision a,b,simps,eps double precision h,s1,s2,d,rc,rp double precision x integer m,i h=(b-a)/2.0 s1=(f(a)+f(b)+4*f(a+h))*h/3.0 rc=f(a+h) rp=f(a)+f(b) m=2 do rp=rp+2*rc x=a-h/2.0 rc=0.0 do i=1,m x=x+h rc=rc+f(x) enddo s2=(rp+4.0*rc)*h/6.0 d=s2-s1 if(abs(s2)>=1.0)d=d/s2 if(abs(d) h=h/2.0 m=m+m enddo simps=s2 end function simps function f1(x) implicit none double precision x,f1 f1=1.0/(1.0+x) end function f1 function f2(x) implicit none double precision x,y,f2 y=1.0/3.0 f2=0.5*((1.0+x*x)**y-(1.0-x*x)**y) end function f2 function f3(x) implicit none double precision x,f3 f3=4.0/(1.0+x*x) end function f3 function simps(a,b,eps,f)表示函数f(x)在(a,b)上积分 其它几个程序和函数是用来调试的 不知是否对你有用 |
6楼2009-09-24 09:05:53
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555230
- 散金: 27828
- 红花: 2148
- 沙发: 66666
- 帖子: 1602255
- 在线: 65200.9小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
2楼2009-09-22 23:02:32
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555230
- 散金: 27828
- 红花: 2148
- 沙发: 66666
- 帖子: 1602255
- 在线: 65200.9小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
★ ★
wangen994(金币+1,VIP+0):辛苦了，欢迎常来 9-23 08:29
heyansheng12345(金币+1,VIP+0): 9-23 16:16
wangen994(金币+1,VIP+0):辛苦了，欢迎常来 9-23 08:29
heyansheng12345(金币+1,VIP+0): 9-23 16:16
3楼2009-09-22 23:05:28
heyansheng12345
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 634.1
- 散金: 54
- 红花: 4
- 帖子: 71
- 在线: 4.9小时
- 虫号: 543364
- 注册: 2008-04-10
- 性别: GG
4楼2009-09-23 16:17:05
