for i = 1:nA
i
    S = 0;
    for k = 1:i-1
        S = S + L(i,k)^2;
    end %k
    L(i,i) = sqrt(A1(i,i) - S);
     
    for jj = i+1:nA
        S = 0;
        for k = 1:i-1
            S = S + L(jj,k) * L(i,k);
        end %k 
        L(jj,i) = (A1(jj,i) - S) / L(i,i);
    end %jj
    
    for b = 1:nE
        S = 0;
        for k = 1:i-1
            S = S + L(i,k) * V(k,b);
        end %k
        V(i,b) = (B(i,b) - S) / L(i,i);
    end %b
end %i