Метод Монте Карло (неправильные значения)

#include "iostream"
#include "cmath"
#include <time.h>
using namespace std;
 
//метод прямоугольников, трапеций, Монте-Карло
 
double Func(double x) 
{
    return pow(x,2)*sin(x);
}
 
double Method_Pr(double a, double b, int n) 
{
    int i;
    double result, h;
    result = 0;
    h = (b - a) / n; 
    for (i = 0; i < n; i++)
    {
        result += Func(a + h * (i + 0.5)); 
    }
    result *= h;
    return result;
}
 
double Method_Trap(double a, double b, int n)
{
    double result, h;
    result = 0;
    h = (b - a) / n;
    double the_result = (Func(a) + Func(b)) / 2;
    while (a<b-h)
    {   
        a = a + h;
        result += Func(a);
    }
    result = h*(the_result + result);
    return result;
}
 
 
double Method_Monte_Carlo(double a, double b, int n)
{
    srand(time(NULL));
    double t, x, z;
    double result;
    int i, HIT=0;
    double minz,maxz;
    maxz = Func(a);
    minz = Func(a);
 
    for (x=a; x <= b; x+=(b-a)/n) 
    {
        if (Func(x) > maxz)
            maxz = Func(x);
        if (Func(x) < minz)
            minz = Func(x);
    }
    for (i = 0; i < n; i++) 
    {
        x = (double)a+(b-a)*rand();
        z = (double)(minz+(maxz-minz))*rand();
        
        if (z < Func(x))
            HIT++;
    }
 
    result = (double)HIT*(maxz-minz)*(b-a) / n;
    return result;
}
 
int main()
{
    double integral;
    integral = Method_Pr(1, 6, 100);
    printf("The value of the integral(pr) is: %lf \n", integral);
    integral = Method_Trap(1, 6, 1000);
    printf("The value of the integral(trap) is: %lf \n", integral);
    integral = Method_Monte_Carlo(1, 6, 1000);
    printf("The value of the integral(mc) is: %lf \n", integral);   
    system("pause");
    return 0;
}

x = (double)a+(b-a)*rand()/RAND_MAX;
z = (double)(minz+(maxz-minz))*rand()/RAND_MAX

double Method_Monte_Carlo(double a, double b, int n)
{
    srand(time(NULL));
    double maxz=Func(a);
    double minz=maxz;
    for (double x=a; x <= b; x+=(b-a)/n) 
    {
        double y=Func(x);
        if (y > maxz)
            maxz = y;
        if y < minz)
            minz = y;
    }
   const int N=100000; // N - это количество испытаний
   int HIT=0;
   for (i = 0; i < N; i++) 
   {
        x = a+(b-a)*rand()/RAND_MAX;
        z = minz+(maxz-minz)*rand()/RAND_MAX;
        double y=Func(x);
        if (y>0 && z<y)
            HIT++;
       else if(y<0 && z>y)
           HIT--;
    }
 
    result = (maxz-minz)*(b-a)*HIT/ N;
    return result;
}