The modulus is not small enough. Add one 9 to make it 0.99 and the errors will increase....
Yes, that's the problem. Using jacobi_sn(u+elliptic_kc(m),m) seems to give correct...
These bug fixes and others will appear in Maxima 5.39, which should be on the horizon....
I'm sorry for the delay. No, that was just a specification for that particular implementation...
I thought I'd check and sure enough, I forgot to specify that the q(x,y) is the nome...
R(x) is the culprit, 6th line from the bottom. If declared like this: R(x) := jacobi_cd(N...
R(x) is the culprit, 6th line from the bottom. If declared like this: R(x) := jacobi_cd(N...
Jacobi elliptic functions and complex arguments