31 SAVE aln2, asnhcs, nterms, xmax, sqeps, first
32 DATA aln2 /0.6931471805 5994530942e0/
33 DATA asnhcs( 1) / -.1282003991 1738186e0 /
34 DATA asnhcs( 2) / -.0588117611 89951768e0 /
35 DATA asnhcs( 3) / .0047274654 32212481e0 /
36 DATA asnhcs( 4) / -.0004938363 16265361e0 /
37 DATA asnhcs( 5) / .0000585062 07058557e0 /
38 DATA asnhcs( 6) / -.0000074669 98328931e0 /
39 DATA asnhcs( 7) / .0000010011 69358355e0 /
40 DATA asnhcs( 8) / -.0000001390 35438587e0 /
41 DATA asnhcs( 9) / .0000000198 23169483e0 /
42 DATA asnhcs(10) / -.0000000028 84746841e0 /
43 DATA asnhcs(11) / .0000000004 26729654e0 /
44 DATA asnhcs(12) / -.0000000000 63976084e0 /
45 DATA asnhcs(13) / .0000000000 09699168e0 /
46 DATA asnhcs(14) / -.0000000000 01484427e0 /
47 DATA asnhcs(15) / .0000000000 00229037e0 /
48 DATA asnhcs(16) / -.0000000000 00035588e0 /
49 DATA asnhcs(17) / .0000000000 00005563e0 /
50 DATA asnhcs(18) / -.0000000000 00000874e0 /
51 DATA asnhcs(19) / .0000000000 00000138e0 /
52 DATA asnhcs(20) / -.0000000000 00000021e0 /
56 nterms =
inits(asnhcs, 20, 0.1*r1mach(3))
57 sqeps =
sqrt(r1mach(3))
63 IF (y.GT.1.0) go
to 20
66 IF (y.GT.sqeps)
asinh = x*(1.0 +
csevl(2.*x*x-1., asnhcs,nterms))
function inits(OS, NOS, ETA)
OCTAVE_EXPORT octave_value_list or N dimensional array whose elements are all equal to the base of natural logarithms The constant ex $e satisfies the equation log(e)
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)
OCTAVE_EXPORT octave_value_list return the value of the option it must match the dimension of the state and the relative tolerance must also be a vector of the same length tem it must match the dimension of the state and the absolute tolerance must also be a vector of the same length The local error test applied at each integration step is xample roup abs(local error in x(i))<
octave_value sqrt(void) const