Mathematica: NDSolve in 3 dimensions

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
2 replies · 3K views
Niles
Messages
1,834
Reaction score
0
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:
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?


Niles.
 
on Phys.org
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.
 
It works now, there was a problem with my IV. Just had to pick a correct value.Niles.