GSL
GSL (GNU Scientific Library)
Hace ya un tiempo que estoy utilizando la librería GSL y la verdad que me ha parecido una librería fantástica con una buena documentación. Aún así, sobre las ODE voy a dar un par de apuntes.
Código algo más legible Vamos a tomar el mismo ejemplo que en la documentación, un oscilador de Van der Pol con la ecuación Ecuación del oscilador de Van der Pol y el sistema Sistema Van del Pol Para hacer el código más legible lo haremos más cercano a nuestras ecuaciones int func (double t, const double y[], double yp[], void* params) { (void)(t); /* Evitar warning no usado */ double mu = *(double*)params; double u = y[0]; double v = y[1]; double ut = v; double vt = -u + mu*v*(1 - u*u); yp[0] = ut; yp[1] = vt; return GSL_SUCCESS; }