1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
|
%sn(1) Scilab Function %sn(1)
NAME
%sn - Jacobi 's elliptic function
CALLING SEQUENCE
[y]=%sn(x,m)
PARAMETERS
x : a point inside the fundamental rectangle defined by the elliptic
integral; x is a vector of complex numbers
m : parameter of the elliptic integral (0<m<1)
y : result
DESCRIPTION
Jacobi 's sn elliptic function with parameter m: the inverse of the ellip-
tic integral for the parameter m.
The amplitude am is computed in fortran and the addition formulas for
elliptic functions are applied
EXAMPLE
m=0.36;
K=%k(m);
P=4*K; //Real period
real_val=0:(P/50):P;
plot(real_val,real(%sn(real_val,m)))
xbasc();
KK=%k(1-m);
Ip=2*KK;
ima_val1=0:(Ip/50):KK-0.001;
ima_val2=(KK+0.05):(Ip/25):(Ip+KK);
z1=%sn(%i*ima_val1,m);z2=%sn(%i*ima_val2,m);
plot2d([ima_val1',ima_val2'],[imag(z1)',imag(z2)']);
xgrid(3)
SEE ALSO
%asn, %k
AUTHOR
F. D.
|
