double precision function dsindg (x)
c june 1977 edition.   w. fullerton, c3, los alamos scientific lab.
      double precision x, raddeg, dsin
c      external dsin
      data raddeg / 0.0174532925 1994329576 9236907684 886 d0 /
c
      dsindg = dsin (raddeg*x)
c
      if (dmod(x,90.d0).ne.0.d0) return
      n = dabs(x)/90.d0 + 0.5d0
      n = mod (n, 2)
      if (n.eq.0) dsindg = 0.d0
      if (n.eq.1) dsindg = dsign (1.0d0, dsindg)
c
      return
      end