MathGroup Archive 2007

[Date Index] [Thread Index] [Author Index]

Search the Archive

NDSolve and Systems of DEs

  • To: mathgroup at smc.vnet.net
  • Subject: [mg74177] NDSolve and Systems of DEs
  • From: "e-touch\ at libero\.it" <e-touch at libero.it>
  • Date: Tue, 13 Mar 2007 04:03:54 -0500 (EST)

Hello everybody!
In order to better understand the behaviour of NDSolve, I decided to try to=
 solve a more simple system of equations...take a look here

eq1--d/dt p[x,t]==== k dd/dx2 p[x,t]
eq2--d/dt n[x,t]==== k dd/dx2 n[x,t]
eq3--d/dx El[x,t]====p[x,t]-n[x,t]
with:
p[x,0]====10*UnitStep[-x]; p[100,t]====0; p[-100,t]====10;
n[x,0]====10*UnitStep[x]; n[100,t]====10; n[-100,t]====0;
El[x,0]====0;

Eq1 and Eq2 are simple diffusion equations and El[x,t] it's just a function=
 that I tried to couple with the other two...
All I get are error messages...
"NDSolve::pdord: Some of the functions have zero differential order so the =
\
equations will be solved as a system of differential-algebraic equations."
"NDSolve::bcart: Warning: An insufficient number of boundary conditions hav=
e \
been specified for the direction of independent variable x. Artificial \
boundary effects may be present in the solution."
"NDSolve::icfail: Unable to find initial conditions which satisfy the resid=
ual \
function within specified tolerances.  Try giving initial conditions for bo=
th \
values and derivatives of the functions."

The strange thing is that if I derive the funcion El[x,t] respect to time i=
nstead of x...so with
eq3--d/dt El[x,t]====p[x,t]-n[x,t];
El[x,0]====0;
...it works!!!
Could anyone explain this to me???
Thanks everybody
Matteo



------------------------------------------------------
Leggi GRATIS le tue mail con il telefonino i-mode=99 di Wind
http://i-mode.wind.it13mar07



  • Prev by Date: Re: Problem with Which
  • Next by Date: Re: Integrate
  • Previous by thread: Re: Replacement question
  • Next by thread: Re: NDSolve and Systems of DEs