Решение определенного интеграла методом Монте Карло

#pragma hdrstop
#include <stdio.h>
#include <time.h>
#include <math.h>
#include <stdlib.h>
 
int main(int argc, char* argv[])
{
 double s=0,x,Integral, a, b;
 int n,i;
 srand(time(NULL));
 printf("a:\n");  //нижний предел
 scanf("%lf",&a);
 printf("b:\n");  //верхний предел
 scanf("%lf",&b);
 printf("n:\n");  //число испытаний
 scanf("%d",&n);
 for (i=1; i<=n; i++) {
    x=a+(b-a)*rand()/RAND_MAX;
    s=sin(2*x);
}
 Integral=((b-a)*s)/n ;
 printf("Integral=%f\n",Integral);
 return 0;
}

#include <iostream>
#include <random>
#include <cmath>
///////////////////////////////////////////////////////////////////////////////
bool    val_in_segment
    (
        double  val,
        double  L,
        double  R
    )
{
    return      L       <=  val
            &&  val     <=  R;
}
///////////////////////////////////////////////////////////////////////////////
double  get_monte_carlo_definite_integral_sin_2x
    (
        double  x_min,
        double  x_max,
        int     n
    )
{
    static  std::default_random_engine          gen;
    std::uniform_real_distribution< double >    x_distr         ( x_min,    x_max );
    const   double                              Y_MIN           { -1    };
    const   double                              Y_MAX           { 1     };
    std::uniform_real_distribution< double >    y_distr         ( Y_MIN,    Y_MAX );
    int                                         hits_counter    {};
 
    for( int  i{}; i < n; ++i )
    {
        double  x       =   x_distr     ( gen );
        double  y       =   y_distr     ( gen );
        double  f_x     =   sin         (2 * x);
 
        bool    is_successful_hit   =   f_x     >   0
                                            ?   val_in_segment( y,  0,      f_x )
                                            :   val_in_segment( y,  f_x,    0   );
 
        if( is_successful_hit )
        {
            f_x     >   0
                ?   ++hits_counter
                :   --hits_counter;
        }//if
    }//for
 
    double  S   =       ( x_max     -   x_min )
                    *   ( Y_MAX     -   Y_MIN );
 
    return      hits_counter
            *   S
            /   n;
}
///////////////////////////////////////////////////////////////////////////////
int     main()
{
    for(;;)
    {
        std::cout   <<  std::endl
                    <<  "a = ";
 
        double  a{};
        std::cin    >>  a;
 
        std::cout   <<  "b = ";
        double  b{};
        std::cin    >>  b;
 
        std::cout   <<  "n = ";
        double  n{};
        std::cin    >>  n;
 
        std::cout   <<  get_monte_carlo_definite_integral_sin_2x
                            (
                                a,
                                b,
                                n
                            )
 
                    <<  std::endl;
    }//for
}