// Generate four groups of data
k=0;
for (i=0;i<4;i++) {
count[i]=(int)(NPT*frac[i]);
for (j=0;j<count[i];j++) {
d0=sigma[i][0]*ndev.dev();
d1=sigma[i][1]*ndev.dev();
x[k][0]=means[i][0]+d0*vec1[i][0]+d1*vec2[i][0];
x[k][1]=means[i][1]+d0*vec1[i][1]+d1*vec2[i][1];
k++;
}