function J = linode_DFDX_aut(x, p) 

x1 = x(1,:);
p1 = p(1,:);

J = zeros(4,4,numel(x1));
J(1,2,:) = 1;
J(2,1,:) = -p1;
J(2,2,:) = -1;
J(2,4,:) = 1;
J(3,4,:) = 1;
J(4,3,:) = -1;
end
