[CPP] 수치해석 : Rectangle Method 구분구적분 (적분 프로그래밍)

#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;

}