program sqrt2 implicit none real(16) :: x, y, z, y2, y3, f, e integer :: i, j, k REAL(16), PARAMETER :: sqrt05 = SQRT(0.5q0) DO j = 1, 1000 CALL RANDOM_NUMBER(x) x = x * 1000.0q0 ! print *, 'x=', x f = fraction(x) i = exponent(x) ! print *…

PROGRAM test IMPLICIT NONE INTEGER :: i REAL :: x, s, c, t c = 0.0 s = 0.0 DO i = 1, 100 x = 0.01 - c t = s + x c = (t - s) - x s = t PRINT *, s, c END DO PRINT *, s == 1.0, s STOP END PROGRAM test

-LOG(2)=LOG(1/2)=LOG(1-1/2)=1/1(1/2)+1/2(1/2)^2+1/3(1/2)^3+....確認テスト PROGRAM test IMPLICIT NONE INTEGER :: i, j, k, m, n REAL :: x, y n = 10 PRINT *, 2**n * LOG(2.0) x = 0.0 y = 0.0 DO k = 1, n m = 1 DO i = 1, n - k m = MOD(2 * m, k) …