Mathematica 9 is now available
Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2003
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

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: [mg44053] Re: Can someone tell me why NDsolve isn't working here?
  • From: mizrandir at hotmail.com (Mizrandir)
  • Date: Sun, 19 Oct 2003 01:10:55 -0400 (EDT)
  • References: <bmqplo$rhb$1@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

There is an error in the equations, but after correcting it, I'm still
getting nowhere. The corrected equations are:

R = 1;
M = 1;
L = 1;
g = 9.81;
NDSolve[{Rx[t] + Fx[t] == 2M x''[t], Ry[t] + Fy[t] - 2M g == 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 g == (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 g/2, Rx[0] == 0, Fy[0] == 5M g/2, Fx[0] == 0}, {X,
\[Theta],
    Ry, Rx, Fy, Fx}, {t, 0, 1}]

miz.


  • Prev by Date: Re: two's complement function
  • Next by Date: Re: Can someone tell me why NDsolve isn't working here?
  • Previous by thread: Re: Can someone tell me why NDso
  • Next by thread: Re: Can someone tell me why NDsolve isn't working here?