MathGroup Archive 2011

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

Search the Archive

Re: Problems with DSolve

  • To: mathgroup at smc.vnet.net
  • Subject: [mg119672] Re: Problems with DSolve
  • From: DrMajorBob <btreat1 at austin.rr.com>
  • Date: Fri, 17 Jun 2011 00:07:47 -0400 (EDT)
  • References: <201106160801.EAA13245@smc.vnet.net>
  • Reply-to: drmajorbob at yahoo.com

a) This is kind of like asking  "Who's on first?"  and mathematica  
answers:  "Who's on first."

No, it's kind of like entering invalid syntax and getting it back  
unevaluated.

b) That is NOT the output returned for your input, and

c) Evaluating your input gives an error message (in a different window)  
saying:

DSolve::deqn: Equation or list of equations expected instead of -0.408  
x[t]-(0.355 +(0.645 t)/(223.06 +t)) x[t]+0.311 y[t]+0.046  
z[t]+(x^\[Prime])[t]
  in the first argument {-0.408 x[t]-(0.355 +(0.645 t)/(223.06 +t))  
x[t]+0.311 y[t]+0.046 z[t]+(x^\[Prime])[t],<<4>>,0}
. >>

d) As the error message says, you need equations, not just things you want  
to be zero (if that's what they are). This is closer to the right syntax:

sol = DSolve[{x'[t] - (0.355 + (0.645 t/(t + 223.06))) x[t] +
      0.311 y[t] + 0.046 z[t] - 0.149 x[t] - 0.259 x[t] == 0,
    y'[t] 0.149 x[t] - 0.311 y[t] == 0,
    z'[t] 0.259 x[t] - 0.046 z[t] == 0, x[0] D, y[0] 0, z[0] 0}, {x[t],
     y[t], z[t]}, t]

but the last three items in your "equation" list were x[0] D, y[0] 0, and  
z[0] 0. I can't guess what you intended them to be, so I haven't changed  
them into equations.

Make all the equations into... umm... equations, and you'll have better  
luck.

e) In your listed output (which Mathematica, I guarantee, did NOT return,  
for ANY input), note that

x = A2[t] - 0.408 x[t] + (-0.355 - (0.645 t)/(223.06 + t)) x[t] +
   0.311 y[t] + 0.046 z[t]

ALSO is not an equation, since = is Set, and == is Equal.

I suspect this is part of your actual input.

Bobby

On Thu, 16 Jun 2011 03:01:05 -0500, David King <drdavidbking at gmail.com>  
wrote:

> Hello Mathematica Forum,
>
> I am a new forum contributer, so it is nice to meet you.  I was hoping
> at least one of you could help me out.  I am attempting to find an
> exact solution to a system of differential equations
>
> I  type in the following Dsolve  function:
>
>  sol=DSolve[{x'[t] -(0.355+(0.645t/(t+223.06)))x[t]+0.311y[t]
> +0.046z[t] - 0.149x[t] -0.259x[t],y'[t] 0.149x[t] -
> 0.311y[t],z'[t] 0.259x[t] -0.046z[t],x[0] D,y[0] 0,z[0] 0},
> {x[t],y[t],z[t]},t]
>
> So far, Mathematica  returns the following:
>
> Out[7]:=  DSolve[{x=A2[t] -0.408 x[t]+(-0.355-(0.645 t)/(223.06 +t)) x[t]
> +0.311 y[t]+0.046 z[t],y=A2[t] 0.149 x[t]-0.311 y[t],z=A2[t] 0.259
> x[t]-0.046 z[t],x[0] D,y[0] 0,z[0] 0},{x[t],y[t],z[t]},t]
>
> Which is what my question is.... This is kind of like asking  "Who's on
> first?"  and mathematica answers:  "Who's on first."  Which is what my
> question is....  I am pretty sure that my syntax is correct.  Anyway,
> I was hoping that you might be able to knock this one out of the
> park.
>
> I have a numerical solution, just not an exact solution.
>


-- 
DrMajorBob at yahoo.com


  • Prev by Date: Re: Sorting into bins
  • Next by Date: Re: Nudge? (Coarse and Fine Manipulation)
  • Previous by thread: Problems with DSolve
  • Next by thread: Re: Problems with DSolve