OPT=optimset; OPT.LargeScale='off';
[x,Y,c,d] = fsolve(@my2deq,[1; 2],OPT),

f=inline('[p(1)*p(1)+p(2)*p(2)-1; 0.75*p(1)^3-p(2)+0.9]','p');
[x,Y] = fsolve(f,[1; 2],OPT)

[x,Y,c,d] = fsolve(f,[-1,0]',OPT); x, Y, kk=d.funcCount

syms x y;
X=solve('x^2+y^2-1=0','75*x^3/100-y+9/10=0')

x1=X.x

y1=X.y

x=x1(2); y=y1(2); [eval('x^2+y^2-1') eval('75*x^3/100-y+9/10')]