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

end;

TextColor(14); WriteLn('Времяравно ',time:3:3); TextColor(15); WriteLn(U,'Времяравно ',time:3:3);

PrintArray;

{ Начало рассчета }

time:=t0;

Repeat

time:=time+tau;

WriteLn; TextColor(14); WriteLn('Времяравно ',time:3:3); TextColor(15);

WriteLn(U,'Времяравно ',time:3:3);

{ 1 этап. Решается методом скалярной прогонки }

a_progonka:=(-2*v-a*h)/(2*SQR(h)); b_progonka:=(SQR(h)+2*v*tau-b*tau*SQR(h))/(SQR(h)*tau);

c_progonka:=(a*h-2*v)/(2*SQR(h));

Alfa[0]:=Alfa_1/(Alfa_1-Betta_1*h); Betta[0]:=Betta_Zero(time);

For m:=1 to Kol_voX-1 do

begin

Alfa[m]:=-c_progonka/(a_progonka*Alfa[m-1]+b_progonka);

Betta[m]:=(Fun_F(x0+m*h,time+tau,a,b,v)+U_m[m]/tau-a_progonka*Betta[m-1])/

(a_progonka*Alfa[m-1]+b_progonka);

end;

U_[Kol_voX]:=U_End(time,Alfa[Kol_voX-1],Betta[Kol_voX-1]);

For m:=Kol_voX-1 downto 1 do U_[m]:=Alfa[m]*U_[m+1]+Betta[m];U_[0]:=Alfa[0]*U_[1]+Betta[0];

{ 2 этап, часть первая. Решается методом скалярной прогонки }

a_progonka:=(-2*v-a*h)/(2*SQR(h)); b_progonka:=(2*SQR(h)+2*v*tau-b*tau*SQR(h))/(SQR(h)*tau);

c_progonka:=(a*h-2*v)/(2*SQR(h));

Alfa[0]:=Alfa_1/(Alfa_1-Betta_1*h); Betta[0]:=Betta_Zero(time);

For m:=1 to Kol_voX-1 do

begin

Alfa[m]:=-c_progonka/(a_progonka*Alfa[m-1]+b_progonka);

Betta[m]:=(Fun_F(x0+m*h,time+tau/2,a,b,v)+2*U_m[m]/tau-a_progonka*Betta[m-1])/