dirmul

dirmul([0,1],[0,1])
dirdiv([0,1,1,1,1],[1])
dirmul(vector(10,n,moebius(n)),vector(10,n,1))
dirmul([1,2,3,4],[1,2,3,4])
dirdiv([1,2,3,4],[1,2,3,4])
dirdiv([1,2,3,4],2*[1,2,3,4])

direuler(p=1.5,10.5,1/((1-X)*(1-p*X)))
direuler(p=2,10,1/((1-X)*(1-p*X)), -1)
direuler(p=2,10,1/((1-X)*(1-p*X)), 1)
direuler(p=2,10,1/((1-X)*(1-p*X)), 2)
direuler(p=2,10,1/((1-X)*(1-p*X)), 3.5)
direuler(p=2,10,-1/(X^2-1))
direuler(p=2,10,(X^2-1)/(X^3+X-1))

direuler(p=I,10.5,1/((1-X)*(1-p*X)))
direuler(p=2.5,I,1/((1-X)*(1-p*X)))
direuler(p=2,10,1/((1-X)*(1-p*X)), I)
direuler(p=2,1,1/((1-X)*(1-p*X)), 1)
direuler(p=2,2,1/((1-X)*(1-p*X)), 2)
direuler(p=2,2,1/((1-X)*(1-p*X)), 5)

direuler(p=2,10,1/2)
direuler(p=2,10,1/(x-2))