function Jf = fcnJCP5p19(x,gamma,delta)
Jf = [(gamma-1)*x(2)-1 (gamma-1)*x(1); -x(2) -x(1)-2*x(2)+delta-1];
