// Updated (13-12-06)

Ti=0.501;
A = [Ti 0 0 0 1-Ti];
abs(roots(A))