function table = chebyfit(y,opd)

cf = cheb(y);
yfit = cf(:,1:2:13);
for k=1:7
    c = yfit(:,1:k)\opd;
    table(k,1:k)=c';
end