% P_JACOBI01  Jacobi polynomial on [0,1] of degree k
%
function y=p_jacobi01(x,k)
global abj
if k==0
  y=ones(size(x)); return
elseif k==1
  y=x-abj(1,1); return
else
  pm1=ones(size(x)); p=x-abj(1,1);
  for j=2:k
    pm2=pm1; pm1=p;
    p=(x-abj(j,1)).*pm1-abj(j,2)*pm2;
  end
  y=p;
end
