alpha = 150;

N = 100;
n = 4;
x0 = 5;

for i=1:N
    if(i == 1)
        x(i) = 1/2*(x0+alpha/x0^(n-1));
    else
        x(i) = 1/2*(x(i-1)+alpha/x(i-1)^(n-1));
    end
end

plot(x);
x(length(x))