for i=1:10
    h=10^(-i);
    fp1=(exp(h)-1)/h;
    fp2=sinh(h)/h;
    err1=abs(1-fp1);
    err2=abs(1-fp2);
    fprintf(1,'%11.10f %11.10f %11.10f %11.10f %40.39f \n', h,fp1,fp2,err1,err2);
end