function druga()
    A = [1 2 3; 2 3 1; 3 2 1]
    B = [1 1; 2 -1; 2 -4]

    # Imamo torej B * X = A
    # Ker je sistem predoločen uporabimo metodo min kvadratov: B' * B * X = B' * A
    # Nato Q, R = QR(B' * B) in R * X = Q' * (B' * A)
    # X = inv(R) * Q' * (B' * A)

    LEVA = B' * B
    DESNA = B' * A
    Q, R = qr(LEVA)

    Q_DESNA = Q' * DESNA

    X = R \ Q_DESNA  # oz. resi z vstavljanjem
    return X
end