 	int xdim = 3; /* data dimension */
 
 	/* solution - for generating the problem */
-	int solution_K = 4;
+	int solution_K = 3;
 	int solution_xmem = 0;
 	
-	double solution_theta[12] = {
+	double solution_theta[9] = {
 		 0.0, 0.0, 0.0,		/* K=1,n=1,2,3:  mu */
 		 1.0, 0.0, 0.0,		/* K=2,n=1,2,3 */
-		 0.0, 1.0, 0.0,		/* K=3,n=1,2,3 */
-		 0.5, 0.5, 0.2		/* K=4,n=1,2,3 */
+		 0.0, 1.0, 0.0		/* K=3,n=1,2,3 */
 	};
 
-	double solution_xstart[4] = {
-		0.0, 0.0, 0.0, 0.0	/* n=1,2,3,4 */
+	double solution_xstart[3] = {
+		0.0, 0.0, 0.0	/* n=1,2,3 */
 	};
 
-	double noise_covariance[12] = {
+	double noise_covariance[9] = {
 		0.05, 0.05, 0.01,  
 		0.05, 0.05, 0.01,  
-		0.05, 0.05, 0.01,  
-		0.05,  0.05, 0.05
+		0.05, 0.05, 0.01  
 	};
 
 	/* model parameters */

int solution_get_cluster_id(int t, int T){
	int id_cluster = 0;
	double coeff = (double)T/4.0;
	if(t >= coeff && t < 2*coeff){
		id_cluster = 1;
	}
	if(t >= 2*coeff && t < 3*coeff){
		id_cluster = 2;
	}
	if(t >= 3*coeff){
		id_cluster = 3;
	}


	int *xmem = new int(num);

	int i;
	for(i=0; i<num; i++ ){
		K[i] = Kforall;
		epssqr[i] = epssqrforall;
		xmem[i] = 0;
