fsolve is not returning the right solutions

My code is below. The correct solutions must be (1,0,1) and (-1,0,1) but I keep getting weird answers regardless of my initial guess

function solveeqs()
x0 = [1 1 1];
x = fsolve(@eqns,x0);
x;
eqns(x0)
end

function F = eqns(z)
x1 = z(1);
x2 = z(2);
x3 = z(3);
F(1) = x1 - x1*x2 - x2.^3 + x3*x1.^2 + x3*x2.^2 - x3 - x3*x1 + x3*x1*x2 + x3*x2.^3;
F(2) = x1 - x3*x1 + x3*x2 + 2*x3*x1*x2;
F(3) = 1 + 2*x2.^2 + x3.^2;
end

Answer 1

[Alex Sha]

Hi, if you think the solutions should be (1,0,1) or (-1,0,1) , then one equation:

F(3) = 1 + 2*x2.^2 + x3.^2;

should be:

F(3) = 1 + 2*x2.^2 - x3.^2;

Answer 2

[Matt J]

The solutions you've cited definitely do not solve the equations you've provided,

F=eqns([-1,0,1])

F = 1×3

0 0 2

F=eqns([1,0,1])

F = 1×3

0 0 2

function F = eqns(z)
x1 = z(1);
x2 = z(2);
x3 = z(3);
F(1) = x1 - x1*x2 - x2.^3 + x3*x1.^2 + x3*x2.^2 - x3 - x3*x1 + x3*x1*x2 + x3*x2.^3;
F(2) = x1 - x3*x1 + x3*x2 + 2*x3*x1*x2;
F(3) = 1 + 2*x2.^2 + x3.^2;
end