Hello,
I am trying to run maple file for k=1, k=2 ....
My do-loop did not work. What I would like to do is to solve the ODE then have new inc and solve it again.
I should have U(405) then U(770) ....
> restart;
> with(DEtools); with(plots);
> lambda := 0.1; delta := .5;
tau := 40;
> for k to 3 do 365*k end do;
> ode := diff(U(t), t) = -lambda(t)*U(t)*U(t);
> inc1 := U(0) = 100;