function [initialPop]=InitializePop(n,p)
for i=1:p
    initialPop(i,1)=rand(1);
    sumx=0;
    for j=2:n
        initialPop(i,j)=rand(1);
        sumx=sumx+initialPop(i,j);
    end
    initialPop(i,n+1)=initialPop(i,1);
    gx=1+9*(sumx/(n-1));
    initialPop(i,n+2)=gx*(1-(initialPop(i,n+1)/gx)^2);
end