%%
tol = 10;

err = 1000;

u = 0.4;

i = 1 ;

while err > tol
    
    [a, b] = eu_implicit(@(x, y) 3*exp(-201*x) - 5*exp(-x), 201, -2, 0, 1, u);
    
    err = b;
    
    u = 1/i;
    
    i = i + 1;
    
end

disp(['The number of steps is ', num2str(u)]);


