С Головная программа

DIMENSION A(3,3),C(3)

WRITE (6,*)’Введите коэффициенты при Х по строкам’

READ (5,*) ((A(I,J),J=1,3),I=1,3)

Введём: 2. 3. 4. -1. -2. 3. 1. -3. 0.

WRITE (6,*)’Введите массив свободных членов’

READ (5,*) C

Введём: 5. 4. 6.

CALL SIMQ(A,C,3,IER)

Идет обращение к стандартной программе SIMQ. В этой СП будет получено решение, а результат этого решения будет находиться в массиве С.

С Анализ кода ошибки

IF(IER.EQ.0) WRITE (6,*)’Решение: X1, X2, X3=’, (C(I),I=1,3)

IF(IER.EQ.1) WRITE (6,*)’Нет решения. Проверьте исходные данные’

IF(IER.EQ.2) WRITE (6,*)’Неверное количество уравнений’

PAUSE

STOP

END

Задача:

Вычислить значение определенного интеграла:

Воспользоваться программой из БСП: QG5(XL,XU,FCT,Y)

где XL – левый предел интеграла;

XU – правый предел интеграла;

FCT – внешняя подпрограмма-функция, которая описывает решаемое уравнение ;

Y – результат.

Эта программа с автоматическим выбором шага и для своей работы требует использование подпрограммы-функции, в которой описывается решаемое уравнение . Начинаем с написания подпрограммы – функции.