ode_o1.h

#ifndef _ODE_O1_H_
#define _ODE_O1_H_

// a * df/dt + b * f = c
// (f(t + dt) - f(t)) / dt = g(f(t)) => f(t + dt) = f(t) + dt * g(f(t)).

double advance(double (*g)(double, void *), double f0, double dt, void *state);

#endif