For now the following seems sufficient for what I need to do: sage: import sage.libs.mpmath.all as mpmath sage: w=mpmath.call(mpmath.besselk,5*I,200,prec=100) sage: w 1.1515969255360280546036312474e-88 sage: type(w) <type 'sage.rings.complex_number.ComplexNumber'>
thanks for your help. Eliot On Feb 18, 7:27 am, Dan Drake <[email protected]> wrote: > 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 > < 1KViewDownload -- To post to this group, send email to [email protected] To unsubscribe from this group, send email to [email protected] For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org
