generate_points <-
function(W, unlabeled, labeling) {
	n <- nrow(W);
	N <- length(labeling[-unlabeled]);
	N_pos <- length(which(labeling[-unlabeled] == 1));
	pos_rate <- N_pos/N;					
	p <- length(unlabeled);
	# generating random labels for unlabeled nodes according to the distribution B(p, N_pos/N) 			
	y_rand <- generate_labels(p, pos_rate);		
	labeling[unlabeled] <- y_rand;			
	labeling[labeling <=0] <- 0;
	pos_vect <- W %*% labeling;	
	neg_vect <- W %*% (1-labeling)
	res <- list(pos_vect = pos_vect, neg_vect = neg_vect);
	return (res);
}