Péndulo accionado por las fuerzas de marea

Ecuación diferencial de segundo orden

Resolver la ecuación diferencial de segundo orden

$\begin{array}{l} \frac{d^{2} θ}{d t^{2}} = \frac{G M \cdot d \cdot sin θ}{2 l} (\frac{1}{r_{1}^{3}} - \frac{1}{r_{2}^{3}}) \\ r_{1}^{2} = {(R + l \cdot \cos θ)}^{2} + {(l \cdot sin θ)}^{2} = R^{2} + l^{2} + 2 l R \cos θ \\ r_{2}^{2} = {(R - l \cdot \cos θ)}^{2} + {(l \cdot sin θ)}^{2} = R^{2} + l^{2} - 2 l R \cos θ \end{array}$

Con las condiciones iniciales siguientes: en el instante t=0, θ=π/6, dθ/dt=0.

public class Oscilador extends RungeKutta{
    double lon; 
    double r; 
    public Oscilador(double lon, double r, double h){
      super(h);
      this.lon=lon;
      this.r=r;
    }
    public double f(double x, double v, double t){
         double r1=1.0+lon*lon/(r*r)-2*lon*Math.cos(x)/r;
         double r2=1.0+lon*lon/(r*r)+2*lon*Math.cos(x)/r;
         double y=-(1.543e-6*Math.sin(x)/(2*lon*r*r))*
(1.0/Math.pow(r1, 1.5)-1.0/Math.pow(r2, 1.5));
         return y;
    }

} 
//Objetos de la clase Oscilador
 	Estado estado=new Estado(0.0, Math.PI/6, 0.0);
	Oscilador oscilador=new Oscilador(lon, r, 10.0);
      
  //rutina que calcula la trayectoria paso a paso
     oscilador.resolver(estado);