On Thu, 18 Feb 2010 at 04:02AM -0800, eliot brenner wrote: > In the meantime, any suggestions for a work-around (besides write the > entire program in PARI or C)?
I would use mpmath, which (for the one example I tried) seems to be agreeing with Pari: sage: import mpmath sage: mpmath.besselk(5*I, 200) mpc(real='1.151596925536028e-88', imag='0.0') sage: mpmath.mp.dps = 50 sage: mpmath.besselk(5*I, 200) mpc(real='1.1515969255360280546036312473791321981420160724815173e-88', imag='0.0') Try mpmath.besselk? for more info, and see the mpmath documentation on the web -- it's quite nice. Dan -- --- Dan Drake ----- http://mathsci.kaist.ac.kr/~drake -------
signature.asc
Description: Digital signature
