Re: Solving simple equations
- To: mathgroup at smc.vnet.net
- Subject: [mg83206] Re: Solving simple equations
- From: Jean-Marc Gulliet <jeanmarc.gulliet at gmail.com>
- Date: Thu, 15 Nov 2007 05:32:48 -0500 (EST)
- Organization: The Open University, Milton Keynes, UK
- References: <fheh0u$mm7$1@smc.vnet.net>
Holger wrote:
> I'm trying to solve these two simple equations but it doesn't work.
>
> meqn = {
> x'[t]==beta(x[t]+(Subscript[P,f]-p[t])),
> p'[t]==(Tanh[Subscript[a,1](x[t]+(Subscript[P,f]-p[t]))+
> Subscript[a,2]x[t]]-x[t])
> };
>
> eqp = NSolve[
> meqn/.{
> p'[t]->0,x'[t]->0},{p[t], x[t]
> }
> ]
>
> I guess I'm doing a mistake somewhere. Does anyone have idea where the mistake is?
So you want to solve *numerically* the system meqn (note that I have
modified the subscripted stuff).
In[1]:= meqn = {x'[t] == beta (x[t] + (Pf - p[t])),
p'[t] == (Tanh[a1 (x[t] + (Pf - p[t])) + a2 x[t]] - x[t])};
meqn /. {p'[t] -> 0, x'[t] -> 0}
Out[2]= {0 == beta (Pf - p[t] + x[t]),
0 == Tanh[a2 x[t] + a1 (Pf - p[t] + x[t])] - x[t]}
To do so, you must provide some numeric values for *all* the parameters,
that is beta, Pf, a1, and a2, (though in the above example the value of
beta does not matter and beta can be discarded before solving).
Of course some (many?) other things might be wrong, but since you did
tell us neither what you mean by "not working" nor reported any error
messages returned by by Mathematica, it is hard to tell more (perhaps
DSolve or NDSolve might be more appropriate).
Regards,
--
Jean-Marc