генерации случайной величины

#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <conio.h>
#pragma hdrstop
void print_mas(double *a[], int n)
{
   printf("\текст массива \n");
   for(int i = 0; i < n; i++)
   {
      printf("m[x] = %.4lf\tm[y] = %.4lf\t",a[i][0], a[i][1]);
      if((i+1) % 2 == 0)
       printf("\n");
   }
   printf("\n");
}
int main(int argc, char* argv[])
{
  double**m,x,y,s,kx,ky,gj,k;
  int n = 0;
  printf("input count mass: "); scanf("%i", &n);
  if(n <= 0)
  {
    printf("error! n <= 0");
    return 0;
  }
  m = new double* [n];
  randomize();
  for(int i=0; i<n; i++)
  {
      m[i]= new double[2];         
      do
      {
        x  = (rand() % 100) / 1000.0;
        y  = (rand() % 100) / 1000.0;
        s = x*x + y*y;
      }
      while((s == 0) || (s > 1));
      m[i][0] = x * sqrt((-2 * log(s))/s)*0.1;
      m[i][1] = y * sqrt((-2 * log(s))/s)*0.1;
  }
  print_mas(m,n);
   for(int i = 0; i < n; i++)
    {
    kx+=m[i][0];
    ky-=m[i][1];
    }
    kx/=n;
    ky/=n;
    for(int i = 0; i < n; i++)
      {
      k+=(m[i][0]-kx)/(m[i][1]-ky) ;
      }
     k/=n;
    for(int i = 0; i < n; i++)
      {
      gj+=(m[i][0]-kx)*(m[i][0]-kx);
      }
      gj/=n;
  for(int i = 0; i < n; i++)
    delete m[i];
  delete []m;
  printf("k = %.4lf\t",k);
  printf("\nM(x) = %.4lf\t",kx);
  printf("\nM(y) = %.4lf\t",ky);
  printf("\n6 = %.4lf\t",gj);
  printf("\nend work program. Press any key...");
  getch();
  return 0;
}