%Zwei gleich grosse Matrizen x,y werden als Vektoren interpretiert
%und davon wird das Skalarprodukt gebildet.
function s=skalarpr(x,y);
  [zx sx]=size(x);
  [zy sy]=size(y);
  if zx==zy & sx==sy
    n=zx*sx;
    xv=reshape(x,n,1);
    yv=reshape(y,n,1);
    s=xv'*yv; else
   disp('skalarpr: matrices must have the same size');
  end; %if 
return;
