Question about using Newton's Method to solve a system of equations

[Ax_xiom]

So the formula used to solve non-linear equations using the Newton-Raphson method is this $$ X_{k+1} = X_k - J^{-1}(X_k)F(X_k) $$ and instead of finding $J^{-1}(X_k)$, we solve for the change that will be applied to $X_k$ ($\Delta X_k$) using this relation $$J(X_k) \Delta X_k = -F(X_k)$$ and the next iteration will be this $$ X_{k+1} = X_k + \Delta X_k $$ My question is that wouldn't you need to solve for $J^{-1}(X_k)$ anyways when doing that? If so, why do we do that step in the first place?

 

[pasmith]

We don't solve linear systems numerically by finding an inverse matrix; that is incredibly inefficient. Instead we use other methods, like Gaussian elimination or LU decomposition.

 

[Ax_xiom]

We don't solve linear systems numerically by finding an inverse matrix; that is incredibly inefficient. Instead we use other methods, like Gaussian elimination or LU decomposition.

So what happens comptationally during the $J(X_k) \Delta X_k = -F(X_k)$ step? Do you express the problem as a series of linear equations and solve them?

Edit: I think that might be the case. I was under the impression that using Gaussian manipulation to solve a system of equations would give you the inverse matrix for free, but I don't think that is the case