% Aitken-Steffensen method
% given function
f=inline('(10/(x+4))^(1/2)');
p0=1.5; iter=0;
for n=1:200
    p1=f(p0);
    p2=f(p1);
    p=p0-(p1-p0)^2/(p2-2*p1+p0);
    err=abs(p-p0);
    if err>=10^(-10)
        iter=iter+1;
        p0=p;
    else
        break
    end
end
answer=p, iterations=iter