N = [5.2,2.4,0,-1.7;2.4,7.4,2.3,2.7;0,2.3,6.3,0;-1.7,2.7,0,4.4];
D = det(N)
DJ = diag(N)'