#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#define N 32

typedef struct particle
{
	float r[N][3];   //position
	float v[N][3];   //velocity
}Particle;

float sq(float a)
{
	return a*a;
}

void init_r(Particle *p) 
{

	int i,j;
	float r2;
	for(i=0;i<N;){
		p->r[i][0]=2*drand48()-1.0;
		p->r[i][1]=2*drand48()-1.0;
		p->r[i][2]=2*drand48()-1.0;
		r2=sqrtf(sq(p->r[i][0]) + sq(p->r[i][1]) + sq(p->r[i][2]));
	
		p->v[i][0]=0.0;
		p->v[i][1]=0.0;
		p->v[i][2]=0.0;
		if (r2 <= 1.0) i++;
	}
}

int main(int argc,char** argv)
{
	int i;
	Particle p;
	init_r(&p);
	for(i=0;i<N;i++)
	{
		printf("%f,%f,%f \n",p.r[i][0],p.r[i][1],p.r[i][2]);
	}
	return 0;
}

-0.815405,-0.025566,0.053501 
-0.091133,-0.533643,0.662584 
0.863463,0.136119,0.112189 
-0.495280,-0.403606,0.751962 
0.063114,0.840522,0.030862 
0.248813,-0.747075,0.235100 
0.541857,-0.626928,-0.006557 
0.019001,-0.368122,0.403405 
-0.265507,0.752625,-0.258011 
0.050505,-0.144331,0.085486 
0.335957,-0.657716,0.661558 
0.620317,-0.071578,0.138370 
0.024406,0.173021,0.477517 
-0.332941,-0.214344,-0.646145 
-0.250082,-0.946428,-0.070469 
0.570686,-0.115530,-0.424192 
0.706706,0.051944,0.456833 
-0.467371,0.566968,-0.063912 
-0.606005,-0.105646,0.456880 
0.599943,-0.713566,0.220921 
-0.673943,-0.150365,-0.040292 
-0.382420,0.885252,-0.223951 
-0.161893,-0.407260,-0.243535 
-0.162567,-0.723104,-0.238307 
-0.010950,0.765044,0.547660 
-0.316169,-0.086376,0.176243 
-0.193679,0.170966,-0.617678 
0.363142,0.623592,0.284430 
-0.341315,0.475107,0.646965 
0.442522,0.315965,-0.282451 
0.220730,-0.169993,0.629908 
0.202082,0.028182,0.260429