# L(n,x) is Laguerre polynomials; laguerre(a_,b_,x_):=if(a<0 and a+b==-1, exp(x)); laguerre(a_,0,x_):=laguerre(a,x); laguerre(1,b_,x_):=b+1-x; laguerre(0,b_,x_):=1; laguerre(n_,x_):=sum((-1)^k*binomial(n,k)*x^k/k!,k,0,n); laguerre(0,x_):=1; laguerre(1,x_):=1-x; laguerre(-1,x_):=exp(x);