function y = orth_polyval(x,oc,ab)

m=length(x);
n=length(ab);
op = ones(n+1,m);  

op(2,:) = (x'-ab(1,1));  % Orth. Pol. P_1

for k=3:n+1 % Pol Grad 2...Deg
    op(k,:) = (x'-ab(k-1,1)).*op(k-1,:) - ab(k-1,2)*op(k-2,:);
end

y = (op'*oc);
end
