Use an adapted Heun's method to get an evaluation of the local error for
each step.

Signed-off-by: Raphael Catolino <raphael.catolino@gmail.com>