       Re: Differential Equation: Not getting result

• To: mathgroup at smc.vnet.net
• Subject: [mg127335] Re: Differential Equation: Not getting result
• From: Murray Eisenberg <murray at math.umass.edu>
• Date: Wed, 18 Jul 2012 01:37:11 -0400 (EDT)
• Delivered-to: l-mathgroup@mail-archive0.wolfram.com
• Delivered-to: mathgroup-newout@smc.vnet.net
• Delivered-to: mathgroup-newsend@smc.vnet.net
• References: <20120717053145.BA8F96877@smc.vnet.net>

```First, it's hard to help you when you _still_ have not learned basic
Mathematica syntax. Did you, for example, look at examples for DSolve or
NDSolve to see how you denote the derivative? It should be

x'[t]

and NOT (x^')[t].

Next, there is no need for NDSolve here, as the differential equation
can be solved exactly (symbolically):

r = 1; K = 1;
eqn = x'[t] - r x[t] (1 - x[t]/K) == 0

DSolve[eqn, x[t], t]

The result, in InputForm, is:

{{x[t] -> E^t/(E^t + E^C)}}

(Did you try that? If not, why not??)

(Of course that is a "generic" solution, depending on ONE parameter
C, that does not include the trivial solution identically zero.)

Next, look at the first initial condition you want, x == 1/2.

DSolve[{eqn, x == 1/2}, x[t], t]

The result, again in InputForm, is:

{{x[t] -> E^t/(1 + E^t)}}

(You'll see a DSolve warning there about inverse functions.)

x'==0 as well. In fact, if you try it...

DSolve[{eqn, x == 1/2, x' == 0}, x[t], t]

... then you get an empty solution set!

Going back to just one initial condition, for x...

sol = First @ DSolve[{eqn, x == 1/2}, x[t], t]

... you'll get a plot from:

ParametricPlot[Evaluate[{x[t] /. sol, D[x[t] /. sol, t]}], {t, 0, 50}]

Finally, suppose you really did want, for reasons unknown, to use
NDSolve instead of solve. Still, you cannot use both initial conditions,
just the first.

On 7/17/12 1:31 AM, Rahul Chakraborty wrote:
> Dear all,
>
>   For the following code i'm not getting result. kindly tell me where i'm making mistake.
>
> Clear[x];
> r=1;
> K=1;
> eqn= (x^')[t]-r  x[t]  (1-x[t]/K)==0//Simplify;
> sol=NDSolve[{eqn,x==1/2,x'==0},x[t],{t,0,50}][]
> ParametricPlot[Evaluate[{x[t]/.sol,D[x[t]/.sol,t]}],{t,0,50},Frame->True,AxesLabel->{"x",Overscript[x,"."]},AspectRatio->1]
> ParametricPlot[Evaluate[{t,x[t]/.sol}],{t,0,50},Frame->True,AxesLabel->{"t","x"},AspectRatio->1]
>
>
> Regards,
> rc
>

--
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

```

• Prev by Date: Re: IsoWeek Function
• Next by Date: Re: Differential Equation: Not getting result
• Previous by thread: Differential Equation: Not getting result
• Next by thread: Re: Differential Equation: Not getting result