Oscilaciones anarmónicas

Procedimientos numéricos

Ecuación diferencial de segundo orden

Resolver la ecuación diferencial de segundo orden

$\frac{d^{2} x}{d t^{2}} + ω_{0}^{2} x - ω_{0}^{2} s x^{2} = 0 ω_{0}^{2} = \frac{k}{m}$

con las condiciones iniciales: En el instante t=0, la posición inicial de la partícula es x₀ (a la derecha del origen) y parte del reposo, su velocidad inicial es dx/dt=0.

public class Oscilador extends RungeKutta{
    double frecAngular; 
    double parametro;
    public Oscilador(double frecAngular, double parametro, double h){
      super(h);
      this.parametro=parametro;
      this.frecAngular=frecAngular;
    }
    public double f(double x, double v, double t){
         return (-frecAngular*frecAngular*x+parametro*frecAngular*frecAngular*x*x);
    } 
} 
//Objetos de la clase Oscilador
 	Estado estado=new Estado(0.0, amplitud, 0.0);
	Oscilador oscilador=new Oscilador(Math.sqrt(k), parametro, 0.05);
      
  //rutina que calcula la trayectoria paso a paso
     oscilador.resolver(estado);