function [S] = SkalarProdukt(a,b)

N = length(a);  % Länge des Vektors a
M = length(b);  % Länge des Vektors b

if N~=M
    fprintf('\nFEHLER! a und b müssen dieselbe Länge haben!\n\n');
    S = [];
    return;
else
    S = 0; % Wert des Skalarprodukts von a und b
    
    for j=1:N
        S = S + a(j)*b(j);
    end
end