version 1.39, 2010/07/16 07:55:06 |
version 1.41, 2010/09/29 00:33:23 |
|
|
%% $OpenXM: OpenXM/src/asir-doc/exp/exp-ja.texi,v 1.38 2010/02/10 05:49:58 takayama Exp $ |
%% $OpenXM: OpenXM/src/asir-doc/exp/exp-ja.texi,v 1.40 2010/09/28 07:24:34 takayama Exp $ |
\input texinfo |
\input texinfo |
@iftex |
@iftex |
@catcode`@#=6 |
@catcode`@#=6 |
|
|
@title $B<B83E*;EMM$N4X?t(B |
@title $B<B83E*;EMM$N4X?t(B |
@subtitle Risa/Asir $B<B83E*;EMM4X?t@bL@=q(B |
@subtitle Risa/Asir $B<B83E*;EMM4X?t@bL@=q(B |
@subtitle 1.0 $BHG(B |
@subtitle 1.0 $BHG(B |
@subtitle 2010 $BG/(B 2 $B7n(B |
@subtitle 2010 $BG/(B 9 $B7n(B |
|
|
@author by Risa/Asir committers |
@author by Risa/Asir committers |
@page |
@page |
Line 134 ChangeLog $B$N9`L\$O(B www.openxm.org $B$N(B cvswe |
|
Line 134 ChangeLog $B$N9`L\$O(B www.openxm.org $B$N(B cvswe |
|
* dcurrenttime:: |
* dcurrenttime:: |
* nd_gr nd_gr_trace nd_weyl_gr nd_weyl_gr_trace (module):: |
* nd_gr nd_gr_trace nd_weyl_gr nd_weyl_gr_trace (module):: |
* nd_gr nd_gr_trace nd_weyl_gr nd_weyl_gr_trace (option):: |
* nd_gr nd_gr_trace nd_weyl_gr nd_weyl_gr_trace (option):: |
* tk_pf.rkn:: |
* tk_pfn.rkn:: |
|
* tk_pfn.graph:: |
|
* tk_rk.runge_kutta_4:: |
|
* tk_rk.runge_kutta_4_linear:: |
@end menu |
@end menu |
|
|
@node $B%/%*!<%H(B,,, $B<B83E*;EMM$N4X?t(B |
@node $B%/%*!<%H(B,,, $B<B83E*;EMM$N4X?t(B |
|
|
@item @code{OpenXM/src/asir-contrib/packages/src/tk_pfn.rr} 1.1, 1.2 |
@item @code{OpenXM/src/asir-contrib/packages/src/tk_pfn.rr} 1.1, 1.2 |
@item $B$3$N%b%8%e!<%k$NA0?H$O(B @code{tk_pf2.rr} $B$G$"$k(B. $B$3$l$OFHN)JQ?t$,(B 2 $B8D$N>l9g$G$"$k(B. |
@item $B$3$N%b%8%e!<%k$NA0?H$O(B @code{tk_pf2.rr} $B$G$"$k(B. $B$3$l$OFHN)JQ?t$,(B 2 $B8D$N>l9g$G$"$k(B. |
@end itemize |
@end itemize |
|
|
|
|
|
@node tk_pfn.graph,,, $B<B83E*;EMM$N4X?t(B |
|
@subsection @code{tk_pfn.graph} |
|
@findex tk_pfn.graph |
|
|
|
@comment --- $B4X?t$N4JC1$J@bL@(B --- |
|
@comment --- @itemx $B$OJ#?t$KBP$7$F@bL@$r0l$D$D$1$k>l9g$K;H$&(B --- |
|
@table @t |
|
@item tk_pfn.graph(@var{Pf},@var{Dom},@var{Iv},@var{Step}) |
|
:: 2$BJQ?t(B Pfaffian $BJ}Dx<0$r(B Runge-Kutta $BK!$G2r$$$F%0%i%UI=<($9$k(B. |
|
@end table |
|
|
|
@table @var |
|
@item return |
|
$B%j%9%H(B $B%j%9%H$NMWAG$O0J2<$N7A<0(B [[x$B$NCM(B,y$B$NCM(B],Y_1$B$NCM(B,Y_2$B$NCM(B, ...]. |
|
[x$B$NCM(B, y$B$NCM(B] $B$O(B [0,0],[0.2,0],[0.4,0], ... , [0,0.2],[0.2,0.2], ... $B$N$h$&$K(B |
|
y $B$NCM$,30B&%k!<%W(B, x $B$NCM$,FbB&%k!<%W$N7A<0$GA}$($k(B. |
|
@item Pf, Dom, Iv, Step |
|
@var{Pf} $B$O(B Pfaffian $BJ}Dx<0$N78?t9TNs%j%9%H(B. $BFHN)JQ?t$O(B x, y $B$G8GDj(B. |
|
@var{Dom} $B%j%9%H(B. $B2r$/$Y$-NN0h(B. |
|
@var{Iv} $B%j%9%H(B. $BNN0h$N:8C<$G$N=i4|CM(B. |
|
@var{Step} $B9o$_I}(B. |
|
@end table |
|
|
|
@itemize @bullet |
|
@item tk_pf2.rr, mt_graph.rr $B$r(B import $B$7$F$*$/I,MW$,$"$k(B. |
|
@item $B$3$N4X?t$OO"N)(BPfaffian$BJ}Dx<0(B dY/dx = Pf[0] Y, dY/dy = Pf[1] Y $B$r?tCME*$K2r$$$F%0%i%UI=<($9$k(B. |
|
@item @var{Dom} $B$O(B [[xmin,xmax],[ymin,ymax]] $B$N7A<0(B. |
|
@item $BNc$O%=!<%9%3!<%I(B (@code{OpenXM/lib/asir-contrib/tk_pfn.rr} )$B$N(B @code{tk_pfn.testgraph1()}, @code{tk_pfn.testgraph2()} $B$r;2>H(B. |
|
@item option $B$H$7$F$O(B fit=1 $B$,$"$k(B. Z$B<4$rE,59D4@0$9$k(B. |
|
@item @var{Dom} $B$NC<$O%0%i%UI=<($N;~$K0lIt%+%C%H$5$l$k$N$GCm0U(B. |
|
@end itemize |
|
|
|
@example |
|
[1355] import("tk_pf2.rr"); import("mt_graph.rr"); import("tk_pfn.rr"); |
|
[1590] tk_pfn.testgraph1(); |
|
|
|
$B$3$3$G(B testgraph1() $B$O0J2<$N$H$*$j(B. |
|
def testgraph1() { |
|
/* tk_bess2.bess2pf(1/2); */ |
|
Pf= [[[ 0, (1)/(x), 0 ], |
|
[ -x, (2*x^2+1)/(x), -2*x ], |
|
[ -y, 0, 0 ]], |
|
[[ 0, 0, (1)/(y) ], |
|
[ -x, 0, 0 ], |
|
[ -x, (1/2)/(x), (-1/2)/(y) ]]]; |
|
/* tk_bess2.bess2Iv(1/2,[0.5,1.5]); */ |
|
Iv = [0.105994,-0.651603,-0.760628]; |
|
Dom=[[0.5,1.5],[1.5,9]]; |
|
Step = 0.5; |
|
return tk_pfn.graph(Pf,Dom,Iv,Step | fit=1); |
|
} |
|
|
|
|
|
@end example |
|
|
|
|
|
@comment --- ChangeLog $B$r=q$/(B. $BF05!(B. $B%=!<%9%3!<%I$N0LCV(B. $BJQ99F|;~(B $B$J$I(B CVS$B%5!<%P$r8+$k$?$a(B |
|
@comment --- openxm $B$N30It$+$i$N4sM?$b=R$Y$k(B. Credit. |
|
@noindent |
|
ChangeLog |
|
@itemize @bullet |
|
@item $B$3$N4X?t$O(B 2010-08 $B$K:G=i$NHG$,=q$+$l$?(B. |
|
@item @code{OpenXM/src/asir-contrib/packages/src/tk_pfn.rr} 1.8 |
|
@end itemize |
|
|
|
@node tk_rk.runge_kutta_4,,, $B<B83E*;EMM$N4X?t(B |
|
@subsection @code{tk_rk.runge_kutta_4} |
|
@findex tk_rk.runge_kutta_4 |
|
|
|
@comment --- $B4X?t$N4JC1$J@bL@(B --- |
|
@comment --- @itemx $B$OJ#?t$KBP$7$F@bL@$r0l$D$D$1$k>l9g$K;H$&(B --- |
|
@table @t |
|
@item tk_rk.runge_kutta_4(@var{Eq},@var{X},@var{Y},@var{X0},@var{Y0},@var{Terminal},@var{Step}) |
|
:: 4$B<!$N(B Runge-Kutta $BK!$K$h$kHyJ,J}Dx<0$N?tCM6a;w2r(B |
|
@end table |
|
|
|
@table @var |
|
@item return |
|
$B%j%9%H(B $B%j%9%H$NMWAG$O0J2<$N7A<0(B [@var{X}$B$NCM(B,Y_1$B$NCM(B,Y_2$B$NCM(B, ...]. |
|
@var{X} $B$NCM$O8:$C$F$$$/(B. $B$h$C$F%j%9%H$N@hF,$,(B @var{Terminal} $BIU6a$G$N(B @var{Y} $B$NCM(B. |
|
@item Eq, X, Y, Step |
|
@var{Eq} $B$O(B $BJ}Dx<0$N1&JU(B. $B%j%9%H(B. Y[0]'=Eq[0], Y[1]'=Eq[1], ... $B$G$"$k(B. |
|
@var{X} $BFHN)JQ?tL>(B. |
|
@var{Y} $B%j%9%H(B. $B=>B0JQ?t$N%j%9%H(B. |
|
@var{Step} $B9o$_I}(B. |
|
@item X0, Y0, Terminal |
|
@var{X0} $B=PH/E@$N(B @var{X} $B$NCM(B. |
|
@var{Y0} $B=PH/E@$G$N(B @var{Y} $B$N=i4|CM(B. |
|
@var{Terminal} @var{X} $B$N=*CeE@(B. |
|
@end table |
|
|
|
@itemize @bullet |
|
@item taka_runge_kutta.rr $B$r(B import $B$7$F$*$/I,MW$,$"$k(B. |
|
@item $B$3$N4X?t$OO"N)>oHyJ,J}Dx<0(B Y[0]'=Eq[0], Y[1]'=Eq[1], ... $B$r?tCME*$K2r$/(B. |
|
@item $BNc$O%=!<%9%3!<%I(B (@code{OpenXM/lib/asir-contrib/src/taka_runge_kutta.rr} )$B$N(B |
|
@code{tk_rk.test4()} $B$r;2>H(B. |
|
@end itemize |
|
|
|
@example |
|
[1355] import("taka_runge_kutta.rr"); |
|
[1590] tk_rk.test4(); |
|
|
|
$B$3$3$G(B test4() $B$O0J2<$N$H$*$j(B. $B?6F0$NJ}Dx<0(B, y0'=y1, y1'=-y0 (y0''+y0=0). $BEz$O(B y0=cos(x) |
|
taka_plot_auto $B$O2<J}8~$G(B y $B$,@5(B. |
|
|
|
def test4() { |
|
A=runge_kutta_4([y1,-y0],x,[y0,y1],0,[1,0],3.14*2,0.1); |
|
taka_plot_auto(A); |
|
return(A); |
|
} |
|
|
|
@end example |
|
|
|
|
|
@comment --- ChangeLog $B$r=q$/(B. $BF05!(B. $B%=!<%9%3!<%I$N0LCV(B. $BJQ99F|;~(B $B$J$I(B CVS$B%5!<%P$r8+$k$?$a(B |
|
@comment --- openxm $B$N30It$+$i$N4sM?$b=R$Y$k(B. Credit. |
|
@noindent |
|
ChangeLog |
|
@itemize @bullet |
|
@item $B$3$N4X?t$O(B 2000 $BG/Be$NA0H>$K:G=i$NHG$,=q$+$l$?(B. 2010$BG/(B Pfaffian $B$N?tCM2r@O$N0Y$K:FEY@0Hw(B |
|
@item @code{OpenXM/src/asir-contrib/packages/src/taka_runge_kutta.rr} 1.17 |
|
@end itemize |
|
|
|
@node tk_rk.runge_kutta_4_linear,,, $B<B83E*;EMM$N4X?t(B |
|
@subsection @code{tk_rk.runge_kutta_4_linear} |
|
@findex tk_rk.runge_kutta_4_linear |
|
|
|
@comment --- $B4X?t$N4JC1$J@bL@(B --- |
|
@comment --- @itemx $B$OJ#?t$KBP$7$F@bL@$r0l$D$D$1$k>l9g$K;H$&(B --- |
|
@table @t |
|
@item tk_rk.runge_kutta_4_linear(@var{P},@var{X},@var{Y},@var{X0},@var{Y0},@var{Terminal},@var{Step}) |
|
:: 4$B<!$N(B Runge-Kutta $BK!$K$h$kHyJ,J}Dx<0$N?tCM6a;w2r(B. $B@~7AJ}Dx<0@lMQ(B. |
|
@end table |
|
|
|
@table @var |
|
@item return |
|
$B%j%9%H(B $B%j%9%H$NMWAG$O0J2<$N7A<0(B [@var{X}$B$NCM(B,Y_1$B$NCM(B,Y_2$B$NCM(B, ...]. |
|
@var{X} $B$NCM$O8:$C$F$$$/(B. $B$h$C$F%j%9%H$N@hF,$,(B @var{Terminal} $BIU6a$G$N(B @var{Y} $B$NCM(B. |
|
@item P, X, Y, Step |
|
@var{P} $B$O(B $BJ}Dx<0$N1&JU(B. $B%j%9%H(B. Y'=P Y $B$G$"$k(B. $B=>B0JQ?t(B @var{Y} $B$OITMW(B. |
|
@var{X} $BFHN)JQ?tL>(B. |
|
@var{Y} $B%j%9%H(B. $B=>B0JQ?t$N%j%9%H(B. $B=>B0JQ?t$O<+F0@8@.$5$l$k(B. $B;H$o$l$F$$$J$$$N$G6u%j%9%H$G$h$$(B. |
|
@var{Step} $B9o$_I}(B. |
|
@item X0, Y0, Terminal |
|
@var{X0} $B=PH/E@$N(B @var{X} $B$NCM(B. |
|
@var{Y0} $B=PH/E@$G$N(B @var{Y} $B$N=i4|CM(B. |
|
@var{Terminal} @var{X} $B$N=*CeE@(B. |
|
@end table |
|
|
|
@itemize @bullet |
|
@item taka_runge_kutta.rr $B$r(B import $B$7$F$*$/I,MW$,$"$k(B. |
|
@end itemize |
|
|
|
@example |
|
[1355] import("taka_runge_kutta.rr"); |
|
[1590] A=tk_rk.runge_kutta_4_linear([[0,1],[-1,0]],x,[ ], 0, [1,0], 3.14*2, 0.1); |
|
[1591] taka_plot_auto(A); |
|
|
|
$B?6F0$NJ}Dx<0(B, y0'=y1, y1'=-y0 (y0''+y0=0). $BEz$O(B y0=cos(x) $B$r2r$$$F$$$k(B. |
|
taka_plot_auto $B$O2<J}8~$G(B y $B$,@5(B. |
|
|
|
@end example |
|
|
|
|
|
@comment --- ChangeLog $B$r=q$/(B. $BF05!(B. $B%=!<%9%3!<%I$N0LCV(B. $BJQ99F|;~(B $B$J$I(B CVS$B%5!<%P$r8+$k$?$a(B |
|
@comment --- openxm $B$N30It$+$i$N4sM?$b=R$Y$k(B. Credit. |
|
@noindent |
|
ChangeLog |
|
@itemize @bullet |
|
@item 2010$BG/(B Pfaffian $B$N?tCM2r@O$N0Y$K:F@0Hw(B. |
|
@item @code{OpenXM/src/asir-contrib/packages/src/taka_runge_kutta.rr} 1.17 |
|
@end itemize |
|
|
|
|
|
|
@node fj_simp.simplify,,, $B<B83E*;EMM$N4X?t(B |
@node fj_simp.simplify,,, $B<B83E*;EMM$N4X?t(B |
@subsection @code{fj_simp.simplify} |
@subsection @code{fj_simp.simplify} |