#topoly() convert to polymonial; topoly(bernoulli(n_,x_) ) := sum(binomial(n,k)*bernoulli(n-1,k)*x^k,k,0,n); topoly(harmonic(-5,x_)) := (-1/12)*x^2+5/12*x^4+1/2*x^5+1/6*x^6; topoly(harmonic(-4,x_)) := (-1/30)*x+1/3*x^3+1/2*x^4+1/5*x^5; topoly(harmonic(-3,x_)) := x^4/4+x^3/2+x^2/4; topoly(harmonic(-2,x_)) := x^3/3+x^2/2+x/6; topoly(harmonic(-1,x_)) := (x^2+x)/2; topoly(hermite(n_,x_)) :=if(n>1, expand(2x*hermite(n-1,x)-2*(n-1)*hermite(n-2,x)) );