10/05/21 02:12:18
>>124
エラーチェック無し、暴走考慮せず、記述通りで無駄あり、Fortran90だ。
PROGRAM unko
IMPLICIT NONE
REAL :: x, y1, y2, a, h
x = 0.0
DO
y1 = (x + 1.0) * (x - 3.0)
y2 = 2.0 * x - 2.0
h = y1 / y2
IF (ABS(h) < 1.0e-7) EXIT
a = x - h
x = a
END DO
WRITE(*, *) 'result=', h, ': x=', x
STOP
END PROGRAM unko
実行結果
result= 0.0000000E+00 : x= -1.000000