I$\text{'}$m unsure how to solve this problem, help would be greatly appreciated.

Consider the following systems of equations:

$y=-x^2-Ax+0\mathrm{.}75$

$y+Bxy=x^2$

where constants A and $B$ are parameters. Develope a function file that will solve for $x$ and $y$ given the values of the two parameters in the system using the Newton$-$Raphson iteration. Your function should accept the following inputs $($in order$)$:

Ascalar value for $A$

A scalar value for $B$

A column vector of initial guesses for $x_0$ and $y_0$

A stopping critierion for the iteration

Your function should return the following outputs $($in order$)$:

A $2-$element column vector of the solutions $x$ and $y$

The Jacobian matrix evaulated at the initial guesses $x_0$ and $y_0$

The Euclidean norm for the residuals vector associated with the solution.$|$

The number of iterations required for covergence with default maximum at $50.$

Note: Use an analytical formulation of the Jacobian matrix and do not rearrange the equations before computing partial derivatives if you want your code to pass the tests. You can use the referenced NRsys$.$m in your solution code.

Function $o.$

function $,$ norm$\_$res, iter$]=$ P$1\_$solver$(A,B,,$ es $)$