{процедура для розрахунку інтеграла за квадратурної
формулою Чебишева}
procedure cheb (var y: aa; var ich: real);
var i: integer;
Begin
ich: = 0;
for i: = 1 to n do
ich: = ich + y [i] * h;
end;
{процедура виведення таблиці}
procedure tabl;
var i: integer;
Begin
writeln (' ___________________________________ '); p> writeln ('| i | t | x | y | ');
writeln (' ___________________________________ '); p> writeln ('| 1 | ', K: 9:6,' | ', x [1]: 9:6,' | ', y [1]: 9:6,' | ');
writeln ('| 2 |', l: 9:6, '|', x [2]: 9:6, ' | ', Y [2]: 9:6,' | ');
writeln ('| 3 | ', Z: 9:6,' | ', x [3]: 9:6,' | ', y [3]: 9:6,' | ');
writeln ('| 4 | ', L: 9:6,' | ', x [4]: ​​9:6,' | ', y [4]: ​​9:6,' | ');
writeln ('| 5 | ', K: 9:6,' | ', x [5]: 9:6,' | ', y [5]: 9:6,' | ');
writeln (' ___________________________________ '); p> end;
Begin
clrscr;
writeln ('П Р О Г Р А М М А Д Л Я В И Ч І З Л Е Н І Я ');
writeln ('О П Р Е Д Е Л Е Н Н О Г О І Н Т Е Г Р А Л А ');
writeln;
writeln ('Введіть кордону інтегрування a, b: ');
readln (a, b);
vvod (a, b, x);
h: = (b-a)/n;
writeln ('h =', h: 9:6);
form (x, y);
cheb (y, ich);
tabl;
writeln ('I =', ich: 8:6);
end.
Висновок результату:
В
П Р О Г Р А М М А Д Л Я В И Ч І З Л Е Н І Я
О П Р Е Д Е Л Е Н Н О Г О І Н Т Е Г Р А Л А
Введіть кордону інтегрування a, b:
0 1.5708
h = 0.314160
____________________________
| i | t | x | y |
____________________________
| 1 | -0.832498 | 0.131556 | 0.131177 |
| 2 | -0.374541 | 0.491235 | 0.471716 |
| 3 | 0.000000 | 0.785400 | 0.707108 |
| 4 | -0.374541 | 0.508765 | 0.487099 |
| 5 | -0.832498 | 0.868444 | 0.763325 |
____________________________
I = 0.804383
