function[c]=produse(s1,s2)
[m,n]=size(s1);
[x,y]=size(s2);
k=n;
 if m~=y
   error('Matricele s1 si s2 nu se pot inmulti!')
 end
   for i=1:n
      for j=1:n
       c[i,j]=0;
for k=1:n
       c[i,j]=c[i,j]+s1[i,k]*s2[k,j];
      end
   end
end
end