24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

heyansheng12345

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 634.1
散金: 54
红花: 4
帖子: 71
在线: 4.9小时
虫号: 543364
注册: 2008-04-10
性别: GG

[交流] 【求助】求助FORTRAN求积分程序

A等于0到X上的积分，用Fortran写出来的程序。

回复此楼

» 收录本帖的淘帖专辑推荐

Deform模拟

» 猜你喜欢

职称评审没过，求安慰已经有49人回复
26申博自荐已经有3人回复
A期刊撤稿已经有4人回复
垃圾破二本职称评审标准已经有17人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
三无产品还有机会吗已经有6人回复

1楼 2009-09-22 21:10:24

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

senlia

木虫 (小有名气)

应助: 2 (幼儿园)
贵宾: 3.521
金币: 2178.1
红花: 1
帖子: 223
在线: 6.7小时
虫号: 849005
注册: 2009-09-16
性别: GG
专业: 计算数学与科学工程计算

★ ★
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) s1=s2
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)上积分其它几个程序和函数是用来调试的
不知是否对你有用