function res = angle_between_vectors(v1, v2) res = acos( (v1'*v2)/ (norm(v1)*norm(v2)));