Buenas noches.
No tengo conocimiento en wxMaxima, pero estoy tratando de resolver un sistema de ecuaciones diferenciales con este programa. Me aparece un error que no sé a qué se refiere. Por favor, les envío lo que tengo escrito y al final, el error que me aparece, para ver si tienen alguna idea que me pueda ayudar a solucionar ese error. Muchas gracias.
(%i4) B:5;
k1:0.1;
k2:0.05;
k3:0.08;
(B) 5
(k1) 0.1
(k2) 0.05
(k3) 0.08
(%i5) eq1:'-diff(C(t),t,1)+k1*B-k2*Z(t)*C(t);
(eq1) mminus('diff(C(t),t,1))-0.05*C(t)*Z(t)+0.5
(%i6) eq2:'-diff(D(t),t,1)+k2*Z(t)*C(t);
(eq2) mminus('diff(D(t),t,1))+0.05*C(t)*Z(t)
(%i7) eq3:'-diff(Z(t),t,1)-k2*Z(t)*C(t)-k3*Z(t);
(eq3) mminus('diff(Z(t),t,1))-0.05*C(t)*Z(t)-0.08*Z(t)
(%i8) eq4:'-diff(P(t),t,1)+k3*Z(t);
(eq4) mminus('diff(P(t),t,1))+0.08*Z(t)
(%i9) eqs:[eq1,eq2,eq3,eq4];
(eqs) [mminus('diff(C(t),t,1))-0.05*C(t)*Z(t)+0.5,mminus('diff(D(t),t,1))+0.05*C(t)*Z(t),mminus('diff(Z(t),t,1))-0.05*C(t)*Z(t)-0.08*Z(t),mminus('diff(P(t),t,1))+0.08*Z(t)]
(%i21) dependent_variables:'[C(t),D(t),Z(t),P(t)];
(dependent_variables) [C(t),D(t),Z(t),P(t)]
(%i11) initial_conditions:[0,0,1,0];
(initial_conditions) [0,0,1,0]
(%i12) t_range:[t,0,50,0.1];
(t_range) [t,0,50,0.1]
(%i20) data: rk(eqs,dependent_variables,initial_conditions,t_range);
rk: variable name expected; found: C(t)
-- an error. To debug this try: debugmode(true);
No tengo conocimiento en wxMaxima, pero estoy tratando de resolver un sistema de ecuaciones diferenciales con este programa. Me aparece un error que no sé a qué se refiere. Por favor, les envío lo que tengo escrito y al final, el error que me aparece, para ver si tienen alguna idea que me pueda ayudar a solucionar ese error. Muchas gracias.
(%i4) B:5;
k1:0.1;
k2:0.05;
k3:0.08;
(B) 5
(k1) 0.1
(k2) 0.05
(k3) 0.08
(%i5) eq1:'-diff(C(t),t,1)+k1*B-k2*Z(t)*C(t);
(eq1) mminus('diff(C(t),t,1))-0.05*C(t)*Z(t)+0.5
(%i6) eq2:'-diff(D(t),t,1)+k2*Z(t)*C(t);
(eq2) mminus('diff(D(t),t,1))+0.05*C(t)*Z(t)
(%i7) eq3:'-diff(Z(t),t,1)-k2*Z(t)*C(t)-k3*Z(t);
(eq3) mminus('diff(Z(t),t,1))-0.05*C(t)*Z(t)-0.08*Z(t)
(%i8) eq4:'-diff(P(t),t,1)+k3*Z(t);
(eq4) mminus('diff(P(t),t,1))+0.08*Z(t)
(%i9) eqs:[eq1,eq2,eq3,eq4];
(eqs) [mminus('diff(C(t),t,1))-0.05*C(t)*Z(t)+0.5,mminus('diff(D(t),t,1))+0.05*C(t)*Z(t),mminus('diff(Z(t),t,1))-0.05*C(t)*Z(t)-0.08*Z(t),mminus('diff(P(t),t,1))+0.08*Z(t)]
(%i21) dependent_variables:'[C(t),D(t),Z(t),P(t)];
(dependent_variables) [C(t),D(t),Z(t),P(t)]
(%i11) initial_conditions:[0,0,1,0];
(initial_conditions) [0,0,1,0]
(%i12) t_range:[t,0,50,0.1];
(t_range) [t,0,50,0.1]
(%i20) data: rk(eqs,dependent_variables,initial_conditions,t_range);
rk: variable name expected; found: C(t)
-- an error. To debug this try: debugmode(true);