function [N] = solutions_Final_p1d()

    sol = exp(-1);

    [t0, y0] = ode45(@solutions_Final_p1a, [0 1], 0);
    e = abs(y0(end) - sol);
    
    N = length(t0);
    
    [~, y] = solutions_Final_p1c(N);
    while abs(y(end)-sol)>e
        N = N + 1;
        [~, y] = solutions_Final_p1c(N);
    end
    


end

