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



proc arctanh(x);
   retp(1/2*ln((1+x)./(1-x)));
endp;