! 例題 4-3  三角関数 ［倍精度計算 ］
PROGRAM ex4_3

  IMPLICIT NONE
  REAL ( 8 ), PARAMETER :: pi = 3.141592653589793D0
  REAL ( 8 ) :: deg, rad, s, c

  PRINT *, 'Input angle value in degree: '
  READ *, deg

  rad = pi * deg / 180.0D0

  s = SIN( rad )
  c = COS( rad )

  PRINT *, 'sin = ', s, ', cos = ', c
  PRINT *, 'cos^2 + sin^2 = ', s ** 2 + c ** 2

END PROGRAM ex4_3