主题:[求救]FORTRAN 程式問題~~
[color=0000FF]请网民求救我…………感恩不尽!
我需要用蒙地卡罗计算 PI值,使绝对误差小于0.005,将计算过程
的结果写到档案 results.dat 中。
程序在下面!我写的不够好 不完整~
绝对值 误差 小于0.005 怎写呢
还有一些 问题!在 注解后面..然后 蒙地卡罗是?[/color] program one
real x,y,r,pi,temppi
integer I,hits,hits_sucess
PI=4.0*atan(1.0) [color=008000]!这是什么??ATAN..?[/color] open(2,file='result.dat',status="unknown")
write(2,*)' hits x y r pi abs.error'
write(*,*)'enter the number to hit..'
read(*,*)hits
hits_sucess=0
DO I=1,hits
x=ran( )
y=ran( )
r=1
r=sqrt(x*x+y*y)
IF (r.le.1.0)then
hits_sucess=hits_sucess+1
temppi=4.0*real(hits_sucess)/real(hits)
write(2,*) I,x,y,tempi,abs(pi-temppi)
write(*,*) I,x,y,tempi,abs(pi-temppi)
end if
enddo
stop
end
我需要用蒙地卡罗计算 PI值,使绝对误差小于0.005,将计算过程
的结果写到档案 results.dat 中。
程序在下面!我写的不够好 不完整~
绝对值 误差 小于0.005 怎写呢
还有一些 问题!在 注解后面..然后 蒙地卡罗是?[/color] program one
real x,y,r,pi,temppi
integer I,hits,hits_sucess
PI=4.0*atan(1.0) [color=008000]!这是什么??ATAN..?[/color] open(2,file='result.dat',status="unknown")
write(2,*)' hits x y r pi abs.error'
write(*,*)'enter the number to hit..'
read(*,*)hits
hits_sucess=0
DO I=1,hits
x=ran( )
y=ran( )
r=1
r=sqrt(x*x+y*y)
IF (r.le.1.0)then
hits_sucess=hits_sucess+1
temppi=4.0*real(hits_sucess)/real(hits)
write(2,*) I,x,y,tempi,abs(pi-temppi)
write(*,*) I,x,y,tempi,abs(pi-temppi)
end if
enddo
stop
end