Please write a numerical solver using the multi$-$dimensional Newton$-$Raphson method to solve the following system of equations with starting values x$0=2$ and y$0=0\mathrm{.}5.$

Please write your own codes of the Newton$-$Raphson method, rather than using the MATLAB function fsolve$().$

Please run the solver for $10$ iterations and demonstrate how the values of x and y converge.

Particularly, please plot the values of x$,$ y$,$ and the norm of $[$dx$,$ dy$]$ as a function of the iteration number.

In MATLAB, the function of $$norm$($v$)$ returns the Euclidean norm of the vector v$.$

f$1($x$,$y$)=$ x$2+$ xy$25=0$

f$2($x$,$y$)=$ y $+3$xy $7=0$