h:x1^2+x2^2;

dh:jacobian([h],[x1,x2,x3]);

D:nullspace(dh);