dd6ea8c50b5b52a55f8876dd04f246494e687bb2
1 <?xml version="1.0" encoding="UTF-8"?>
3     <refnamediv>
4         <refname>%sn</refname>
5         <refpurpose>ヤコビ楕円関数</refpurpose>
6     </refnamediv>
7     <refsynopsisdiv>
8         <title>呼び出し手順</title>
9         <synopsis>[y]=%sn(x,m)</synopsis>
10     </refsynopsisdiv>
11     <refsection>
12         <title>パラメータ</title>
13         <variablelist>
14             <varlistentry>
15                 <term>x</term>
16                 <listitem>
17                     <para>
18                         楕円積分により定義される基本矩形の中の点;
19                         <literal>x</literal> は複素数ベクトル
20                     </para>
21                 </listitem>
22             </varlistentry>
23             <varlistentry>
24                 <term>m</term>
25                 <listitem>
26                     <para>
27                         楕円積分のパラメータ (<literal>0&lt;m&lt;1</literal>)
28                     </para>
29                 </listitem>
30             </varlistentry>
31             <varlistentry>
32                 <term>y</term>
33                 <listitem>
34                     <para>結果</para>
35                 </listitem>
36             </varlistentry>
37         </variablelist>
38     </refsection>
39     <refsection>
40         <title>説明</title>
41         <para>
42             パラメータ<literal>m</literal>を指定した
43             ヤコビの楕円関数:
44             パラメータ<literal>m</literal>の楕円積分の逆.
45         </para>
46         <para>
47             ゲインamはFortran形式で計算され,
48             楕円関数用の加算式が適用されます.
49         </para>
50     </refsection>
51     <refsection>
52         <title>例</title>
53         <programlisting role="example"><![CDATA[
54 m=0.36;
55 K=%k(m);
56 P=4*K; //実周期
57 real_val=0:(P/50):P;
58 plot(real_val,real(%sn(real_val,m)))
59 clf();
60 KK=%k(1-m);
61 Ip=2*KK;
62 ima_val1=0:(Ip/50):KK-0.001;
63 ima_val2=(KK+0.05):(Ip/25):(Ip+KK);
64 z1=%sn(%i*ima_val1,m);z2=%sn(%i*ima_val2,m);
65 plot2d([ima_val1',ima_val2'],[imag(z1)',imag(z2)']);
66 xgrid(3)
67  ]]></programlisting>
68     </refsection>
69     <refsection role="see also">
70         <title>参照</title>
71         <simplelist type="inline">
72             <member>