To solve non-linear equations using matrix method, they must still be linearized.
The coefficients of the linearized equations (specifically the partial derivatives taken with respect to the unknown variables) will be gathered to form the Jacobian Matrix.
Let X be the matrix of unknown corrections dx and dy, and K be the matrix of constants.
