Численные методы решения обыкновенных дифференциальных уравнений и систем: Методические указания к курсовой работе, страница 7

Задаваясь начальными значениями x₀=a , y₀  , последовательно вычислим для всех i=0,1, … , n-1

           (2.8)

                                                 (2.9)

x _i+1=x_i+h

где функция  f (x,y)  вычисляется по формуле (в данном примере):

                                                                      (2.10)

Алгоритм решения задачи

Исходными  данными являются: a, b, h, y0 В случае необходимости следует также  ввести параметры. В данном случае - это коэффициент c . Блок-схема решения задачи представлена на рис.2.4.

В разделе описаний следует описать необходимые переменные и функции, в соответствии с расчетными формулами. Их исходные обозначения, имена в программе, типы и назначение перечислим в табл. 2.1.

Таблица 2.1

Структура данных

Исходные обозначения величин	Обозначения величин в программе	Тип	Пояснения
a=x0, b	a, b	real	Начальное и конечное значение независимой переменной x (заданы)
y0	y0
h	h	real	Щаг изменения переменной x (задан)
xгр=b-h/2	xgr	real	Граничное значение переменной x
xi	x	real	Значение x в узле
yi	y	real	Значение искомой функции y и узловой точке
f(x, y)	f	Функция типа real	Правая часть уравнения
k1, k2, k3, k4	k1, k2, k3, k4	real	Приращения функции
c	c	real	Параметр уравнения