[問題]品管上fortran語法上的小問題
先跟各位大大前輩問聲好
只是小妹我們這組team實在是沒辦法了
只好求助與此
希望較懂得fortran語法或迴圈的能幫個忙
我們要寫出模擬一萬筆資料的失控點
前一百筆符合常態 每次抽五個樣本出來
從第一百零一筆後改變分配
並求出失控那點的次數
最後求出一萬筆失控點的平均(品管上稱之為ARL)
下列是我們的語法
(但因為在前一百筆時想看一萬個R-BAR就錯了!因此不列失控後的語法)
parameter(m=10000)
parameter(n=5)
integer iseed,iperm(n)
real u0,s0,sya(n),x(n),z
real rbar,a,u1,s1,no,arl,sumr
external rnset,rnnof,svrgp,rnget
open(unit=2,file='c:\fortran\www.out',status='old')
u0=0
s0=1
sumr=0
Z=100
iseed=1582585565
do 110 j=1,m
do 111 i=1,z
do 112 k=1,n
call rnset(iseed)
x=u0+s0*rnnof()
call rnget(iseed)
112 continue
call svrgp (n,x,sya,iperm)
r=sya(n)-sya(1)
sumr=sumr+r
111 continue
rbar=sumr/z
write(2,*) rbar
110 continue
end
如果說不看continue110後
只列出一萬筆中的第一筆r-bar是正確的
但我們現在遇到的問題就是加入了110continue後
出現的10000個rbar是會累加的
我們真的不知道該怎麼辦了
有哪位大大了解能幫個忙嗎
謝謝您了
--
我不是人
我是隻貓
而且還是兩隻小鬼的媽
我不笨也不呆
只是常常騎車撞到牆
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.136.12.166
※ 編輯: zeromagi 來自: 140.136.12.166 (06/05 22:08)
討論串 (同標題文章)
完整討論串 (本文為第 1 之 2 篇)：