| 查看: 2088 | 回复: 7 | ||||
| 当前主题已经存档。 | ||||
heyansheng12345银虫 (小有名气)
|
[交流]
【求助】求助FORTRAN求积分程序
|
|||
| A等于0到X上的积分，用Fortran写出来的程序。 |
回复此楼» 收录本帖的淘帖专辑推荐
 Deform模拟 |
» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有23人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有22人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有14人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有11人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有10人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有12人回复
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
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555230
- 散金: 27828
- 红花: 2148
- 沙发: 66666
- 帖子: 1602255
- 在线: 65200.9小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
5楼2009-09-23 22:34:56
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
wuchenwf
荣誉版主 (职业作家)
- 应助: 0 (幼儿园)
- 贵宾: 3.433
- 金币: 19427.4
- 散金: 10
- 红花: 4
- 帖子: 3566
- 在线: 1037小时
- 虫号: 398569
- 注册: 2007-06-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
7楼2009-09-24 22:47:46
jove1782
木虫 (正式写手)
小木虫风雅协会代理秘书长
- 应助: 3 (幼儿园)
- 贵宾: 0.006
- 金币: 5540.8
- 散金: 27
- 红花: 3
- 帖子: 982
- 在线: 134.9小时
- 虫号: 570051
- 注册: 2008-06-06
- 性别: GG
- 专业: 电化学
8楼2009-09-27 11:24:01
