function table = zernfit(y,opd)

yfit = zern(y);
for k=1:7
    c = yfit(:,1:k)\opd;
    table(k,1:k)=c';
end