function s = mydot(x,y)
% Euclidean inner product
s = 0; for j=1:length(x), s = s+x(j)*y(j); end