%% Get cdf results from data file of simulations
function result=getsimu(filename)
    fin = fopen(filename, 'r');
    status = 0;
    result = [];
    while(~feof(fin))
        str = fgetl(fin);
        if(status==0)
            Trhop = strread(str);  % line of T, rho, p
            str = fgetl(fin);
            cds = strread(str);% line of cluster distribution function
            result = [Trhop, cds];
            status=1;
        elseif(status==1)
            Trhop = strread(str);
            status=2;
        else
            cds = strread(str);
            result = [result; Trhop, cds];
            status=1;
        end
    end
    fclose(fin);
end