636
L_x=sqrt(t(x) %*% x) or L_x=sqrt(sum(x^2)) theta=acos(t(x) %*% y) / (L_x * L_y)*360 / (2*pi) ### solve(cbind(x1,x2,x3),rep(0,3)) A-1=solve(A)