Метод прямоугольника

#include "stdafx.h"
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <conio.h>
#include <math.h>
 
#define n_MAX 100000L
 
double Karlo();
double MRect();
 
//////////////////////////////////////////////////////////////////////
double f1(double x){return 1+sqrt(1-x*x);}
double f2(double x){return 1-sqrt(1-x*x);}
 
//////////////////////////////////////////////////////////////////////
double main(int argc, char* argv[])
{
//  printf("Hello World!\n");
//  printf("\n%x",RAND_MAX);
//  getch();
 
    double z,s;
    z=Karlo();
    printf("\nz=%f",z);
    s=MRect();
    printf("\ns=%f",s);
    printf("\n");
    return 0;
}
 
//////////////////////////////////////////////////////////////////////
 
double Karlo(){
    long k,n;
    long t;
    double x,y,s;
    double a,b,c,d;
    double r=1; //radiys
    //------------------
    a=-1; b=1; c=0; d=2;
    srand(time(&t)%37);
    n=k=0;
    while (n<n_MAX){
        x=(double)rand()/RAND_MAX;
        y=(double)rand()/RAND_MAX;
        x=(b-a)*x+a;
        y=(d-c)*y+c;
        if (((x*x+pow(y-1,2)<=1) && (x>0) && (y<exp(-x))) ||
            ((x*x+y*y<=1) && (x<0) && (y>exp(-x)))) k++;
        n++;
    }
    s=(double)k/n*(b-a)*(d-c);
    return s;
}
/////////////////////////////////////////////////////////////////////////
double MRect(){
    double a,a2,b,b2,x,dx,dy,s1,s2;
    int i,n;
    n=10000;
    a=-1; b=0;
    dx=(b-a)/n;
    s1=0; x=a; i=0;
    while (i<n) {
        dy=f1(x)-exp(-x);
        if (dy>0) {s1=s1+dy*x;}
        x=x+dx; i++;
    }
    a2=0; b2=1;
    dx=(b-a)/n;
    s2=0; x=a2; i=0;
    while (i<n) {
        dy=f2(x)-exp(-x);
        if (dy>0) {s2=s2+dy*x;}
        x=x+dx; i++;
    }
    return s1+s2;
 
}