function [f J]=NRP_1(x)
 f(1) = 3*x(1)-cos(x(2)*x(3))-0.5;
 f(2) = x(1)^2-81*(x(2)+0.1)^2+sin(x(3))+1.06;
 f(3) = exp((-x(1)*x(2)))+20*x(3)+(10*pi-3)/3;
 if nargout<2, return, end
 J(1,1) = 3.0; J(1,2) = sin(x(2)*x(3))*x(3); J(1,3) = sin(x(2)*x(3))*x(2);
 J(2,1) = 2.0*x(1); J(2,2) = -162.0*(x(2)+0.1); J(2,3) = cos(x(3));
 J(3,1) = -exp((-x(1)*x(2)))*x(2); J(3,2) = -exp((-x(1)*x(2)))*x(1);
 J(3,3) = 20.0;
end