function [fun,jac] = f(x,par)
fun = [1-par(1)*x(1)^2+x(2);par(2)*x(1)];
jac = [-2*par(1)*x(1), 1;par(2), 0];


