一维蒙特卡洛积分（C++实现）

#include<iostream>
#include"math.h"
using namespace std;

double suiji()
{
    return  fabs(2.0 * ( (double) rand() / (double) RAND_MAX ) - 1.0 );
}

int main()
{
double result;
int m;
double a,b;
cout<<"a=";
cin>>a;
cout<<"b=";
cin>>b;
cout<<"m=";
cin>>m;

double sum=0.0;
for(int i=0;i<m;i++)
{
double s=a+(b-a)*suiji();
sum+=s;
}
result=((b-a)/m)*sum;
cout<<result<<endl;
}