/*
FOR I THRU 4 DO FOR J THRU 4 DO (
N[I,J]:RATSIMP(
SUM(SUM(UE[K,L]*DIFF(DIFF(LG[I,J],OMEGA[K]),OMEGA[L]),K,1,4),L,1,4)-
SUM(SUM(SUM(SUM(
UE[A,B]*UG[L,P]*DIFF(LG[I,L],OMEGA[A])*DIFF(LG[J,P],OMEGA[B]),A,1,4),B,1,4),P,1,4),L,1,4)));
*/
/* IF DIAGMETRIC */

FOR I THRU 4 DO (
NN[I,I]:RATSIMP(
SUM(UE[K,K]*DIFF(LG[I,I],OMEGA[K],2),K,1,4)-
SUM(SUM(UE[A,A]*UG[P,P]*DIFF(LG[I,P],OMEGA[A])^2,A,1,4),P,1,4)));