MaplePrimes - answers and comments on Question, dsolve system with boundary conditions

boundary condition

Robert Israel — Sun, 21 Mar 2010 22:46:24 Z

Note that using EDP[2] your IBC[2] reduces to i(100)*(9*s(100)-4) = 0. There are really two possibilities: i(100) = 0 or s(100) = 4/9. Since your IBC[1] says s(100) = 1, you can't have s(100) = 4/9, so the second boundary condition is really i(100) = 0.
And then dsolve(..., numeric) will happily solve your system, but of course the solution is s(z) = 1, i(z) = 0.

Ah ok! I've changed my

Huricane — Mon, 22 Mar 2010 21:50:33 Z

Ah ok! I've changed my boundary conditions and i got : proc(x_rkf45) ... end proc

How can i now see my fonction ?

Thanks for your help !

Thanks a lot !

Huricane — Tue, 23 Mar 2010 18:24:30 Z

Thanks a lot !

seeing the function

Robert Israel — Tue, 23 Mar 2010 03:10:30 Z

odeplot from the plots package. For example:

> Sol:= dsolve({EDP, IBC}, numeric); 
  plots[odeplot](Sol, [[z,s(z)],[z,i(z)]], z=0..100, colour=[red,blue], legend = [s(z), i(z)]);