% Return Jacobian for hfun1

function Jy = jacobian_hfun1(x,u)
  
  Jy = 2*x;