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