c     
c     Compute approximations of sinhyp(x)
c
      real err1,err2,sinh1,sinh2,sinh3,x
      integer i,n
      n=8
      x=1.0
      do 10 i=1,n
         x=x/10.0
         sinh1=(exp(x)-exp(-x))/2  
         sinh2=sinh(x)  
         err1=abs((sinh1-sinh2)/sinh2 )
         sinh3=((x*x/120.0+1/6.0)*x*x+1)*x
         err2=abs((sinh3-sinh2)/sinh2)
         write(*,200)x,sinh1,sinh2,err1,err2
 10   continue
 200  format(t5,5e14.4)


      end