function J=lvjac(x,y)

J=[-4*x+1-y   -x ;...
   -y         -4*y+3-x];


end