Процедура вычисления кратного интеграла методом Монте-Карло

#include <stdio.h>
#include <time.h>
#include <math.h>
#include <conio.h>
#include <stdlib.h>
 
double func(double a, double b, double c)
{
    return (a*b*c)/sqrt(pow(a,2)+pow(b,2)+pow(c,2));
}
 
int debug(int x)
{
    if((x>-32768) && (x<32767));
    else
   {
    printf("\n \t \t Vy vveli nevernyj simvol. \n \t \t Neobhodimo bylo vvesti celoe chislo. \n \n");
    getch();
    exit(1);
    }
    return 0;
}
 
void main()
{
    double x, y, z;
   double a, b, c;
   unsigned int n;
   int i;
   double sum=0;
   double S=0;
   double min=0;
 
   srand(time(NULL));
    printf("\n \t Vvediteniz. granicu integrirovanija (x): ");
    scanf("%lf",&x);
   debug(x);
 
    printf("\t Vvedite granicu integrirovanija (y): ");
    scanf("%lf",&y);
   debug(y);
 
   printf("\t Vvedite  granicu integrirovanija (z): ");
    scanf("%lf",&z);
   debug(z);
 
 
    inv: printf("\n \t Vvedite chislio ispytanij \n");
    printf("\t N = ");
    scanf("%d",&n);
   if((n<0) || (n>65535))
   {
    printf("\n \t Vy vveli nepravilnoe chislo, kolichestvo mozhet byt ot 0 do 65535 \n");
    goto inv;
   }
 
 
    printf("\n \t / / / \n");
    printf("\t | | | (x*y*z)/sqr(pow(x,2)+pow(y,2)+pow(z,2))dxdydz \n");
    printf("\t / / / \n");
   printf("\t   D   \n");
 
   // вычисления двойного интеграла
   for(i=0;i<n;i++)
   {
 
    a=x*rand()/RAND_MAX;//генерируем x в интервале
      b=y*rand()/RAND_MAX;//генерируем y в интервале
      c=z*rand()/RAND_MAX;//генерируем z в интервале
 
      if (sqrt(pow(a,2)+pow(b,2)+pow(c,2))<=1)
    {
           sum+=func(a, b, c);
        }
         printf("\n \n \t a,b,c = %8.3lf %8.3lf %8.3lf\n", a,b,c);
   }
   S=(((x-0)*(y-0)*(z-0))/n)*sum;
   printf("\n \n \t Integral = %8.3lf \n", S);
 
   getch();
}