/* arcsinh */
new; cls;
print arcsinh(0.1);



proc arcsinh(x);
   retp(ln(x+sqrt(x^2+1)));
endp;