! double precision
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