(b/sqrt(a^2 + b^2) | -a/sqrt(a^2 + b^2) -sin(t) | cos(t))

line | T^^(t) = (b/sqrt(a^2 + b^2), -a/sqrt(a^2 + b^2)) circle | T^^(t) = (-sin(t), cos(t))

line | x(t) = b t y(t) = a (-t) - c/b circle | x(t) = a cos(t) y(t) = a sin(t)

line | a x + b y + c = 0 circle | x^2 + y^2 = a^2

line | r(θ) = -c/(a cos(θ) + b sin(θ)) circle | r(θ) = a

-(b λ)/sqrt(a^2 + b^2) - (a sin(t))/sqrt(a^2 + b^2) + (b cos(t))/sqrt(a^2 + b^2) + λ^2 - λ cos(t)

λ_1 = (sqrt(a^2 + b^2) cos(t) - sqrt(4 a sqrt(a^2 + b^2) sin(t) + (a^2 + b^2) cos^2(t) - 2 b sqrt(a^2 + b^2) cos(t) + b^2) + b)/(2 sqrt(a^2 + b^2))

λ_2 = (sqrt(a^2 + b^2) cos(t) + sqrt(4 a sqrt(a^2 + b^2) sin(t) + (a^2 + b^2) cos^2(t) - 2 b sqrt(a^2 + b^2) cos(t) + b^2) + b)/(2 sqrt(a^2 + b^2))

v_1 = (-(-sqrt(a^2 + b^2) cot(t) + b csc(t) - csc(t) sqrt((4 a sqrt(a^2 + b^2) - 2 b sqrt(a^2 + b^2) cot(t) + a^2 cos(t) cot(t) + b^2 cos(t) cot(t) + b^2 csc(t)) sin(t)))/(2 sqrt(a^2 + b^2)), 1)

v_2 = (-(-sqrt(a^2 + b^2) cot(t) + b csc(t) + csc(t) sqrt((4 a sqrt(a^2 + b^2) - 2 b sqrt(a^2 + b^2) cot(t) + a^2 cos(t) cot(t) + b^2 cos(t) cot(t) + b^2 csc(t)) sin(t)))/(2 sqrt(a^2 + b^2)), 1)

(b/sqrt(a^2 + b^2) | -a/sqrt(a^2 + b^2) -sin(t) | cos(t)) = P.D.P^(-1) where P = ((0.5 (sqrt(a^2 + b^2) cot(t) + csc(t) sqrt(4 a sqrt(a^2 + b^2) sin(t) + (a^2 + b^2) cos^2(t) - 2 b sqrt(a^2 + b^2) cos(t) + b^2) - b csc(t)))/sqrt(a^2 + b^2) | 0.5 cot(t) - (0.5 csc(t) (sqrt(4 a sqrt(a^2 + b^2) sin(t) + (a^2 + b^2) cos^2(t) - 2 b sqrt(a^2 + b^2) cos(t) + b^2) + b))/sqrt(a^2 + b^2) 1 | 1) D = ((0.5 (sqrt(a^2 + b^2) cos(t) - sqrt(4 a sqrt(a^2 + b^2) sin(t) + (a^2 + b^2) cos^2(t) - 2 b sqrt(a^2 + b^2) cos(t) + b^2) + b))/sqrt(a^2 + b^2) | 0 0 | (0.5 (sqrt(a^2 + b^2) cos(t) + sqrt(4 a sqrt(a^2 + b^2) sin(t) + (a^2 + b^2) cos^2(t) - 2 b sqrt(a^2 + b^2) cos(t) + b^2) + b))/sqrt(a^2 + b^2)) P^(-1) = 1/sqrt(4 a sqrt(a^2 + b^2) sin(t) + (a^2 + b^2) cos^2(t) - 2 b sqrt(a^2 + b^2) cos(t) + b^2)(1 sqrt(a^2 + b^2) sin(t) | 0.5 (-1 sqrt(a^2 + b^2) cos(t) + 1 sqrt(4 a sqrt(a^2 + b^2) sin(t) - 2 b sqrt(a^2 + b^2) cos(t) + a^2 cos^2(t) + b^2 cos^2(t) + b^2) + 1 b) -1 sqrt(a^2 + b^2) sin(t) | 0.5 (1 sqrt(a^2 + b^2) cos(t) + 1 sqrt(4 a sqrt(a^2 + b^2) sin(t) - 2 b sqrt(a^2 + b^2) cos(t) + a^2 cos^2(t) + b^2 cos^2(t) + b^2) - 1 b))