/b> - Використовуючи масиви x i і y i , вираховує по квадратурної формулою Чебишева наближене значення інтеграла.
Процедура TABL - Це підпрограма, що здійснює висновок таблиці вузлів (аргумент - функція)
При запуску програми потрібно ввести кордону інтегрування.
Після введення кордонів інтегрування використовується процедура VVOD, а потім вираховується і виводитися на екран крок табулювання функції h.
Після цього використовуємо процедури FORM і CHEB.
Отримавши результат, виводимо таблицю (процедура TABL) і інтеграл.
В
В
4. Висновок і висновки.
Таким чином очевидно, що при обчисленні визначених інтегралів за допомогою квадратурних формул, а зокрема за формулою Чебишева не дає нам точного значення, а тільки наближене.
Щоб максимально наблизитися до достовірного значенням інтеграла треба вміти правильно вибрати метод і формулу, за якою вестиметься розрахунок. Так само дуже важливо те, який буде взятий крок інтегрування.
Хоча чисельні методи і не дають дуже точного значення інтеграла, але вони дуже важливі, оскільки не завжди можна вирішити завдання інтегрування аналітичним способом.
В В
В
5.Спісок літератури:
1. Ракітін Т.А., Первушин В.А. "Практичне керівництво з чисельних методів з додатком програм мовою Basic "
2. Крилов В.І. "Наближені обчислення інтегралів" - М.: Фізмат.
3. Демидович і Марон "Основи обчислювальної математики"
4. Копченова і Марон "Обчислювальна математика в прикладах і задачах "
5. Вольвач О.М., Крисевич В.С. Програмування на мові Паскаль для ПЕОМ ЄС. Мінськ.: 1989
6. Зуєв Є.А. Мова програмування Turbo Pascal. М.1992 р.
7. Скляров В.А. Знайомтеся: Паскаль. М. 1988
В
6. Лістинг програми.
Програма написана мовою Tubro Pascal 7.0 для MS-DOS. Нижче наведено її лістинг:
program integral;
uses crt;
const n = 5;
k = -0.832498;
l = -0.374541;
z = 0.0;
type aa = array [1 .. n] of real;
var x, y: aa;
a, b, h, ich: real;
{заповнення х-сов в масив х [5]}
procedure vvod (var a, b: real; var c: aa);
var i: integer;
t: aa;
Begin
t [1]: = k;
t [2]: = l;
t [3]: = z;
t [4]: ​​= l;
t [5]: = k;
for i: = 1 to n-1 do
c [i]: = ((b + a)/2 + (b-a)/2 * t [i]);
for i: = n-1 to n do
c [i]: = 1 - c [n +1- i];
end;
{заповнення y-ков в масиві біля [5]}
procedure form (var x: aa; var y: aa);
var i: integer;
Begin
for i: = 1 to n do
y [i]: = sin (x [i]); {функція}
end;
...
