Anterior

Ecuación diferencial de segundo orden

Resolver la ecuación diferencial de segundo orden

d 2 r d t 2 = L 2 m(m+M) r 3 − M m+M g

por procedimientos numéricos con las condiciones iniciales: con las condiciones iniciales:en el instante t=0, la velocidad radial de la masa deslizante es dr/dt=0, y su distancia al eje r=r0.

public class Sistema extends RungeKutta{
    double mBloque;
    double mAngular;
// la masa de la partícula deslizante es 1 kg
    public Sistema(double mBloque, double mAngular, double h){
      super(h);
      this.mBloque=mBloque;
      this.mAngular=mAngular;
    }
    public double f(double x, double v, double t){     
         double aceleracion=mAngular*mAngular/((1.0+mBloque)*(x*x*x))-
         mBloque*9.8/(1.0+mBloque);
         return aceleracion;
    }
}

Se establece el estado incial 

Estado estado=new Estado(0.0, pos, 0.0);

Se crea un objeto de la clase derivada 

double mAngular=pos*pos*vAngular;
Sistema sis=new Sistema(mBloque, mAngular, 0.005);

Se llama a la función resolver que determina el estado del sistema en el instante t+h conocido el estado en el instante t 

sis.resolver(estado);

Anterior