#include <stdlib.h>
#define RND ((double)rand()/RAND_MAX)

void pseudo_random_vector(int dim,double* a,double* b, double* x) {

	for(int i=0;i<dim;i++) x[i]=a[i]+RND*(b[i]-a[i]);

	}