function [c] = my_prodmm(a,b)
%calcola c=a*b
n=length(a);
c=zeros(n);
%y=zeros(n,1);
%z=zeros(n,1);
for j=1:n
%     for k=1:n
%         y(k)=b(k,j);
%     end
%y=b(:,j);%y=b(1:n,j)
c(:,j)=my_prodmv(a,b(:,j))
% for k=1:n
%     c(k,j)=z(k);
% end
%c(:,j)=z;
end
end