24小时热门版块排行榜

>论坛更新日志 (3209)
>虫友互识 (467)
>导师招生 (119)
>休闲灌水 (113)
>文献求助 (92)
>硕博家园 (58)
>考研 (40)
>论文投稿 (38)
>基金申请 (34)
>考博 (32)
>博后之家 (30)
>论文道贺祈福 (28)
>公派出国 (26)
>教师之家 (22)
>找工作 (20)
>健康生活 (19)

返回列表

当前主题已经存档。

jianchaoyv2651

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 131
帖子: 29
在线:
虫号: 657814
注册: 2008-11-19

[交流] 【求助】fortran程序出错了

一个计算积分的程序如下:
program integral_1
integer,parameter::nmax=12
real*8,parameter::pio2=6.2831852
external Func
real*8::A=0.0,B=pio2,s
write(*,'(1x,a)') 'integral of Func with 2^(n-1) points'
write(*,'(1x,t7,a,t16,a)')'n','Approx.integral'
do i=1,nmax
call TRAPZD(Func,A,B,s,i)
write(*,'(1x,i6,f12.6)')i,s
end do
end program integral_1

Function Func(x)
real*8::x
real*8::Func
Func=sin(x)**4/x**3
end function Func

subroutine TRAPZD(Func,a,b,s,n)
integer,parameter::k1=selected_int_kind(9)
integer(kind=k1)::tnm
real*8::del,sum,x,a,b,s
integer::n
if(n==1)then
  s=0.5*(b-a)*(Func(a)+Func(b))
  else
  tnm=2**(n-1)
  del=(b-a)/tnm
  x=a
  sum=0.0
  do j=2,tnm
x=x+del
sum=sum+Func(x)
  end do
s=0.5*(Func(a)+Func(b)+2.0*sum)*del
end if
end subroutine TRAPZD
运行结果如下:
integral of Func with 2^(n-1) points
   n       Approx.integral
   1 NaN
   2 NaN
   3 NaN
   4 NaN
   5 NaN
   6 NaN
   7 NaN
   8 NaN
   9 NaN
   10 NaN
   11 NaN
   12 NaN
Press any key to continue
我感觉是计算范围越界,请高手指点!!

回复此楼