[Date Index]
[Thread Index]
[Author Index]
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**
| |