function z=func(x0,u)
% С помощью этой функции мы находим значение функции
global a p;
p1=polyval(p,x0(1));
z=[p1-a(1)*u(3);
a(2)*x0(1)*x0(3)-a(3)*x0(2)-a(4)*u(1);
a(5)*u(2)-a(6)*x0(1)*x0(2)];
end