00001 subroutine xdatanh (x, result) 00002 external datanh 00003 double precision x, result, datanh 00004 result = datanh (x) 00005 return 00006 end