lagp.m

求解量子力学的薛定谔方程

function f=lagp(k,p)
%> <lagp.m> calculates lagp(k,p), the coefficients of the Laguerre
%> polynomials L(k,0) ... L(k,p)
%> Call P = lagp(k,p);  
%> input p,k non-negative integers
%> output a (p+1)*(p+1) matrix of coefficients
%> To evaluate the value of the polynomial:
%> Call pol = evalpol(P,x);
%> input  x scalar or row vector 
%> output a (p+1)*length(x) matrix evaluating all the polynomials, 
%> OR 	pol=polyval(P(n,:),x) for each row, 
%> input: x a matrix 
%> output: a matrix of same dimension
%>
%>