blob: ee229ebee16efecd4670d9d85ddea573e541eb93 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
// Demo function to test %sn and %k functions
// %sn is the Jacobi's elliptic function
// %k is the Jacobi's complete integral
function main()
u=[1+2*%i 2+3*%i 3+1*%i]; // the corresponding complex vector
k=0.6; // parameter of elliptic integral
[out]=%sn(u,k);
disp(out);
m=[0.1 0.2 0.3]; // parametrt of Jacobi complete integral
s=%k(m);
disp(s);
endfunction
//Output
//for %sn
// 1.5641568 - 0.2037920i 1.0580372 + 0.0100213i 1.1386855 - 0.4631921i
// for %k
// 1.6124413 1.6596236 1.7138894
|
