syms a b c d e f;
Ia = (d.*f-e.^2)
Ib = (c.*e-b.*f)
Id =(a.*f-c.^2)
Ic= (b.*e-c.*d)
Ie=(b.*c-a.*e)
If=(a.*d-b.^2)
[Ia,Ib,Ic;Ib,Id,Ie;Ic,Ie,If]
expand([Ia,Ib,Ic;Ib,Id,Ie;Ic,Ie,If][a ,b,c;b,d,e;c,e f])/ (- fb^2 + 2bce - dc^2 - ae^2 + ad*f)
syms a b c d e f;
Ia = (d.*f-e.^2)
Ib = (c.*e-b.*f)
Id =(a.*f-c.^2)
Ic= (b.*e-c.*d)
Ie=(b.*c-a.*e)
If=(a.*d-b.^2)
[Ia,Ib,Ic;Ib,Id,Ie;Ic,Ie,If]
expand([Ia,Ib,Ic;Ib,Id,Ie;Ic,Ie,If][a ,b,c;b,d,e;c,e f])/ (- fb^2 + 2bce - dc^2 - ae^2 + ad*f)