% Aitkens Method

clear;
% Generating an arbitrary sequence {pn}

itmax=input('itmax=');
for l=1:itmax
    ps(l)=feval('pn',l);  
end
% Sequence
disp(ps');

% Generating an aitkens seq. from ps

for l=1:itmax-2
    pa(l)=ps(l)-(ps(l+1)-ps(l))^2/(ps(l+2)-2*ps(l+1)+ps(l)); 
 end
 pa(itmax-1)=0;
 pa(itmax)=0;

disp([ps',pa']);