MaplePrimes - answers and comments on Question, find max of values, save solve values

Not sure how to save to a file

Willem Ottevanger — Mon, 18 Mar 2013 11:22:31 Z

but you dont necessarily need to save the data to a file (see attached file)

Kitonum — Mon, 18 Mar 2013 11:34:00 Z

Now all works:

restart;

with(DEtools): with(plots):

k := 0; A := 0.37e-1; B := 0.2e-6;

ode1 := diff(U(t), t) = -(A+r(t)+B*U(t))*U(t);

ic[1] := U(365*k) = 1000;

R := RandomTools:-Generate(distribution(Normal(-0.2e-1, 0.4e-1)), makeproc = true);

r := proc (t) if not type(t, numeric) then return 'procname(args)' end if; R() end proc;

for i to 5 do

sol[i] := dsolve({ode1, ic[1]}, numeric, known = r(t), method = classical[foreuler], stepsize = 0.1e-1):

p[i] := odeplot(sol[i], [[t, U(t)]], t = 365*k .. 365*k+365, colour = red) end do:

display(seq(p[i], i = 1 .. 5));

for i to 5 do x[i] := proc(t) options operator, arrow; rhs(sol[i](t)[2]) end proc:

y[i] := x[i](2) end do:

seq(y[i], i=1..5);

max(seq(y[i], i=1..5));

reemeaaaah — Mon, 18 Mar 2013 18:51:53 Z

Great!

Thank you so much for your replay and the attached file.

reemeaaaah — Mon, 18 Mar 2013 18:55:30 Z

Yes it works!

Thank you so much for your correction and replay.