随机方向法程序及结果

C==============PROGRAMRANDIRC==============DIMENSIONX(25),GX(50)COMMON/ONE/NSR,KWR,ILI,NPE,NFX,NGRREAD(*,*)N,KG,NSRWRITE(*,1001)WRITE(*,1002)N,KG,NSR1001FORMAT(25X,'========PRIMARYDATA========'/)1002FORMAT(5X,'N=',I4,2X,'KG=',I4,2X,'NSR=',I4)CALLMAISUB(N,KG,X,GX