MATLAB: Solving odes (restricted 3 body problems)


I have been trying to solve these equations -

d2x dt2 = 2 dy dt + x− M(x + E) r3 E − E(x−M) r3 M

d2y dt2 = −2 dx dt + y− My r3 E − Ey r3 M

z1 = x, z2 = dx dt , z3 = y, and z4 = dy dt

but keep getting these errors:-

Unrecognized function or variable 'rE'.

Error in secodef (line 4) zprime=[z(2);2*z(4)+z(1)-((M*(z(1)+E))/rE^3)-((E*(z(1)-M))/rM^3);z(4);-2*z(2)+z(3)-((M*z(3))/rE^3)-((E*z(3))/rM^3)];

Error in odearguments (line 90) f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.

Error in ode45 (line 115) odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);

Error in trial1 (line 13) [t,z] = ode45(@secodef,tspan,Z);

what do I do?

Best Answer

Variables rE and rM should be in the function

Related Solutions

MATLAB: Solving odes (resiscted 3 body problem)

Building matrices with [ ] can be picky when it comes to spaces. E.g.,

>> [1;5-3;4]
ans =
     1
     2
     4
>> [1;5 -3;4]
??? Error using ==> vertcat
CAT arguments dimensions are not consistent.
>> [1;(5 -3);4]
ans =
     1
     2
     4

In the second case above, because of that extra space the parser interpreted the "5 -3" as a 1x2 vector and not as a single scalar. To fix this, either get rid of all of your spaces, or use parentheses as shown.

Side Note: Your line 10 seems to be missing some E's and M's when compared to your posted equations. You might check that.

MATLAB: Solving a system of 4 second order differential equations

I believe that re and rm should be functions too

re = @(z) sqrt((z(1)+E)^2+z(3)^2);
rm = @(z) sqrt((z(1)-M)^2+z(3)^2);
zprime = @(t,z) [z(2)
    2*z(4)+z(1)-((M*(z(1)+E))/re(z)^3)-((E*(z(1)-M))/rm(z)^3)
    z(4)
    -2*z(2)+z(3)-((M*z(3))/re(z)^3)-((E*z(3))/rm(z)^3)];
tspan = [0 5e-3];
z0 = [0 1 0 1];
[t, z] = ode45(zprime,tspan,z0);

Pay attention

Best Answer

Related Solutions

MATLAB: Solving odes (resiscted 3 body problem)

MATLAB: Solving a system of 4 second order differential equations

Related Question