=================================================================== RCS file: /home/cvs/OpenXM/src/asir-contrib/packages/doc/Attic/Diff.texi,v retrieving revision 1.5 retrieving revision 1.6 diff -u -p -r1.5 -r1.6 --- OpenXM/src/asir-contrib/packages/doc/Attic/Diff.texi 2000/01/31 11:01:00 1.5 +++ OpenXM/src/asir-contrib/packages/doc/Attic/Diff.texi 2000/02/20 11:30:44 1.6 @@ -1,10 +1,11 @@ -@c $OpenXM: OpenXM/src/asir-contrib/packages/doc/Diff.texi,v 1.4 2000/01/21 12:54:39 okutani Exp $ +@c $OpenXM: OpenXM/src/asir-contrib/packages/doc/Diff.texi,v 1.5 2000/01/31 11:01:00 okutani Exp $ @node Differential equations,,, その他の函数 @section Differential equations ファイル @file{gr}, @file{Matrix} が必要です. @tex -この節で紹介される関数では微分作用素は次のようなリストで表現します. +この節で紹介される関数では微分作用素はリストまたは多項式で表現されます. +リストによる表現は次のようになります. $$ [ [f_{\alpha},[\alpha_{1},\ldots,\alpha_{n}]],\ldots ] $$ これは $$ \sum_{\alpha}f_{\alpha}\partial^{\alpha} $$ @@ -12,11 +13,15 @@ $$ \sum_{\alpha}f_{\alpha}\partial^{\alpha} $$ $$ (\sum_{\alpha^{(i)}}f_{\alpha^{(i)}}\partial^{\alpha^{(i)}})\bullet u = 0 \quad (i = 1,\ldots,s) $$ などのように複数の微分作用素を表現するときは微分作用素のリストを使います. $$ [ [ [f_{\alpha^{(1)}},[\alpha_{1}^{(1)},\ldots,\alpha_{n}^{(1)}]],\ldots ],\ldots,[ [f_{\alpha^{(s)}},[\alpha_{1}^{(s)},\ldots,\alpha_{n}^{(s)}]],\ldots ] ] $$ -例えば微分作用素$x dx + y dy + 1$の場合は +例えば微分作用素$x \partial_{x} + y \partial_{y} + 1$の場合は $$ [ [x,[1,0]],[y,[0,1]],[1,[0,0]] ] $$ -となります. また微分作用素のリストで$x dx + y dy + 1, dx^2 + dy^2$を表すと +となります. また微分作用素のリストで$x \partial_{x} + y \partial_{y} + 1, {\partial_{x}}^{2} + {\partial_{y}}^{2}$を表すと $$ [ [ [x,[1,0]],[y,[0,1]],[1,[0,0]] ],[ [1,[2,0]],[1,[0,2]] ] ] $$ -となります. またこれらの表現法を使うときは変数リストを常に意識している必要があります。 +となります. またこの表現法を使うときは変数リストを常に意識している必要があります. +次に多項式による表現について述べます. 変数$x$に対する微分は$dx$で表現されます. +例えば$x \partial_{x} + y \partial_{y} + 1$については +$$ x*dx+y*dy+1 $$ +と表現されます. @end tex @menu @c * diff_op_hg1:: @@ -26,6 +31,9 @@ $$ [ [ [x,[1,0]],[y,[0,1]],[1,[0,0]] ],[ [1,[2,0]],[1, * diff_op_appell4:: @c * diff_op_selberg2:: @c * diff_op_gkz:: +* diff_op_tosm1:: +* diff_op_toasir:: +* diff_op_fromasir:: * diff_act:: @c * diff_act_hg1:: @c * diff_act_appell1:: @@ -72,6 +80,97 @@ $$ [ [ [x,[1,0]],[y,[0,1]],[1,[0,0]] ],[ [1,[2,0]],[1, [(-a-b-1)*y+c2,[0,1]], [(-a-b-1)*x,[1,0]], [-b*a,[0,0]] ] ] @end example +@node diff_op_tosm1,,, Differential equations +@subsection @code{diff_op_tosm1} +@findex diff_op_tosm1 +@table @t +@item diff_op_tosm1(@var{LL},@var{V}) +:: リスト形式の微分作用素リストを sm1 形式に変換します. +@end table +@table @var +@item return +リスト +@item LL +リスト +@item V +リスト +@end table +@itemize @bullet +@item 微分作用素の係数は整数多項式に変換されます. +@item @code{diff_op_tosm1}の例 +@end itemize +@example +[299] diff_op_tosm1([[[x,[2,0]],[-1,[0,0]]], + [[y,[0,2]],[-1,[0,0]]]],[x,y]); +[ + ( + (1) x) dx^2 + ( + (-1)), + ( + (1) y) dy^2 + ( + (-1))] + +[300] diff_op_tosm1([[[x,[1,0]],[y,[0,1]],[1,[0,0]]], + [[1,[2,0]],[1,[0,2]]]],[x,y]); +[ + ( + (1) x) dx + ( + (1) y) dy + ( + (1)), + ( + (1)) dx^2 + ( + (1)) dy^2] + +[301] diff_op_tosm1([[[1/2,[1,0]],[1,[0,0]]], + [[1/3,[0,1]],[1/4,[0,0]]]],[x,y]); +[ + ( + (6)) dx + ( + (12)), + ( + (4)) dy + ( + (3))] + +[302] diff_op_tosm1([[[1/2*x,[1,0]],[1,[0,0]]], + [[1/3*y,[0,1]],[1/4,[0,0]]]],[x,y]); +[ + ( + (6) x) dx + ( + (12)), + ( + (4) y) dy + ( + (3))] +@end example + +@node diff_op_toasir,,, Differential equations +@subsection @code{diff_op_toasir} +@findex diff_op_toasir +@table @t +@item diff_op_toasir(@var{LL},@var{V}) +:: リスト形式の微分作用素リスト @var{LL} を @code{asir} の多項式に変換します. +@end table +@table @var +@item return +リスト +@item LL +リスト +@item V +リスト +@end table +@itemize @bullet +@item @code{diff_op_toasir}の例 +@end itemize +@example +[303] diff_op_toasir([[[1/2*x,[1,0]],[1,[0,0]]], + [[1/3*y,[0,1]],[1/4,[0,0]]]],[x,y]); +[1/2*x*dx+1,1/3*y*dy+1/4] + +[304] diff_op_toasir([[[x,[1,0]],[y,[0,1]],[1,[0,0]]], + [[1,[2,0]],[1,[0,2]]]],[x,y]); +[x*dx+y*dy+1,dx^2+dy^2] +@end example + +@node diff_op_fromasir,,, Differential equations +@subsection @code{diff_op_fromasir} +@findex diff_op_fromasir +@table @t +@item diff_op_fromasir(@var{D_list},@var{V}) +:: @code{asir} の多項式からリスト形式の微分作用素リストに変換します. +@end table +@table @var +@item return +リスト +@item D_list +リスト +@item V +リスト +@end table +@itemize @bullet +@item @code{diff_op_fromasir}の例 +@end itemize +@example +[305] diff_op_fromasir([1/2*x*dx+1,1/3*y*dy+1/4],[x,y]); +[[[1/2*x,[1,0]],[1,[0,0]]],[[1/3*y,[0,1]],[1/4,[0,0]]]] + +[306] diff_op_fromasir([x*dx+y*dy+1,dx^2+dy^2],[x,y]); +[[[x,[1,0]],[y,[0,1]],[1,[0,0]]],[[1,[2,0]],[1,[0,2]]]] +@end example + @node diff_act,,, Differential equations @subsection @code{diff_act} @findex diff_act @@ -83,7 +182,7 @@ $$ [ [ [x,[1,0]],[y,[0,1]],[1,[0,0]] ],[ [1,[2,0]],[1, @item return 有理式 @item L -リスト +リスト or 多項式 @item F 有理式 @item V @@ -98,6 +197,9 @@ $$ [ [ [x,[1,0]],[y,[0,1]],[1,[0,0]] ],[ [1,[2,0]],[1, [303] diff_act([[1,[1,0]],[1,[0,1]]],x^2+y^2,[x,y]); 2*x+2*y + +[349] diff_act(x*dx+y*dy, x^2+x*y+y^2, [x,y]); +2*x^2+2*y*x+2*y^2 @end example @node diff_act_appell4,,, Differential equations @@ -157,6 +259,9 @@ $$ [ [ [x,[1,0]],[y,[0,1]],[1,[0,0]] ],[ [1,[2,0]],[1, [298] diff_poly_solve([[[x,[1,0]],[-2,[0,0]]],[[y,[0,1]],[-2,[0,0]]]],5,[x,y]); [_33*y^2*x^2,[_33]] + +[356] diff_poly_solve([x*dx+y*dy-3,dx+dy],4,[x,y]); +[-_126*x^3+3*_126*y*x^2-3*_126*y^2*x+_126*y^3,[_126]] @end example @node diff_poly_solve_hg1,,, Differential equations @@ -236,6 +341,9 @@ $$ [ [ [x,[1,0]],[y,[0,1]],[1,[0,0]] ],[ [1,[2,0]],[1, @example [333] diff_rat_solve([[[x,[1]],[1,[0]]]],x,1,[x]); [(_8)/(x),[_8]] + +[361] diff_rat_solve([x*(1-x)*dx^2+(1-3*x)*dx-1],1-x,2,[x]); +[(_180)/(-x+1),[_180]] [350] D = diff_op_appell4(0,0,3,0,[x,y])$ [351] diff_rat_solve(D,x^2,2,[x,y]);