%
%   LPOLY_DEMO1A.M   Called by QUADL to determine Pm*Pn integrals
%
%   File prepared by J. R. White, UMass-Lowell (Aug. 2003)
%
      function f = ifile(x,M,N)
      APm = legendre(M-1,x);   APn = legendre(N-1,x);
      f = APm(1,:).*APn(1,:);   % recall that the 1st row is the Legendre poly
%
%   end of function
