% error: [sn,cn,dn] = ellipj(u,m[,tol]) % % Compute the Jacobi elliptic functions sn(u|m), cn(u|m) and dn(u|m) % for complex argument u and parameter 0 <= m <= 1. % % WARNING: the approximation blows up for abs(u)>20 near m=1. % % tol is accepted for compatibility, but ignored % % Ref: Abramowitz, Milton and Stegun, Irene A % Handbook of Mathematical Functions, Dover, 1965 % Chapter 16 (Sections 16.4, 16.13 and 16.15) % % Example % m = linspace(0,1,200); u=linspace(-10,10,200); % [U,M] = meshgrid(u,m); % [sn, cn, dn] = ellipj(U,M); % imagesc(sn); % % See also: ellipke