Matlab - Resolver sistema de ecuaciones no lineales en Matlab

>> syms x y

>> digits(5)

>> u=x*cos(x)-y-1; v=x-y*sin(y)-0.5;

>> ux=diff(u,x); uy=diff(u,y);

>> vx=diff(v,x); vy=diff(v,y);

>> x=1;  y=1; n=0;

>> x=x-vpa(eval((u*vy-v*uy)/(ux*vy-uy*vx)),5), y=y-vpa(eval((v*ux-u*vx)/(ux*vy-uy*vx)),5), n=n+1

En la 8va, 9na y 10ma iteración me sale este resultado

x=x-vpa(eval((u*vy-v*uy)/(ux*vy-uy*vx)),5), y=y-vpa(eval((v*ux-u*vx)/(ux*vy-uy*vx)),5), n=n+1

x =

2.2067

y =

-2.3105

n =

    10

y en solve me sale así:

>> [x,y]=solve('x*cos(x)-y=1','x-y*sin(y)=0.5','x','y');

>> x=vpa(x,5)

x =

0.70646

>> y=vpa(y,5)

y =

-0.46262

f(x,y)=4x+6y  entonces f(2,3)=4(2)+6(3)=8+18=26

n                 x                           y

   0         1.00000000          1.00000000

   1        -0.18314024         -0.10337291

   2         0.58262365         -0.45266042

   3         0.69885430         -0.45370129

   4         0.70637897         -0.46259366

   5         0.70646428         -0.46261958

   6         0.70646428         -0.46261959