Solución numérica de las ecuaciones del movimiento

Ecuación diferencial de segundo orden

Resolver la ecuación diferencial de segundo orden

$\frac{d^{2} θ}{d t^{2}} = - \frac{g (β + 1)}{l} \frac{sin θ}{{(1 + 4 β (1 + β) {sin}^{2} (\frac{θ}{2}))}^{3 / 2}}$

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

public class Oscilador extends RungeKutta{
    double beta; 
    public Oscilador(double beta,  double h){
      super(h);
      this.beta=beta;
    }
    public double f(double x, double v, double t){
         double r=Math.sqrt(1+4*beta*(1+beta)*Math.sin(x/2)*Math.sin(x/2));
         double y=-9.8299*(beta+1.0)*Math.sin(x)/(beta*6.37e6*r*r*r);
         return y;
    }
}
//Objetos de la clase Oscilador
 	Estado estado=new Estado(0.0, angulo, 0.0);
	Oscilador oscilador=new Oscilador(beta, 10.0);
      
  //rutina que calcula la trayectoria paso a paso
     oscilador.resolver(estado);