想解一个非线性方程组，求非负解，自己试着编了程序，现有如下疑问，特请教高人指点。程序代码如下：function F = myfun1(x)
F=[x(1)*cos(x(1))^2+x(2)*cos(x(2))^2-x(1);x(1)*sin(x(1))^2+x(3)*cos(x(3))^2-x(2);x(2)*sin(x(2))^2+x(3)*sin(x(3))^2-x(3);x(1)+x(2)+x(3)-1]

x0=[1;0;0];
options=optimset('display','iter');
[x,fval]=fsolve('myfun1',x0,options)
结果执行如下：
x =
    0.6807
    0.2944
    0.0248
但是当我设初值为x0=[0;0;1]时，结果为：
x =
   -0.4704
   -0.0971
    1.5674

请问如何设置才能使得在设置可以使得在初值x0=[0;0;1]的情况下，得到非负解？