Sun-Earth 3-body problem

[Gerry]

Hi people,
I am a new user of MatCont.

1) Is anybody of you using MatCont to calculate periodic orbits (halo, lyapunov orbits etc)in a 3 body problem?? I am asking to understand if it is possible using MatCont and where to start from.

2) In a continuation of equilibrium points for the three-body problem with radiation pressure I always have the error message "Evaluation of test functions failed at starting point" with the status in the MatCont window set to Current step size too small.
any suggestion to overcome such a problem?
I am sure I started from an equilibrium point (the lagrangian L1 point) and I tried to change all the tolerances, but nothing happened.

[Gerry]

Here are the equations

x'=xp
y'=yp
z'=zp
xp'= x + 2*yp - (mu_e*(mu_e + x - 1))/((mu_e + x - 1)^2 + y^2 + z^2)^(3/2) + ((mu_e + x)*(mu_e - 1))/((mu_e + x)^2 + y^2 + z^2)^(3/2) - (B0*(mu_e + x)*(mu_e - 1))/((mu_e + x)^2 + y^2 + z^2)^(3/2)
yp'= y - 2*xp - (mu_e*y)/((mu_e + x - 1)^2 + y^2 + z^2)^(3/2) + (y*(mu_e - 1))/((mu_e + x)^2 + y^2 + z^2)^(3/2) - (B0*y*(mu_e - 1))/((mu_e + x)^2 + y^2 + z^2)^(3/2)
zp'= (z*(mu_e - 1))/((mu_e + x)^2 + y^2 + z^2)^(3/2) - (mu_e*z)/((mu_e + x - 1)^2 + y^2 + z^2)^(3/2) - (B0*z*(mu_e - 1))/((mu_e + x)^2 + y^2 + z^2)^(3/2)

where mu_e = 3.0359e-6 and B0 is a parameter that varies the radiation pressure.

Starting from the equilibrium point for B0=0
[0.989990937176541, 0, 0, 0, 0, 0]

and using B0 as free parameter for the equilibirum point continuation, I have the above-mentioned error message.

Thank you for your help.

[cc]

same question, "Evaluation of test functions failed at starting point" ocured when I do the limit cycle continuation