先跟各位大大前辈问声好 只是小妹我们这组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)