24.09.2011, 12:08. Просмотров 508. Ответов 0
Доброго времени суток.
Мне нужно реализовать интегрирование методом трапеций.
Формулу я взял в Википедии
и вот, что я написал:
| 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) и почему формула именно такова, ведь тут все, насколько я понял, основывается на сумме площадей трапеций. Почему не так:
?
Спасибо за внимание.
