蒲丰氏投针问题的模拟过程.doc

1、蒲丰氏投针问题的模拟过程，随机数发生器也是自编的，以供大家参考和提出建议。谢谢。（seed1 和 seed2 最好选择 3 和 5，为了使投针次数达到 1000000，CVF 进行如下设置Project-settings-link- output，将 stack allocations reserve:设为 1000000000） program getpi implicit none real,parameter:a=5,L=4,pi=3.14159 integer:n1,i,counter=0 real,allocatable:R1(:),R2(:)real:theta,x,pi1 wri

2、te(*,*) input the size of the array: read(*,*) n1 allocate(R1(n1) allocate(R2(n1) call random(n1,R1,R2) do i=1,n1x=a*(2*R1(i)-1) theta=pi*R2(i) if(abs(x)0element(k)=rem(temp,base(j);temp=fix(temp/base(j);k=k+1;endres(i,j)= 0;for k=1:length(element)res(i,j)=res(i,j)+element(k)/(base(j)k);endendendresult=res;