restart;
eqns:=[x[1]^3+x[1]^2*x[2]-x[1]*x[3]+6,
       exp(x[1])+exp(x[2])-x[3],
       x[2]^2-2*x[1]*x[3]-4];
with(linalg):
F:=evalf(vector(eqns));
X:=[seq(x[i],i=1..3)];
J:=jacobian(F,X);
with(CodeGeneration);
C(F);
C(J);