% P_MATRIX  The matrix P of order N+1
%    Be sure to append the file fxy_pkpl to the routine Ial.m
%
function y=P_matrix(N)
global abj abl k l
y=zeros(N+1);
abj=r_jacobi01(N+1); abl=r_jacobilog(N+1);
for k=0:N
  for l=0:N
    [Ialg,Ilog]=Ial(N+1);
    y(k+1,l+1)=-Ilog;
  end
end
