Anterior Ecuación diferencial de segundo orden

Resolver la ecuación diferencial de segundo orden

m d 2 x d t 2 =−2k( 1− l 0 x 2 + l 0 2 )x

con las condiciones iniciales siguientes: en el instante inicial t=0, x=A, y v=dx/dt=0, siendo A la amplitud de la oscilación.

public class Oscilador extends RungeKutta{
    double lonMuelle;
    double cte;
     public Oscilador(double lonMuelle, double cte, double h){
      super(h);
      this.lonMuelle=lonMuelle;
      this.cte=cte;
    }
    public double f(double x, double v, double t){     
         double aceleracion=-2*cte*x*(1.0-lonMuelle/
Math.sqrt(x*x+lonMuelle*lonMuelle));
         return aceleracion;
    }
} 
//Objetos de la clase Oscilador
 	Estado estado=new Estado(0.0, x0, 0.0);
	Oscilador oscilador=new Oscilador(lonMuelle, cte, 0.01);
      
  //rutina que calcula la trayectoria paso a paso
     oscilador.resolver(estado);
Anterior