Замечание. 1. В приведенной программе реализованы все три рассматриваемые в настоящем курсе формулы численного интегрирования: п/п PR – вычисление интеграла по

1. В приведенной программе реализованы все три рассматриваемые в настоящем курсе формулы численного интегрирования:

• п/п PR – вычисление интеграла по формуле прямоугольников;
• п/п TR – вычисление интеграла по формуле трапеций;
• п/п SIMPS – вычисление интеграла по формуле Симпсона.

2. Переменная (в программе – kmax) определяет количество отрезков, на которые делится интервал интегрирования для вычисления интеграла с заданной точностью. При этом значение n последовательно удваивается до тех пор, пока изменение приближенного значения интеграла не станет меньше заданного малого числа (в программе – eps), значение которого вводится. При расчете было принято:

kmax =20, eps =0.0001.

3. К данной программе (проекту) необходимо подключить стандартную подпрограмму QATR.FOR.

4. Формальному параметру sint (имени подпрограммы) в подпрограмме-функции aintegral соответствуют фактические параметры pr,tr и simps (имена подпрограмм-функций). Эти параметры, а также фактический параметр f в обращении к стандартной подпрограмме qatr, должны быть перечислены в списке оператора external:

external f,pr,tr,simps.

5. Параметры a0,a1,a2 вычисляются в основной программе и передаются в подпрограмму-функцию f оператором общих областей common:

common a0,a1,a2.