I have a matrix A
A
A = [ 124.6,95.3,42.7 ; 95.3,55.33,2.74 ; 42.7,2.74,33.33 ]
The eigenvalues and vectors:
[V,D] =
To compute all dot products between columns of V:
V
M = squeeze(sum(bsxfun(@times, conj(V), permute(V, [1 3 2]))));
The columns of V (eigenvectors) will be orthogonal if the above M is diagonal.
M