
可以自己写一个积分函数调用就够了,我可以帮你写一个简单的:C为积分下限,d积分上限,eps误差范围,s2返回值,f被积分的函数,n循环次数,这里f已经举了一个例子。
subroutine simp(c,d,eps,s2,f,n)
h=d-c
t1=h/2(f(c)+f(d))
n=1
5 s=00
do k=1,n
s=s+f(c+(k-05)h)
end do
t2=t1/2+h/2s
s2=t2+(t2-t1)/3
if(n/=1) goto 20
15 n=n+n
h=h/2
t1=t2
s1=s2
goto 5
20 if(abs(s2)<=1) e=abs(s2-s1)
if(abs(s2)>=1) e=abs((s2-s1)/s2)
if(e>=eps) goto 15
return
end
function f(x)
a=20
b=10
f=sqrt(a2(sin(x))2+b2(cos(x))2)
end
当然要学习最新版本了,95是77的升级版!77里面很多语言相比95略显笨拙
FORTRAN语言的本身计算能力相当出众,尤其是里面的数组、复数计算等,都是工科数值计算的最佳选择
我可以帮助你,你先设置我最佳答案后,我百度Hii教你。
内部库函数就去编译器的帮助文件里查,从外部导入的就从相应的帮助文件查。一般来说,intel的编译器都集成有imkl库函数,另外免费的库函数有lapack和blsk之类的,都有相应的说明。
对于fortran77,可以做下面的事:
1、程序封装:fortran程序可以直接生成exe文件,这就是封装。
2、fortran库:你的过程或者函数经过调式无误后,其编译生成的obj文件可以直接并入到fortran库中去。
后面的f90、f95是否有同样的功能我没试过 ! 你可以查查有关资料。
欢迎分享,转载请注明来源:内存溢出
微信扫一扫
支付宝扫一扫
评论列表(0条)