Perform two iterations of the multivariate Newton-Raphson method to approximate a solution of the following system of nonlinear equations
f1 (x, y) = x^2 + xy - 10 = 0
f2 (x, y) = y + 3xy^2 - 57 = 0
Use initial values x = 1.5, y = 3.5. Include calculation of errors