x[k][1]=means[i][1]+d0*vec1[i][1]+d1*vec2[i][1];
k++;
}
}
Gaumixmod gmix = new Gaumixmod(x,guess);
for (i=0;i<100;i++) {
flag=gmix.estep();
if (flag < 1.e-6) break;
gmix.mstep();
}
// check for convergence
// System.out.println(" flag: %f\n", flag);
localflag = flag > 1.e-6;