function y = vectdist (a, b)

  y = sqrt (sum ((a - b').^2));

return

