Hello,
I am working with jacobi elliptic functions with both 2d and parametric graphs and am finding they leave the domain when I use these functions under multiplication. I believe the numeric package is running out of memory or the plotter makes some similar error. How do I fix this, or how can it be fixed?
specifically try jacobi_dn()*jacobi_dn() with two frequencies (one about one tenth the other) for times exceeding 48.
This reliably causes the graph to spiral out of control and I was wondering if there was a fix.
Paris Miles-Brenden
2008-08-16
Robert Dodier
2008-08-24
Logged In: YES
user_id=501686
Originator: NO
Thanks for this report. Please post an example in the form of specific Maxima expressions which we can use to reproduce the error. Thanks for your help. Robert Dodier
Paris Miles-Brenden
2008-08-29
Logged In: YES
user_id=2183455
Originator: YES
jacobi_cn(t,.7)*jacobi_cn(.1*t,.9)
wxplot2d([%], [t,0,40])
Raymond Toy
2008-08-30
Logged In: YES
user_id=28849
Originator: NO
This is a bug in the computation of jacobi_cn. Try, eg, jacobi_cn(100,.7)
The algorithms for computing jacobi_cn and jacobi_dn need to be examined. jacobi_sn appears to be ok.
Changing category to floating point instead of plotting.
Raymond Toy
2008-08-30
Raymond Toy
2009-01-01
This should be fixed in CVS, ellipt.lisp, rev 1.30.
Raymond Toy
2009-01-01