Реферат: Численное решение модельного уравнения

a, b, v : real; { Коэффициенты, задаются пользователем }

h, tau : real; { Шаг по X и по времени соответственно }

X,x0 : real; { Конечное и начальное значение X }

m,n,k : word; { Переменные используемые в циклах для расчета }

T,t0 : real; { Конечное и начальное значение времени }

Kol_voX, Kol_voT : word; { Количество разбиений по X и по времени }

U_m,U_,_U_1_2,_U_1 : array [0..200] of real; { Массивы результатов }

z : array [0..200] of real; { Массив точных решений }

Xm : real; { Промежуточный X }

Alfa,Betta : array [0..200] of real; { Массив коэффициентов используемых при скалярной прогонке }

a_progonka, b_progonka, c_progonka, d_progonka : real; { Коэффициенты для скалярной прогонки }

Error : real; { Значение ошибки }

time : real; { Переменная времени }

ch : char; { Код нажатой клавиши }

SelectFunction:word; { Номер выбранной функции }

U : text; { Переменная для вывода результата в файл }

Alfa_1,Alfa_2,Betta_1,Betta_2 : real; { Коэффициенты граничных условий }

Data : word; { Переменная режима ввода начальных данных }

Function Fun_U (Xm,t:real):real; { Функция U (точное решение) }

begin

If SelectFunction=1 then Fun_U:=SQR(Xm)*Xm+SQR(t);

If SelectFunction=2 then Fun_U:=SQR(Xm)*SQR(t)*t+10*Xm*t+SQR(SQR(t))*Xm;

If SelectFunction=3 then Fun_U:=Xm*SIN(Xm*t)-4*SQR(Xm)*COS(t);

If SelectFunction=4 then Fun_U:=t*EXP(Xm);

If SelectFunction=5 then Fun_U:=SIN(Xm)+EXP(t);

end;

Function Fun_F (Xm,t,a,b,v:real):real; { Функция F }

begin

if SelectFunction=1 then Fun_F:=2*t-v*6*Xm+a*3*SQR(Xm)-b*(SQR(Xm)*Xm+SQR(t));