%eks 10.5

n=4;
h=1/n;
y=zeros(n+1,1);
y(1)=1;

for i = 1:n
    k_1=-y(i);
    k_2=-(y(i)+h*k_1/2);
    k_3=-(y(i)+h*k_2/2);
    k_4=-(y(i)+h*k_3);
    y(i+1)=y(i)+h*(k_1+2*k_2+2*k_3+k_4)/6;
end

disp(y(end)-exp(-1))
