Решить методом трапеций или Симпсона определенный интеграл

@gika111 · Регистрация: 04.01.2010

Здавствуйте! Помогите пожалуйста!
Решить методом трапеций или Симпсона определенный интеграл
Заранее спасибо!

@Maxak · 16.07.2013, 19:44

как то так по Сипсону наверное. FormulaValue () - функция которую хотите интегрировать. а, b - начало и конец отрезка, eps - точность относительная вычислений.

C++

                           int N = 0;
        
                           double h = DBL_MAX;
        double eps = 0.01;
                           double a, b; 
                            a = 1;
                           b = 2;   
        double S1 = DBL_MAX;
        double S2 = DBL_MAX;
        do
        {
            
            S2 = S1;
            S1 = 0;
            N+=2;
            h = (a - b)/N;
            double x = a + h;
            for(int n = 1; n< N; n++)
            {
                n++;
                double x1 = (x-h);
                double x2 = (x+h);
                S1 += FormulaValue(x1) + 4*FormulaValue(x)+FormulaValue(x2);
                x +=2.*h;
            }
            S1 = S1*h/3.;
        }
        while((2*abs(S1-S2)/abs(S1+S2)) > eps);
        
             dblVal = (S1+S2)/2.;
    
    return dblVal;

Добавлено через 50 минут
S1 и S2 лучше инициализировать нулями, а то не работает=))

Опции темы

Решить методом трапеций или Симпсона определенный интеграл

Решение