[Matrix-SIG] Power of trigs

Dave Stinchcombe dars@soton.ac.uk
Thu, 19 Mar 1998 11:28:57 +0000 (GMT)


Hi again,
This morning I needed to write a quick and easy root finder, so obviously I
grabbed Numeric Python.

I have a problem though. I wrote the line:
Numeric.power(Numeric.cos(phi),1-2*nus)

where nus is a constant of 0.3. phi is an array holding a list of values
I'm searching through. The problem comes that this line fails when phi is
greater than pi/2. 

Now obviously this is due to the fact I'm trying to calculate a complex (or
at the very least an imaginary) number. But I want to. Can any one give
some suggestions as to how to cope, or in some way make python also believe
in complex numbers.

Yours
Dave