- 116 名前：デフォルトの名無しさん [2006/12/17(日) 09:35:52 ]
- 下記のプログラムについて、fortran77で行うと実行できないのです＞＜
何がおかしいですかね＞＜
よければ、教えていただけませんか？
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
