function [c] = dprod(a, b)
%DPROD Dot product
%   C = DOT(A,B) returns the scalar product of the vectors A and B.
%   A and B must be vectors of the same length.

c = sum(a(:) .* b(:));

end

