function theta = vectorangle(A,B,C);
V1 = A - B; V2 = C - B;
V = dot(-V1,V2);
theta = acos(V/(norm(V1)*norm(V2)));
