Advanced Search

Question: ODEs with pertuabutions at specific time points.

Question:ODEs with pertuabutions at specific time points.

Colin Gillespie 60 Maple
example
December 19 2007
0
Dear All, I have two odes, dX/dt = -X(t)*Y(t) and dY(t)/dt = X(t)*Y(t) I can solve these numerically using the following code: eq_1:=diff(X(t), t)=-X(t)*Y(t); eq_2:=diff(Y(t),t)=X(t)*Y(t); IC:={eq_1, eq_2,X(0)=10,Y(0)=1}; F:=dsolve(IC, {X(t),Y(t)},type=numeric, method=rosenbrock); h := theta->eval(X(t),F(theta)); h(1); However, I would like to add perturbations at specific time points. For example, at t=10, set X=5, or at t=15, Y=0. How should I do this in maple? Note: my actual situation involves a 100 odes. Thanks, Colin
Please Wait...