Re: How to rectify the error for NDSolve ?
- To: mathgroup at smc.vnet.net
- Subject: [mg127148] Re: How to rectify the error for NDSolve ?
- From: Murray Eisenberg <murray at math.umass.edu>
- Date: Mon, 2 Jul 2012 05:29:06 -0400 (EDT)
- Delivered-to: l-mathgroup@mail-archive0.wolfram.com
- References: <201207010611.CAA01041@smc.vnet.net>
- Reply-to: murray at math.umass.edu
So many errors in your code!
(1) As somebody already noted in a reply to your earlier post (before
you showed any code), you forgot to use an uppercase letter in Clear.
And you also forgot to use an uppercase letter in the sine expression on
the right-hand side of the differential equation.
(2) The root cause of the error message you received is because first
you set t = 5 and then a multiple of that; now t is a constant! But then
you use the _same_ symbol t as if it were the variable in your
differential equation. Change either the constant t to some other name
or the variable t in your differential equation et seq. to some other name.
(3) You are not using correct syntax in the definition of eqn: all
function arguments should be enclosed in brackets, not parentheses. And
you forgot to include the argument of t for x at the end of the
left-hand side of the equation. Thus you seem to want:
eqn = x''[\[Tau]] +1/Q x'[\[Tau]] +
(1+Subscript[V,1] Sqrt[2] Sin[2 \[Tau]]) x[t]
==\[CapitalGamma] (Sin[\[Tau]+\[Phi]]);
Note that I removed the parentheses in the right-hand side, as there's
no reason to have them.
(4) Remove Subscript[V,p] and Subscript[V,1] from the Clear expression,
as they are not Symbols (or Strings). If they don't already have values,
you don't need to get "clear" them in any sense. If they do, you have to
use Unset:
Unset[Subscript[V,p]]
Unset[Subscript[V,1]]
(You can Unset several things at once by using Map.)
(5) After all that, at this point I do not understand your differential
equation -- I don't know what you really intend the value of \[Tau] to
be. You had set t = 5, so that after
\[Tau] = Subscript[\[Omega], 0] t;
\[Tau] is also a constant, as is, therefore, the right-hand side of your
differential equation. I suspect that's NOT what you intended.
In fact, I don't understand at all why you even set t to be a constant.
But if you simply fix all the things indicated above and omit the
assignment t = 5, I think Mathematica will not be able to deal with this
differential equation.
However, if you change the term that involving what I said should be
x[t] in place of just x to have variable \[Tau] instead of t -- is that
what you actually intended?? -- I think Mathematica will still balk at
providing a solution.
If, though, you make a change of variable from t to tt ....
neweqn = eqn /. t -> tt/10
... then the following _will_ provide a numerical solution:
NDSolve[{neweqn, x[0] == 0.5, x'[0] == 0}, x[tt], {tt, 0, 250}]
P.S. By any chance, did you retype your code when writing your e-mail
post? If so, that in itself may have introduced errors. You should use
Mathematica's menu selection Edit > Copy As Plain Text to copy, then
past into your message.
On 7/1/12 2:11 AM, Rahul Chakraborty wrote:
> Here is my code:
>
>
> clear[x, t, \[Tau], Q, Subscript[V, p], Subscript[V,
> 1], \[CapitalGamma], \[Phi]];
> Q = 100;
> \[CapitalGamma] = 50;
> Subscript[\[Omega], 0] = 10;
> Subscript[V, p] = 0.5;
> Subscript[C, 0] = 4.06 * 10^-4;
> t = 5;
> \[Tau] = Subscript[\[Omega], 0] t;
> \[Phi] = 90;
> Subscript[V, 1] = Subscript[C, 0]*Subscript[V, p];
> eqn = x'' (\[Tau]) +
> 1/Q x' (\[Tau]) + (1 +
> Subscript[V, 1] Sqrt[2]
> Sin[2 \[Tau]]) x == \[CapitalGamma] (sin[\[Tau] + \[Phi]]);
> sol = NDSolve[{eqn, x[0] == 0.5, x'[0] == 0}, x, {t, 0, 250},
> MaxSteps -> 500000, WorkingPrecision -> 25][[1]]
> ParametricPlot[Evaluate[{x[t], x'[t]} /. sol], {t, 0, 100},
> Frame -> True, PlotRange -> All,
> AxesLabel -> {"t", "\!\(\*OverscriptBox[\(\[Theta]\), \(.\)]\)"}]
>
>
> Output i'm getting is as below:
>
> NDSolve::dsvar: 5 cannot be used as a variable. >>
>
>
> {0.999855 x + Derivative[1][x]/2 + 50 x^\[Prime]\[Prime] ==
> 50 sin[140], x[0] == 0.5, Derivative[1][x][0] == 0}
>
>
> ReplaceAll::reps: {0.999855 x+x^\[Prime]/2+50 x^\[Prime]\[Prime]==50 sin[140],x[0]==0.5,(x^\[Prime])[0]==0} is neither a list of replacement rules nor a valid dispatch table, and so cannot be used for replacing. >>
>
>
> ReplaceAll::reps: {0.999855 x+0.5 x^\[Prime]+50. x^\[Prime]\[Prime]==50. sin[140.],x[0.]==0.5,(x^\[Prime])[0.]==0.} is neither a list of replacement rules nor a valid dispatch table, and so cannot be used for replacing. >>
>
> ReplaceAll::reps: {0.999855 x+0.5 x^\[Prime]+50. x^\[Prime]\[Prime]==50. sin[140.],x[0.]==0.5,(x^\[Prime])[0.]==0.} is neither a list of replacement rules nor a valid dispatch table, and so cannot be used for replacing. >>
>
> General::stop: Further output of ReplaceAll::reps will be suppressed during this calculation. >>
>
>
> Regards,
>
> rahul
>
--
Murray Eisenberg murray at math.umass.edu
Mathematics & Statistics Dept.
Lederle Graduate Research Tower phone 413 549-1020 (H)
University of Massachusetts 413 545-2859 (W)
710 North Pleasant Street fax 413 545-1801
Amherst, MA 01003-9305
- References:
- Re: How to rectify the error for NDSolve ?
- From: Rahul Chakraborty <rahul.6sept@gmail.com>
- Re: How to rectify the error for NDSolve ?