>>390 もう少し自力でがんばれｗ PROGRAM DAIKEIHOU C gakusekibangou NAME C sekibun.f conded at 29 May DIMENSION NJUN(3) DATA NJUN/10,50,100/ A = -1.0 B = 1.0 WRITE(6,100) A,B 100 FORMAT(5X,'SEKIBUN FOR 1/SQR(PI)*EXP(-X*X) A', F5.2, 3X, 'B',F5.2) DO 10 K = 1, 3 N = NJUN(K) H = (B - A) / N S = ( FX(A) + FX(B) ) / 2.0 DO 20 J = 1, N - 1 X = A + H * J S = S + FX(X) 20 CONTINUE S = S * H WRITE(6, *) 'BUNKATSU = ', NJUN(K), 'MENSEKI =', S 10 CONTINUE STOP END C FUNCTION FX(X) PARAMETER (PI = 3.14159265) FX = EXP(- X * X) / SQRT(PI) RETURN END >>392
