くだすれFORTRAN(超初心者用)その2
■ このスレッドは過去ログ倉庫に格納されています
0116デフォルトの名無しさん
2006/12/17(日) 09:35:52何がおかしいですかね><
よければ、教えていただけませんか?
dimension ia(100),vw(12),ivw(12),vw2(100)
data am,ix,n/1.0 ,0 ,100/
call poiss(am,ix,ia,n,vw,ivw,vw2)
write(*,*) (ia(i),i=1,n)
stop
end
subroutine poiss(am,ix,ia,n,vw,ivw,vw2)
dimension ia(*),vw(*),ivw(*),vw2(*)
real*8 x,y
data bm/-1.0/
if (am.ne.bm) then
bm=am
vw(1)=exp(-am)
vw(2)=vw(1)*am
m=2.0*am+10.0
do 10 k=2,m-1
vw(k+1)=vw(k)*am/k
10 continue
■ このスレッドは過去ログ倉庫に格納されています