С Подпрограмма-функция

FUNCTION FUN(X)

В качестве подпрограммы-функции может быть любая функция с любым именем.

FUN=SIN(X)**2–COS(X)**2

RETURN

END

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

EXTERNAL FUN

WRITE (6,*)’Введите Х_нач , Х_кон’

READ (5,*) A,B

Обращаемся к стандартной подпрограмме вычисления определенного интеграла. Фактические параметры А и В, имя функции FUN. Результат вычисления интеграла будет помещен в переменную с именем Y.

CALL QG5(A,B,FUN,Y)

WRITE (6,*)’Вычисленное значение интеграла =’,Y

PAUSE

STOP

END

Поскольку имя подпрограммы-функции FUN находится в списке фактических параметров при обращении к стандартной программе QG5, то имя этой функции должно быть описано в операторе EXTERNAL в самом начале программы. Для компилятора это является признаком, что FUN не простая переменная, а функция с именем FUN.

Задача:

Решить дифференциальное уравнение и построить график для:

на отрезке с шагом dx при .

Использовать стандартную программу RK2(FUN,H,X1,Y1,K,N,VEC)

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

Н – шаг интегрирования;

X1, Y1 – начальные значения;

K – количество шагов между двумя сохраняемыми значениями;

N– количество значений;

VEC – результирующий вектор-столбец длинны N, в котом сохраняются вычисленные значения.

Величина K задаёт точность решения (количество промежуточных вычисленных точек, которые не запоминаются).