# laguerreplus(n,x) is Laguerre plus polynomials;

laguerreplus(a_,b_,x_):=if(a<0 and a+b==-1, exp(x));
laguerreplus(a_,0,x_):=laguerreplus(a,x);
laguerreplus(1,b_,x_):=b+1+x;
laguerreplus(0,b_,x_):=1;

laguerreplus(n_,x_):=sum(binomial(n,k)*x^k/k!,k,0,n);
laguerreplus(0,x_):=1;
laguerreplus(1,x_):=1+x;
laguerreplus(-1,x_):=exp(x);