Anterior

Procedimiento numérico

Ecuación diferencial de segundo orden

Resolver la ecuación diferencial de segundo orden

d 2 θ d t 2 = g R sin⁡θ− μ R ( g cos⁡θ−R ( dθ dt ) 2 )

por procedimientos numéricos con las condiciones iniciales t=0, θ=0, dθ/dt=v0/R

public class Particula extends RungeKutta{
    final double radio=1.0;
    double roza;
    public Particula(double roza, double h){
        super(h);
        this.roza=roza;
    }
    public double f(double x, double v, double t){
        double temp=(9.8/radio)*(Math.sin(x)-roza*Math.cos(x))+roza*v*v;
        return temp;
    }

}

Se establece el estado incial 

    double vAngular=v0/radio;
    Estado estado=new Estado(0.0, 0.0, vAngular);

Se crea un objeto de la clase derivada 

 Particula p=new Particula(roza, 0.0025);

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 

p.resolver(estado);
Anterior