Форум программистов, компьютерный форум CyberForum.ru

Интегрирование методом трапеций - C++

Восстановить пароль Регистрация
 
Остап Бендер
12 / 12 / 2
Регистрация: 17.11.2010
Сообщений: 143
24.09.2011, 12:08     Интегрирование методом трапеций #1
Доброго времени суток.
Мне нужно реализовать интегрирование методом трапеций.
Формулу я взял в Википедии
http://www.cyberforum.ru/cgi-bin/latex.cgi?\int_{a}^{b}f(x)=(\frac{f(a)+f(b)}{2}+\sum_{i=1}^{k-1}f({x}_{i}))h

и вот, что я написал:
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
typedef double (*real)(double x);
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
    : TForm(Owner)
{
}
//---------------------------------------------------------------------------
double f(double x)
{
    return 2*x;
}
//---------------------------------------------------------------------------
double trap(double a, double b, int k, real f)
{
    int i;
    double summa=0;
    double x[100];
    double step = abs(a-b)/k;
    x[0]=a;
    x[k]=b;
    for (i = 0; i < k; i++)
    {
        x[i]=x[i-1]+step;
    }
    for (i = 1; i < k-1; i++)
    {
        summa=summa+f(x[i]);
    }
    return ((f(a)+f(b))/2 + summa)*step;
}
//---------------------------------------------------------------------------
void __fastcall TForm1::BitBtn1Click(TObject *Sender)
{
 double a,b,k;
 a = StrToFloat(Edit1->Text);
 b = StrToFloat(Edit2->Text);
 k = StrToFloat(Edit3->Text);
 Edit4->Text = FloatToStr(trap(a,b,k,*f));
}
//---------------------------------------------------------------------------
Скажите, пожалуйста, все что тут не так (правильно считает, только когда k==1) и почему формула именно такова, ведь тут все, насколько я понял, основывается на сумме площадей трапеций. Почему не так:
http://www.cyberforum.ru/cgi-bin/latex.cgi?\sum_{i=1}^{k}\frac{f({x}_{i})+f({x}_{i-1})}{2}({x}_{i}-{x}_{i-1})
?

Спасибо за внимание.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
24.09.2011, 12:08     Интегрирование методом трапеций
Посмотрите здесь:

C++ Решение интеграла методом трапеций
Решение интеграла методом трапеций C++
Интегрирование методом прямоугольников и трапеций C++
Вычисление интеграла методом трапеций C++
C++ Численное интегрирование (метод трапеций). Как гарантировать нужную точность?
Решение интеграла методом трапеций? C++
C++ Вычисление интеграла методом трапеций
Определенный интеграл методом трапеций C++

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

Текущее время: 20:54. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru