// Generate four groups of data
k=0;
for (i=0;i<4;i++) {
for (j=0;j<(int)(NPT*frac[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++;
}