Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Mathematica: NDSolve in 3 dimensions

  1. Apr 28, 2012 #1
    Hi

    I am trying to solve Newtons equation for a particle in (x, y, z) using NDsolve. Here I what I have so far:

    Code (Text):

    sol = NDSolve[{
        x''[t] == acceleration[x'[t], y[t], z[t]],
        y''[t] == acceleration[y'[t], x[t], z[t]],
        z''[t] == 0,
        x[0] == 2, x'[0] == 0, y[0] == 0, y'[0] == 0, z[0] == 0,
        z'[0] == 250},
       {x[t], y[t], z[t]}, {t, 0, 1}];

    ParametricPlot[{x[t], x'[t]} /. sol, {t, 0, 0.001}, PlotRange -> {{0, 0.35}, {0, 5}}, AspectRatio -> 0.75]
     
    The force acting on the particle only affects (x, y), not z. However, my plot is empty.

    Is my use of NDsolve correct?

    Best,
    Niles.
     
  2. jcsd
  3. Apr 28, 2012 #2
    I found one error. It should be {x, y, z}, not {x[t], y[t], z[t]}. However, the parametric plot still doesn't work.
     
  4. Apr 28, 2012 #3
    It works now, there was a problem with my IV. Just had to pick a correct value.

    Best,
    Niles.
     
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook