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

0.5833
-0.1736 - 0.4603i -0.1736 + 0.4603i 0.2215 + 0.4898i
0.2215 - 0.4898i
R3 =
4.0755 + 4.1517i

0 0

0

0 4.0755 - 4.1517i 0

0

0 0 -1.0755 + 2.7440i

0

0 0 0 -

1.0755 - 2.7440i

A3*U3 - U3*R3

ans =
1.0e-014 *
-0.2998 - 0.4885i -0.2998 + 0.4885i 0.0944 + 0.3553i
0.0944 - 0.3553i
0.1776 + 0.3109i 0.1776 - 0.3109i 0.0111 - 0.1055i
0.0111 + 0.1055i
0.0444 + 0.0888i 0.0444 - 0.0888i -0.0666 - 0.0666i -
0.0666 + 0.0666i
-0.1110 - 0.3109i -0.1110 + 0.3109i -0.1776 + 0.2970i -
0.1776 - 0.2970i

Again, withthe e-014 term this is essentially zero.

[U4, R4] = eig(A4)

U4 =
[ 1,

1]
[ -(-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]
R4 =
[ 1/2*d+1/2*a+1/2*(d^2-2*a*d+a^2+4*c*b)^(1/2),
0]