# 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);