The Museum of HP Calculators

HP Forum Archive 13

[ Return to Index | Top of Index ]

Message #1 Posted by hugh steers on 22 May 2003, 2:11 p.m.

i see luiz has an article on trigs for the 16c (possibly also the 17bii). sorry, i missed this thread.

however, on the subject of inverses, i spent a while some time ago developing a fast arctan (but only single precision accuracy) using a limited pade approximation. not sure whether this is accurate or compact enough to be useful when converted to a calculator program, but here it is:

(sorry about the bad formatting, it gets mangled)

#define M_PI 3.141592654f

#define A0X 0.2388229612f

#define A1X 2.445205396f

#define B1X 3.943529789f

#define A2X 1.314747223f

#define B2X 1.798249626f

float atan6(float x) { /* correct to 3e-7. 4 multiplies/divides */ float x2; float v; x2 = x*x; if (x > 1.0f) v = M_PI/2 - (A0 - A1/(x2 + B1 - A2/(x2 + B2)))/x; else if (x >= 0.0f) v = x*(A0X + A1X/(x2 + B1X - A2X/(x2 + B2X))); else v = -atan6(-x); return v; }

[ Return to Index | Top of Index ]

Go back to the main exhibit hall