From: samah karim on
E = 0.5 *m *(dr/dt)^2 + L^2/2mr^2 - a/r;
L = m r^2 (d theta/dt);
for a given central force field(L,E), an ODE solver is supposed to do the trick. i coded my own euler but it's giving me the wrong result. im supposed to get an ellipse as the orbit. and im not really familiar with the RK4 in MATLAB
From: James Tursa on
"samah karim" <sammouha_007(a)hotmail.com> wrote in message <hmup8f$6f6$1(a)fred.mathworks.com>...
> E = 0.5 *m *(dr/dt)^2 + L^2/2mr^2 - a/r;
> L = m r^2 (d theta/dt);
> for a given central force field(L,E), an ODE solver is supposed to do the trick. i coded my own euler but it's giving me the wrong result. im supposed to get an ellipse as the orbit. and im not really familiar with the RK4 in MATLAB

Can you post your code?

James Tursa
 | 
Pages: 1
Prev: Binary input
Next: GUI output again