alpha = 2;
x0 = 2;

for n=2:5
    f = section4(n,x0,alpha);
    plot([x0 f]);
    [n f(length(f))]
    hold on;
end

