Troubleshooting - A Guide to MATLAB for Beginners and Experienced Users - page 263

[ 0, 1/2*d+1/2*a-

1/2*(d^2-2*a*d+a^2+4*c*b)^(1/2)]

A4*U4 - U4*R4

ans =
[
0,
0]
[c-d*(-1/2*d+1/2*a-1/2*(d^2-2*a*d+a^2+4*c*b)^(1/2))/b+(-
1/2*d+1/2*a-1/2*(d^2-
2*a*d+a^2+4*c*b)^(1/2))/b*(1/2*d+1/2*a+1/2*(d^2-
2*a*d+a^2+4*c*b)^(1/2)), c-d*(-1/2*d+1/2*a+1/2*(d^2-
2*a*d+a^2+4*c*b)^(1/2))/b+(-1/2*d+1/2*a+1/2*(d^2-
2*a*d+a^2+4*c*b)^(1/2))/b*(1/2*d+1/2*a-1/2*(d^2-
2*a*d+a^2+4*c*b)^(1/2))]

simplify(ans)

ans =
[0,0]
[0,0]

(b)

A=[102;-104;-1-15];
clear U1 U2 R1 R2
[U1, R1] = eig(A)

U1 =
-0.8165

0.5774 0.7071

-0.4082 0.5774 -0.7071

-0.4082 0.5774 0.0000

R1 =

2.0000 0 0

0 3.0000 0

0 0 1.0000

B = [5 2 -8; 3 6 -10; 3 3 -7];
[U2, R2] = eig(B)