The code attached below is used to find the four lowest non-zero natural frequencies of a system. It produces the correct values for the natural frequencies and for mode shapes 1, 2 and 4. However, mode shape 3 is incorrect and I’ve been having trouble figuring out why. Any help would be greatly appreciated.

NOTE: matrix A is a 13 by 13 matrix
Relevant MATLAB code