先跟各位大大前輩問聲好 只是小妹我們這組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)