computeCPCModel(covList, m_constants)

function [cpcCovList, LList, B, FList] = computeCPCModel(covList, m_constants)

[converge, B, FList] = falg(covList, m_constants.numVar, m_constants.numGrp);

if converge ~=0
    disp(sprintf('\n!!F-G algorithm did not converge!!\n No CPC Model is possible!'))
    BList = {};
    LList = {};
    cpcCovList = {};
end;

if converge == 0
    for i = 1:m_constants.numGrp 
        LList{i} = diag(FList{i});
        cpcCovList{i} = B * diag(LList{i}) * B';
    end
end