Вычисление интеграла фукции f (x) (методом Симпсона WinWord) — страница 4

  • Просмотров 1102
  • Скачиваний 26
  • Размер файла 128
    Кб

соответствии с приведённым словесным описанием алгоритма решения поставленной задачи разработана блок схема решаемой задачи, которая изображена на рис. 3. В изображенном алгоритме блоки имеют описанное ниже назначение: Блок 1. Начало программы; Блок 2. Очистка экрана;; Блок 3. Запрос на ввод значений А и В; Блок 4. Ввод значений А и В с клавиатуры; Блок 5. Вызов процедуры вывода графика функции на экран; Блок 6. Установка начального

значения счётчика отрезков равным 3; Блок 7. Вычисление значения начального значения интеграла методом трапеций; Блок 8. Запоминание предыдущего значения интеграла, вычисленного методом трапеций, увеличение значения числа отрезков на 2, вычисление следующего значения интеграла методом трапеций; Блок 9. Проверка условия : абсолютное значение разности текущего и предыдущего значений интегрирования меньше чем 0.001, если да, то

выход из цикла, если нет, то переход на блок 8. Блок 10. Вывод результатов, полученных при вычислении интеграла методом трапеций на экран. Блок 11. Установка начального значения счётчика отрезков равным 3; Блок 12. Вычисление значения начального значения интеграла методом Симпсона; Блок 13. Запоминание предыдущего значения интеграла, вычисленного методом Симпсона, увеличение значения числа отрезков на 2, вычисление следующего

значения интеграла методом Симпсона; Блок 14. Проверка условия: абсолютное значение разности текущего и предыдущего значений интегрирования меньше чем 0.001, если да, то выход из цикла, если нет, то переход на блок 13. Блок 15. Вывод результатов, полученных при вычислении интеграла методом Симпсона на экран. Блок 16. Конец программы. 5. Текст программы program tr_s; uses crt,graph; var   a,b:real;  { Границы отрезка }   r,r2:real;  { Предыдущее и текущее

приближенные значения интеграла}   n:integer;   { Счетчик } { Интегрируемая функция } function f(x:real):real; begin   f:=1/(x*ln(x)*0.43429); end; { Метод трапеций } function trap(a,b:real;n:integer):real; var   s:real;      { Полученная сумма }   h:real;      { Шаг }   m:integer;   { Счетчик } begin   h:=(b-a)/(n-1);                   { Определяется шаг }   s:=(f(a)+f(b))/2;                 {

Начальное значение суммы }   for m:=1 to n-2 do s:=s+f(a+m*h); { Суммиование остальных элементов}   trap:=s*h;                        { Возвращается значение интеграла } end; { Метод Симпсона } function simpson(a,b:real;n:integer):real; var   s:real;      { Сумма }   h:real;      { Шаг }   m:integer;   { Счетчик }   mn:integer;  { Очередной множитель } begin