Перевести код из maple в Wolfram Mathematica

> f1 := proc (x1, x2) options operator, arrow; x1-2*x2^2+1 end proc; f2 := proc (x1, x2) options operator, arrow; -x1^2+2*x2-1 end proc;
> F := proc (x1, x2) options operator, arrow; Matrix(2, 1, [f1(x1, x2), f2(x1, x2)]) end proc;
> plots[implicitplot]({f1(x, y) = 0, f2(x, y) = 0}, x = -10 .. 2, y = -10 .. 10, numpoints = 200000);
> W := proc (x1, x2) options operator, arrow; Matrix([[(D[1](f1))(x1, x2), (D[2](f1))(x1, x2)], [(D[1](f2))(x1, x2), (D[2](f2))(x1, x2)]]) end proc;
> with(LinearAlgebra);
> 
> x[0] := Matrix(2, 1, [1, 2]); e := 0.1e-2; k := 0;
> do x[k+1] := evalf(x[k]+Multiply(MatrixInverse(W(x[k][1, 1], x[k][2, 1])), -F(x[k][1, 1], x[k][2, 1]))); if max(abs(x[k+1][1, 1]-x[k][1, 1]), abs(x[k+1][2, 1]-x[k][2, 1])) <= e then break else k := k+1 end if end do;
> x0 := x[k+1];

Quit
 
f1[x_, y_] = Sin[x + 1] - y - 1
f2[x_, y_] = 2*x + Cos[y] - 2
 
ContourPlot[{f1[x, y] == 0, f2[x, y] == 0}, {x, -10, 10}, {y, -10, 10}]
 
FindRoot[{f1[x, y] == 0, 
  f2[x, y] == 0}, {{x, 0, 0, 10}, {y, -1, -1, 10}}]
 
F[x_, y_] = {{f1[x, y]}, {f2[x, y]}};
W[x_, y_] = {{D[f1[x, y], x], D[f1[x, y], y]}, {D[f2[x, y], x], 
   D[f2[x, y], y]}}
 
x[0] = MatrixForm[{{1}, {2}}]
e = 0.0001
k = 0
 
While[True; 
 x[k + 1] = 
  N[x[k] + Inverse[
      W[x[k][[1, 1]], x[k][[2, 1]]]].-F[x[k][[1, 1]], x[k][[1, 2]]]] ;
 If[Max[Abs[x[k + 1][[1, 1]] - x[k][[1, 1]]];
   Abs[x[k + 1][[1, 1]] - x[k][[1, 2]]] <= e, Break[], k = k + 1 ]]]
x0 = x[k + 1]