Descenso de un paracaidista

Procedimiento numérico

Resuelve la ecuación diferencial de segundo orden

$\frac{d^{2} x}{d t^{2}} = - g + \frac{k_{0}}{m} v^{2} \exp (- \frac{x}{λ})$

por procedimientos numéricos con las siguientes condiciones iniciales t=0, x=x₀, (dx/dt)=0

public class Sistema extends RungeKutta{
    double masa;
    double k;
    Sistema(double masa, double k, double h){
      super(h);
      this.masa=masa;
      this.k=k;
    }
    public double f(double x, double v, double t){
         double y=-9.8+k*v*v*Math.exp(-x/7482.2)/masa;
         return y;
    }
}

Se establece el estado incial

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

Se crea un objeto de la clase derivada

 Sistema sis=new Sistema(masa, k, 0.5);

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);