HP HP-15C - Continuous Memory

60 Section 2: Working With \7j]
The HP-15C readily handles this integral. Don't worry about
\/l — u2 as u approaches 1 because the figures lost to roundoff
aren't needed.
Application
The following program calculates the values of four special
functions for any argument x: ^
i cx
P(x) = J e~* /2dt (normal distribution
2ir ~°° function)
Q(x) = 1 — P(x) = J e~f /2dt (complementary normal
distribution function)
2 Cx 2
erf(^) = ——I e t dt (errorfunction)
^^:*/o
2 r°° ,2
erfc(x) = 1 — erf(^) = —j=J e l dt (corriplementary error
V77- x function)
The program calculates these functions using the transformation
u — e'1 whenever|x\ 1.6 .
The function value is returned in the X-register, and the
uncertainty of the integral is returned in the Y-register. (The
uncertainty of the function value is approximately the same order
of magnitude as the number in the Y-register.) The original
argument is available in register R0.
^
The program has the following characteristics:
• The display format specifies the accuracy of the integrand in
the same way as it does for |7T| itself. However, if you specify
an unnecessarily large number of display digits, the
calculation will be prolonged.
• Small function values, such as Q(20), P(-20), and erfc(10), are ^
accurately computed as quickly as moderate values.
'-