function [ baller, er ] = eu_implicit(func, a, y0, t0, tf, h)

t = t0;

y = y0;

N = (tf-t0)/h;

true = -1.839396746023524;

for j = 1:N
    
    y = (y + h*(-1000*exp(-(t+h))))/(1 + a*h);

    er = abs(((true - y)/true)*100);

    t = t + h;
end

baller = y;

