function J=lvjac2a(x,y,b)

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


end