version 1.8, 2009/02/05 08:43:59 |
version 1.9, 2009/02/10 23:51:19 |
|
|
/* $OpenXM: OpenXM/src/asir-contrib/packages/src/taka_runge_kutta.rr,v 1.7 2009/02/05 08:36:33 takayama Exp $ */ |
/* $OpenXM: OpenXM/src/asir-contrib/packages/src/taka_runge_kutta.rr,v 1.8 2009/02/05 08:43:59 takayama Exp $ */ |
/* From misc/200205/runge-kutta.rr */ |
/* From misc/200205/runge-kutta.rr */ |
|
|
/* They have not yet been registered in names.rr */ |
/* They have not yet been registered in names.rr */ |
Line 69 def taka_runge_kutta_replace(V,F,Y,N,X,Xk,Rule_vector) |
|
Line 69 def taka_runge_kutta_replace(V,F,Y,N,X,Xk,Rule_vector) |
|
for (J=0; J<N; J++) { |
for (J=0; J<N; J++) { |
V[I] = subst(V[I],Y[J],Rule_vector[J]); |
V[I] = subst(V[I],Y[J],Rule_vector[J]); |
} |
} |
|
V[I] = eval(V[I]*exp(0)); |
} |
} |
} |
} |
|
|