function J = DG(X)
u = X(1); v = X(2);

j_11 = -3 * u^2; j_12 = 1; j_21 = 2 * u; j_22 = 2 * v;

J = [j_11, j_12; j_21, j_22];