! table of sin and cos
PROGRAM ex5_3

  IMPLICIT NONE
  INTEGER :: k
  REAL, PARAMETER :: rad = 3.1415926E0 / 180.0E0
  REAL :: r, s, c

  PRINT '(1X,A)', 'Šp“x      sin       cos       check'

  DO k = 0, 90, 10
    r = k * rad
    s = SIN( r )
    c = COS( r )
    PRINT '(1X, I3, 2X, 3F10.5)', k, s, c, s ** 2 + c ** 2
  END DO

END PROGRAM ex5_3