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