El péndulo

Procedimientos numéricos

Ecuación diferencial de segundo orden

Resolver la ecuación diferencial de segundo orden

$\frac{d^{2} θ}{d t^{2}} + 2 γ \frac{d θ}{d t} + ω_{0}^{2} \sin θ = f \cos (ω t)$

con las condiciones iniciales: t=0, dθ/dt=v₀, θ=x₀

public class Oscilador extends RungeKutta{
    final double w0=1.0;
    double w;
    double fuerza;
    double gamma;
    public Oscilador(double w, double fuerza, double gamma, double h){
      super(h);
      this.w=w;
      this.fuerza=fuerza;
      this.gamma=gamma;
    }
    public double f(double x, double v, double t){
         return (-w0*w0*Math.sin(x)-gamma*v+fuerza*Math.cos(w*t));
    }
} 
//Objetos de la clase Oscilador
 	Estado estado=new Estado(0.0, x0, v0);
	Oscilador oscilador=new Oscilador(w, fuerza, gamma, 0.02);
      
  //rutina que calcula la trayectoria paso a paso
     oscilador.resolver(estado);