function JF=fnlb9p5j(x);
JF=[-x(2)-x(1).^2-2*x(1).*(3+x(1)),3-x(1);2*(2+x(1))./(1+(2+x(1)).^2).^2,0];
