function solve_nonl
clear all;clc
format long
global w
a=0:0.1:10;
n=length(a);
x0=[0.1 0.1 -1];
for i=1:n
    w=a(i);
    xresult=fsolve(@nonl,x0);
    xx(i)=xresult(1);
    yy(i)=xresult(2);
    zz(i)=xresult(3);
end
[a' xx' yy' zz']

plot(a,xx,'k',a,yy,'b',a,zz,'r')

function f=nonl(x)
global w
f(1)=3*x(1)-cos(x(2)*x(3))-1/2;
f(2)=x(1)*x(1)-81*(w*x(2)+0.1)^2+sin(x(2))+1.06;
f(3)=exp(-x(1)*x(2))+20*x(3)+(10*pi-3)/3;