!-----  Ex4_3: Double precision   -----
      REAL(8) :: deg, rad, s, c
      REAL(8), PARAMETER :: pi = 3.141592653589793_8
      PRINT*,'Input angle variable in degree:'
      READ*, deg;  rad = pi*deg/180.0_8
      s = SIN(rad);  c = COS(rad)
      PRINT*, 'SIN=', s, ', COS=', c
      PRINT*, 'COS^2+SIN^2=', s**2+c**2
      END