At 02:55 PM 9/29/00 +0100, Daniel wrote:
>Scott,
>
>thanks for the Knuth's books tip! I will check these out.
>
>I'm currently working on the algorithm. I will keep you
>up to date when I've made any progress.

The standard method for x^y (x, y real) in high-level languages is via
scaling to a number between 0 and 1, then using a rational Chebyshev
expansions to calculate exp(a) and log(b). (x^y = e^(y ln x)

If x and y are integer, or y is integer, much simpler methods exist.

================================================================
Robert A. LaBudde, PhD, PAS, Dpl. ACAFS  e-mail: ral@lcfltd.com
Least Cost Formulations, Ltd.            URL: http://lcfltd.com/
824 Timberlake Drive                     Tel: 757-467-0954
Virginia Beach, VA 23464-3239            Fax: 757-467-2947

"Vere scire est per causas scire"
================================================================

--
http://www.piclist.com hint: PICList Posts must start with ONE topic:
"[PIC]:" PIC only "[EE]:" engineering "[OT]:" off topic "[AD]:" ad's