RE: HP48-HP50G : Error functions ERF(z), FRESNEL integrals, PROBIT, etc.

[Gil]

New version 11b

All the explanations relative to the argument were somewhat changed.


Subroutine inv¦z was completed (in bold) for special cases :

inv¦z
\\<<
DO \\-> p z
\\<< z DUP ERF 4 ROLLD 3 DROPN p - DUP \\-> f
\\<< '2/\\v/\\pi*EXP(-z^2)+f*z' / - \\->NUM DUP DUP RE ABS 9 > SWAP IM ABS 9 > OR
IF
THEN ROT DROP \"
Failure
 as
Interm.Newton
result (
a
,
b
)=
\" SWAP 1 RND + \", &
|
a
| or |
b
| \\>=9 not OK in 
ERF
\" + POP DOERR \"To force a sol, try,
for max.z, 
99
 instead
of 
9
, 
here
 & in 
ERF
\" DROP
END p SWAP z
\\>>
\\>>
UNTIL OVER - DUP RE ABS .0000000001 < SWAP IM ABS .0000000001 < AND
END UNROT DROP2
\\>>