C/C++ - Найти интеграл методами правых прямоугольников и Симпсона

Внимание, важное сообщение: Дорогие Друзья!
В ноябре далекого 2001 года мы решили создать сайт и форум, которые смогут помочь как начинающим, так и продвинутым пользователям разобраться в операционных системах. В 2004-2006г наш проект был одним из самых крупных ИТ ресурсов в рунете, на пике нас посещало более 300 000 человек в день! Наша документация по службам Windows и автоматической установке помогла огромному количеству пользователей и сисадминов. Мы с уверенностью можем сказать, что внесли большой вклад в развитие ИТ сообщества рунета.

Но... время меняются, приоритеты тоже. И, к сожалению, пришло время сказать До встречи! После долгих дискуссий было принято решение закрыть наш проект. 1 августа форум переводится в режим Только чтение, а в начале сентября мы переведем рубильник в положение Выключен

Огромное спасибо за эти 24 года, это было незабываемое приключение. Сказать спасибо и поделиться своей историей можно в данной теме.

С уважением, ваш призрачный админ, BigMac...

C/C++ - Найти интеграл методами правых прямоугольников и Симпсона

kadett46

Пользователь

Сообщения: 58
Благодарности: 0

Профиль | Отправить PM | Цитировать

Составить программу вычисления определенного интеграла методами правых прямоугольников и Симпсона с помощью вызова двух различных функций и оценить погрешность.

СРОЧНО! Заранее благодарю.
Вот что-то есть, не знаю правильно ли

Код:

#include <vcl.h>
#include <iostream>
#include <conio.h>
#include <math.h>

using namespace std;

double f(double x, double dx){

    return ((exp(x)*dx)/x);
}

double f1(double x, double dx1){

    return ((exp(x)*dx1)/x);
}

int main(){

    const int N = 1000; //число точок
    double pogr;
    
    double a, b;
    do{
        cout << "Enter the interval of integration through the gap: ";
        cin >> a >> b;
    }while(a >= b);

    double dx1 = (b - a) / (N - 1);
    double sum1 = 0;
    for (double x = a; x < (b + dx1 / 2); x += dx1)
        sum1 += dx1 * f1(x,dx1);

    double dx = (b - a) / (N - 1);
    double sum = 0;
    for (double x = a; x < (b + dx / 2); x += dx)
        sum += dx * f(x,dx);
    cout << "Method rectangles: " << sum << endl;

    sum = 0;
    for (int i = 1; i < N - 1; ++i){
        double x = a + i * dx;
        if (i % 2 == 0)
            sum += 2 * f(x,dx);
        else
            sum += 4 * f(x,dx);
    }

    pogr=((dx1 / 3) * (a + sum + b))-(sum1);

    cout << "Simpson's method: " << (dx / 3) * (a + sum + b) << endl;

    cout << "Error: " << pogr << endl;

    //cout << "The exact value: " << b * b * b - a * a * a << endl;

    getch();
}

Отправлено: 19:16, 11-12-2014

lxa85

Необычный

Сообщения: 4466
Благодарности: 995

Профиль | Сайт | Отправить PM | Цитировать

Показательная интегральная функция...
раз
два
три
четыре
пять

-------
- Я не разрешаю тебе быть плохой! Потому что плохие люди совершают плохие поступки. А это нехорошо!
(Из наставлений 5 летней девочки своей младшей сестре)

Отправлено: 03:19, 12-12-2014 | #2