Multiplicação de polinômios no Matlab sem o conv

Implementando, com auxílio dos conhecimentos do post anterior, a multiplicação de dois polinômios (que, normalmente, é feita pela função conv):

1. function D=multiplicarPolinomios(A,C)
2. % Entre com A e C na forma de vetor coluna e obtenha D, onde:
3. % A representa a0^N+...+aN
4. % B representa b0^M+...+bM
5. %e a saida eh a0*b0^(M+N)+...+aNbM
6.  
7. N=length(C);
8. B=zeros(0);
9. for i=1:N,
10.     Aux=[zeros(i-1,1);A;zeros(N-i,1)];
11.     B=[B,Aux];
12.     i=i+1;
13. end
14. D=B*C;

Precisei implementar essa função por extenso para calcular um sistema de equações, que não estava na forma explícita.