Subject: about rk when introducing a function parameter.
From: Xiang Liu
Date: Tue, 12 Aug 2008 18:05:57 +0800
load(dynamics);
k(x):= if oddp(floor(x/3)) then 1 else 0;
eq:[2.3-k, -x +k];
sol:rk(eq, [x,y],[3,5],[t, 0, 10, 0.3]);
sol[2];
----
I got the result(below), why k was not be taken into computing? I just
began my maxima journey yesterday. Say thanks in advance.
[0.3,0.05*(-2*k+4*(2.3-k)+4.6)+3,0.05*(4*(k-0.15*(2.3-k)-3)+2*k-0.3*(2.3-k)-6)+5]