第7章-蒙特卡罗方法-(附录)25页.doc

第7章 附录7.2.1 均匀分布随机数例题7.2.1计算程序! rand1.forprogram rand1implicit nonereal rinteger n,c,x,i open(5,file=rand1.txt)n = 32768c = 889x = 13 do i = 1,1000 x = c*x-n*int(c*x/n) r = real(x)/(n-1) write(5,(f8.5) rend doend!rand2.for!program rand2implicit noneinteger, parameter : n=1000integer ix,ireal ropen(5,file=rand2.txt)ix=32765do i=1,n call rand(ix,r) write(5,(f8.6) rend doend