#include <iostream>
#include <math.h>
#define PI 3.141592
using namespace std;
/**********************
**** Function ****
***********************/
void Test_Function(void);
double Rectangle_Method(double First, double Second, int Number);
/**********************
**** Main ****
***********************/
int main(void)
{
Test_Function();
return 0;
}
void Test_Function(void)
{
cout <<"**************** TEST ****************"<<endl;
cout <<"Value of Intergral 0 to 10 -> (x^5-x^4)dX "<<endl;
cout <<"Numerical Methods (TEST) : "
<< Rectangle_Method_for_Test(0, 10, 10000000) <<endl;
cout <<"True : "<<146666.66666<<endl;
cout <<"**************** TEST ****************"<<endl;
cout <<endl;
}
double Rectangle_Method_for_Test(double First, double Second, int Number) // 처음 값, 마지막 값, 등분 수
{
double Sum=0.0;
double *Division=new double[Number+1];
double x, y, dX;
x = (Second - First) / Number;
for(int iPoint=0; iPoint<=Number; iPoint++)
{
Division[iPoint]=x*iPoint;
if(iPoint==Number)
cout <<"Last : Division["<<iPoint<<"] : "
<<Division[iPoint]<<endl;
}
for(int iActive=0; iActive<Number; iActive++)
{
dX=Division[iActive]+(Division[iActive+1]-Division[iActive])/2;
y=(dX*dX*dX*dX*dX - dX*dX*dX*dX);
Sum=Sum+(x*y);
}
delete []Division;
return Sum;
}