NAME
        IEEEDPSinh -- compute the hyperbolic sine of a floating point number
   SYNOPSIS
          x   = IEEEDPSinh(  y  );
        d0/d1              d0/d1
        double  x,y;
   FUNCTION
        Compute hyperbolic sine of y in IEEE double precision
   INPUTS
        y - IEEE double precision floating point value
   RESULT
        x - IEEE double precision floating point value
   BUGS
   SEE ALSO
        ieeedpcosh(), ieeedptanh()