function C = CreateC(E)
C=zeros(8,8);
C(1,1)=sqrt(E(1,1));
for j=2:8
    for k=1:j-1
        for l=1:k-1
            C(j,k)=C(j,k)+C(j,l)*C(k,l);
		end
        C(j,k)=(E(j,k)-C(j,k))/C(k,k);
    end
    C(j,j)=0;
    for l=1:j-1
        C(j,j)=C(j,j)+C(j,l)*C(j,l);
    end
    C(j,j)=sqrt(E(j,j)-C(j,j));
end