function [ baller, er ] = eu_explicit(func, dfunc, y0, t0, tf, h)

t = t0;

y = y0;

N = (tf-t0)/h;

for j = 1:N
    
    y = y + dfunc(t,y)*h;

    er = abs(((func(t) - y)/func(t))*100);

    t = t + h;
end

baller = y;

