function [euestimate, realval, error] = e_explicit(f, df, y0, t0, tf, h)
