A=[1,1,1];
L=[1,2,2];
H=[1.003,0.501,-0.503];
w=7;
P=diag(2./L);
P1=inv(P);
N=A*P1*A';
Ka=7/N*8;
V=P1*A'*Ka;
h=V'+H