MathGroup Archive 2003

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

Search the Archive

Re: Can someone tell me why NDsolve isn't working here?

  • To: mathgroup at smc.vnet.net
  • Subject: [mg44084] Re: Can someone tell me why NDsolve isn't working here?
  • From: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>
  • Date: Wed, 22 Oct 2003 03:24:35 -0400 (EDT)
  • Organization: Universitaet Leipzig
  • References: <bmqplo$rhb$1@smc.vnet.net>
  • Reply-to: kuska at informatik.uni-leipzig.de
  • Sender: owner-wri-mathgroup at wolfram.com

Hi,

 Ry[t] + Fy[t] == 0

is a algebraic equation and only Mathematica 5.0 can solve
mixed diferential/algebraic equations.

Regards
  Jens

Mizrandir wrote:
> 
> Hello, Im trying to solve a system of 6 differential equations with
> NDSolve but I can't get Mathematica to do it. I wrote:
> 
> R = 1;
> M = 1;
> L = 1;
> NDSolve[{Rx[t] + Fx[t] == 2M x''[t], Ry[t] + Fy[t] == 0,
>     Fx[t] R == -M R x''[t], -Rx[t] ==
>       M x''[t] + (M L \[Theta]''[t] Cos[-\[Theta][t]])/
>           2 + (M L \[Theta]'[t]^2 Sin[-\[Theta][t]])/2, -Ry[
>           t] == (M L \[Theta]''[t] Sin[-\[Theta][t]])/
>           2 - (M L \[Theta]'[t]^2 Cos[-\[Theta][t]])/
>           2, -(L Cos[\[Theta][t]] Rx[t])/2 + (L Sin[\[Theta][t]])/
>           2 == (M L^2 \[Theta]''[t])/12, x[0] == 0, \[Theta][0] ==
> Pi/2,
>     Ry[0] == -M/2, Rx[0] == 0, Fy[0] == 5M/2, Fx[0] == 0}, {x,
> \[Theta], Ry,
>     Rx, Fy, Fx}, {t, 0, 1}]
> 
> And Mathematica gave me the following error:
> 
> NDSolve::"ndode": "Input is not an ordinary differential equation."
> 
> Does anybody know what am I doing wrong. Can I use Mathematica to
> solve this?
> 
> Thanks in advance, miz.


  • Prev by Date: Help understanding Accuracy[]
  • Next by Date: [no subject]
  • Previous by thread: Re: Can someone tell me why NDsolve isn't working here?
  • Next by thread: ask for help on list operation