function g = grad_fex1(x)
% gradient of fex1
g=[2*x(1); 4*x(2)];
