clear all
A=[1 1 1]; %coefficienti di P(x) di grado length(A)-1
B=[1 1 1]; % monomio di grado 2 (length(B)-1)
k_A=find(A | not(A));  % vettore degli indici di A (tutti)
k_B=find(B | not(B));  % vettore degli indici B (tutti)
gradoDiA=length(A)-1;
gradoDiB=length(B)-1;
gradoDiC=gradoDiA+gradoDiB;
%B=A(2:end).*k(1:end-1) %Coefficienti di P'(x)
C=zeros(1,gradoDiA+gradoDiB+1);
for k = 0:gradoDiB
    M=zeros(1,k+1);    %inizializzazione monomio di grado k
    M(end)=B(k+1);     % monomio di grado k di B 
    D=prodottoPperM(A,M);
    C(1:length(D))=C(1:length(D))+D;   % Calcolo dei coefficienti di P(x) * B(x)
end

fprintf("I coefficienti dl polinomio P sono:"); disp(C)