function table = taylorfit(y,opd)

yfit =  [ ones(size(y)) y.^2 y.^4 y.^6 y.^8 y.^10];
sz=size(yfit);
for k=1:sz(2)
    c = yfit(:,1:k)\opd;
    table(k,1:k)=c';
end

Published with MATLAB® 7.6