program main implicit double precision(a-z) double precision AI1c COMMON/CS/EQ,GAMA real*8::X(6),w(6) CHARACTER*22 INP,OUT,DATA WRITE(*,'(/A$)') ' INPUT FILE: ' READ(*,'(A)') INP DATA=INP OPEN(4,FILE=INP,STATUS='OLD',ACCESS='SEQUENTIAL') WRITE(*,'(/A$)') ' OUTPUT FILE: ' READ(*,'(A)') OUT OPEN(10,FILE=OUT,ACCESS='SEQUENTIAL') X=(/-0.93246951,-0.66120939,-0.23861919,0.23861919, + 0.66120939,0.93246951/) w=(/0.17132449,0.36076157,0.46791393,0.46791393,0.36076157, + 0.17132449/) READ(4,*) N AFA=0.0 AI1=0.0 u=1 u=-1100 DO 10 I=1,35 10 CONTINUE AI1=1/SQRT(2*3.14159265358979)*2.718281828**(-X(I)**2/2) x(I)=AFA*u FXI=1/SQRT(2*3.14159265358979)*2.718281828**(-X(I)**2/2) Fu=AFA*FXI DO 20 I=1,6 AI1=AI1+w(i)*Fu20 CONTINUE AIAFA=AFA*AI1 FAI=0.5*(AIAFA+1) WRITE(10,*) 'AFA=',AFA,' ','FAI=',FAIC AFA=AFA+0.1 IF(AFA.LT.3.6) GOTO 100C END各位大神能帮我看看么，倒是可以运行，就是结果不对啊，用高斯勒让德公式求标准正态分布，到底错在哪啊？
