Приложение б

(справочное)

Образец оформления титульного листа

Министерство образования и науки Украины

Севастопольский национальный технический университет

Кафедра технической кибернетики

Курсовая работа

по дисциплине «Вычислительные методы»

Студент (ф.и.о.)

Группа _______________________

№ уравнения_____№ метода_____

Руководитель (должность, ф.и.о.)

Дата защиты________________Оценка___________

Подписи членов комиссии

________________

________________

Севастополь- 2003г.

Приложение в

(справочное)

Процедура интегрирования

Procedure RK4 (var t: real; h: real; n: integer; var x: vector_n; right: type_right);

Назначение Выполнение одного шага интегрирования системы обыкновенных дифференциальных уравнений x'=f(x(t),t) методом Рунге-Кутты четвёртого порядка. Основные локальные переменные h2,h6- коэффициенты Рунге-Кутты; f1,f2,f3,f4 – массивы значений Fi (3.2) правой части уравнения системы (2.4) для i=1,2,3,4 соответственно; xr- массив значений векторных аргументов правой части уравнения системы при вычислении Fi , i=1,2,3,4 Формальные параметры t- независимая переменная. При входе в процедуру- значение в текущем узле интегрирования, при выходе - в следующем. h- шаг интегрирования; n- порядок системы; x- массив решений системы. При входе в процедуру- решение в текущем узле интегрирования, при выходе - в следующем; right- имя процедуры вычисления правой части уравнения системы, где type_right- имя процедурного типа.

var x,xr,f1,f2,f3,f4: vector n;

i: integer;

h2,h6: real;

begin

h2:= 0.5*h;

h6:=0.16666666666666*h;

{Первый этап}

right (t,x,f1);

for i:=1 to n do xr[i]:=x[i]+h2*f1[i];

{Второй этап}

right (t+h2,xr,f2);

for i:=1 to n do xr[i]:=x[i]+h2*f2[i];

{Третий этап}

right (t+h2,xr,f3);

for i:=1 to n do xr[i]:=x[i]+h*f3[i];

{Четвёртый этап}

t:= t+h;

right (t,xr,f4);

for i:=1 to n do x[i]:=x[i]+h6*(f1[i]+2*(f2[i[+f3[i])+f4[i])

end.