If anyone could help with a matlab code working with newton raphson to solve this

In dimensionless units, the continuous stirred tank reactor equations can be expressed as: dtdx1=1x1exp(x2)x1dtdx2=1x2+exp(x2)x1x2 Assume that =100,=5,=149.39 and =0.553. - You need to implement a Newton-Raphson multivariable algorithm. Find the steady state solutions. Your tolerance should be 1012. - Plot your error (norm of solution or/an function residual) as a function of the number of iterations. - Compare your N-R solution with a Matlab solution using a root finding function available in Matlab. Turn in your code for N-R. Hint: you should have 3 solutions. Steady state means no time dependence