D=[3 2 1 0
    2 4 -2 1
    10 17 21 -1
    5 14 -3 0]
b=[1;9;-8;2]
determinant=det(D)
transponeeritud=D'
poord=inv(D)
transvektor=b'
Q=D^2+D.*D-4.5*D'
korrutis=Q*b
(D(2,3)-D(3,3)^2)/(D(1,4)+D(3,2))*(D(2,1)+Q(1,2))-b(1)+b(2)-2*b(3)