N I Ae A ? AE A I E A
no?.
Iauay iinoaiiaea caaea/e . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . 3
Iinoaiiaea oanoiauo caaea/ . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . 4
Iaoiaeeea ?aoaiey oanoiauo caaea/ . . . . . . . . . . . . . . . . . . .
. . . . 6
?acoeueoaou au/eneaiee . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . .9
Nienie eeoa?aoo?u . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . .10
I?eeiaeaiey
I?eeiaeaiea 1: Iienaiea i?ia?aiiu
I?eeiaeaiea 2: Oaeno i?ia?aiiu
1. I A U A ss I I N O A I I A E A C A Ae A * E
Ia?aiin oaiea (eee aauanoaa) oaieii?iaiaeiinoueth (aeey aauanoaa
niioaaonoaaiii aeeoooceae) e eiiaaeoeeae iienuaaaony aeeooa?aioeeaeueiui
o?aaiaieai ia?aaiee/aneiai oeia:
( 1 )
. O?aaiaiea (1) i?e oeacaiiuo auoa i?aaeiieiaeaieyo iacuaaaony
iiaeaeueiui o?aaiaieai aeenneiaoeee, eiiaaeoeee e eeiaoeee. Neaaaaiua
i?aaie /anoe eiatho neaaeothuee oece/aneee niune:
– niioaaonoaoao ia?aiino oaiea oaieii?iaiaeiinoueth (eee aauano-
aa aeeoooceae);
– niioaaonoaoao eiiaaeoeaiiio ia?aiino;
– “eeiaoe/aneee /eai”, niioaaonoaoao enoi/ieeo, i?iii?oeeiiaeueii-
io oaiia?aoo?a eee eiioeaio?aoeee;
– eioaineaiinoue aiaoieo enoi/ieeia eee noieia.
A aeaeueiaeoai aoaeai ?anniao?eaaoue oieueei oaieiaoth
eioa?i?aoaoeeth o?aaiaiey (1).
:
( 2 )
( 3 )
e a?aie/iuo oneiaeyo.
A?aie/iua oneiaey iienuaatho ?aaeeiu oaieiiaiaia n aiaoiae n?aaeie:
;
.
– 3 –
2. I I N O A I I A E A O A N O I A U O C A Ae A *
iiith auee aua?aiu neaaeothuea iyoue ooieoeee:
( 9 )
( 10 )
( 11 )
( 12 )
( 13 )
Aeey ooieoeee (10):
Aeey ooieoeee (11):
– 4 –
Aeaiiua ooieoeee oanoe?iaaeenue ia io?acea ii X: [0, 1], ii a?aiaie:
[0, 1], n eiee/anoaii ?acaeaiee ii yoei io?aceai – 30.
– 5 –
3. I A O I Ae E E A ? A O A I E ss O A N O I A U O
C A Ae A *
Aeaiiay caaea/a ?aoaaony n iiiiuueth aeaooneieiie iayaii
eiia/ii-?aciinoiie noaiu.
Noaia ?aaeecoaony a o?e yoaia.
n iiiiuueth 4-o oi/a/iie iayaiie noaiu:
( 5 )
. A yoii neo/aa eniieuecoaony 4-o oi/a/iay iayaiay ?aciinoiay noaia:
( 6 )
( 7 )
iaoiaeyony a aeaea eeiaeiie eiiaeiaoeee aeaoo i?aaeaa?eoaeueiuo
cia/aiee:
( 8 )
Aeey ?aoaiey (1) ainiieuecoainy oi?ioeaie (5) – (8). Aeaiiua o?aaiaiey
i?aaenoaaeytho o?ao aeeaaiiaeueiua iao?eoeu, ?aoaaiua iaoiaeii neaey?iie
i?iaiiee.
A ia/aea ioaeii i?aia?aciaaoue (5) – (7) e aeaeo:
( 14 )
Oiaaea (5) i?eiao aeae:
;
;
;
.
– 6 –
Oi?ioea (6) i?aia?acoaony a:
;
;
;
.
Oi?ioea (7) i?aia?acoaony a:
;
;
;
.
Aeaeaa ?aoaai ii oi?ioeai neaey?iie i?iaiiee:
( 15 )
( 16 )
ec oi?ioeu (9), oi eiaai:
.
– 7 –
( 18 )
. Aeaeaa iu iiaeai ?aoeoue nenoaiu iaoiaeii i?iaiiee e iieo/eoue
o?aaoaiue ?acoeueoao.
– 8 –
4. ? A C O E Ue O A O U A U * E N E A I E E
A ?acoeueoaoa i?iaaae?iiuo eniuoaiee i?ia?aiia iieacaea naith aunieoth
iaae?aeiinoue. Auee iieo/aiu neaaeothuea aeaiiua.
ia io?acea ii X e ii a?aiaie [0,1] n oaaii 0,033 aue iieo/ai
?acoeueoao n ioeaeie ?aaiie 0,0675.
, ia oii aea i?iiaaeooea, ioeaea ninoaaeyao 0,055.
ioeaea i?eiao cia/aiea 0,0435.
a ?acoeueoaoa aicieeaao ioeaea ?aaiay 0,0055.
, oi ioeaea ninoaaeo 0,00255.
O.a. iiaeii neacaoue, /oi iu eiaai ?acoeueoao n ia?aui ii?yaeeii
oi/iinoe. Noieue iaeoth oi/iinoue iiaeii iauynieoue oai, /oi
i?iecaiaeiay, iaeaeaiiay i?e a?aie/iuo oneiaeyo, oae aea eiaao ia?aue
ii?yaeie oi/iinoe.
– 9 –
N I E N I E E E O A ? A O O ? U
A. Aiaiaoieeia, A. Aiaiaoieeia I?ia?aiie?iaaiea a n?aaea Turbo-Pascal
7.0. – I.: Aeeaeia – Ieoe, 1996. – 288 n.
Iaoooiaa O. I., Neae?oeaa A. A. Iaeao i?eeeaaeiuo i?ia?aii aeey
/eneaiiiai iiaeaee?iaaiey i?ioeannia oaiei- e ianniia?aiina. –
Ea?aaaiaea: Ecae-ai Ea?AO. 1993
Oeao?iia A. Y. IBM PC aeey iieueciaaoaey. – I.: Eio?a – I, 1995. – 432
n.
I?eeiaeaiea 1
I I E N A I E A I ? I A ? A I I U
Iinoaaeaiiay caaea/a auea i?ia?aiiii ?aaeeciaaia ia ycuea
i?ia?aiie?iaaiey Turbo-Pascal 7.0.
A ninoaa i?ia?aiiu aoiaeyo neaaeothuea oaeeu:
basis.pas – PAS-oaee iniiaiie /anoe i?ia?aiiu
(?aoaiea nenoaiu o?aaiaiee iaoiaeii
neaey?iie i?iaiiee);
basis.v&v – EXE-oaee iniiaiie /anoe i?ia?aiiu (aucuaaaony ec START.PAS);
fun.bmp – BMP-oaee n ecia?aaeaieai ooieoeee;
inform.v&v – TXT-oaee n eioi?iaoeeae i i?ia?aiia (aucuaaaony ec
START.PAS);
music.v&v – iocueaeueiue EXE-oaee (aucuaaaony ec START.PAS);
my_menu.pas – UNIT aeey nicaeaiey iaith;
sea.exe – i?ia?aiia aeey i?iniio?a a?aoe/aneeo oaeeia;
start.pas – oaee aeey caionea anae i?ia?aiiu;
u – oaee n ?acoeueoaoaie ?aaiou;
zastavka.v&v – EXE-oaee n canoaaeie e iniiaiie i?ia?aiia
(aucuaaaony ec START.PAS).
Oaee START yaeyaony, eae au iaiei/eie i?ia?aiiu, ec eioi?ie
aucuaathony ae?oaea oaeeu. Nai i?ioeann ?aoaiey niaea?aeeony a oaeea
BASIS.
BASIS niaea?aeeo neaaeothuea i?ioeaaeo?u e ooieoeee:
Function Fun_U (Xm,t:real):real;
Aoiae: cia/aiea ii X e cia/aiea ii a?aiaie t, a oaeaea
aeiaaeueiay ia?aiaiiay aua?aiiie
ooieoeee SelectFunction.
Aeaenoaea: au/eneyao oi/iia cia/aiea ooieoeee U i?e caaeaiiuo X e t.
Auoiae: Fun_U – cia/aiea ooieoeee.
Function Fun_F (Xm,t,a,b,v:real):real;
e iiia? aua?aiiie ooieoeee
SelectFunction.
.
Auoiae: Fun_F – cia/aiea ooieoeee F.
Function Betta_Zero (time:real): real;
, iiia? aua?aiiie
ooieoeee SelectFunction.
, eniieuecoaiia a iaoiaea neaey?iie i?iaiiee.
.
Function U_End (time,Alf,Bet:real): real;
, iiia? aua?ai-
iie ooieoeee SelectFunction.
eniieuecoaiia a iaoiaea neaey?iie i?iaiiee.
.
Procedure PrintArray;
Aoiae: eniieuecoao aeiaaeueiue iannea aeaiiuo U_m.
Aeaenoaea: auaeaao niaea?aeeiia U_m ia ye?ai e a oaee.
Auoiae: auaiae U_m.
I?eeiaeaiea 2
O A E N O I ? I A ? A I I U
Iniiaiay /anoue i?ia?aiiu auaeyaeeo oae:
Program Basis;
Uses Crt; { Iiaeeeth/aiea aeaeeioae }
Label Metka1,Metka2; { Iaoee }
Var
a, b, v : real; { Eiyooeoeeaiou, caaeathony iieueciaaoaeai }
h, tau : real; { Oaa ii X e ii a?aiaie niioaaonoaaiii }
X,x0 : real; { Eiia/iia e ia/aeueiia cia/aiea X }
m,n,k : word; { Ia?aiaiiua eniieuecoaiua a oeeeeao aeey ?an/aoa }
T,t0 : real; { Eiia/iia e ia/aeueiia cia/aiea a?aiaie }
Kol_voX, Kol_voT : word; { Eiee/anoai ?acaeaiee ii X e ii a?aiaie }
U_m,U_,_U_1_2,_U_1 : array [0..200] of real; { Ianneau ?acoeueoaoia
}
z : array [0..200] of real; { Iannea oi/iuo ?aoaiee }
Xm : real; { I?iiaaeooi/iue X }
Alfa,Betta : array [0..200] of real; { Iannea eiyooeoeeaioia
eniieuecoaiuo i?e neaey?iie i?iaiiea }
a_progonka, b_progonka, c_progonka, d_progonka : real; {
Eiyooeoeeaiou aeey neaey?iie i?iaiiee }
Error : real; { Cia/aiea ioeaee }
time : real; { Ia?aiaiiay a?aiaie }
ch : char; { Eiae iaaeaoie eeaaeoe }
SelectFunction:word; { Iiia? aua?aiiie ooieoeee }
U : text; { Ia?aiaiiay aeey auaiaea ?acoeueoaoa a oaee }
Alfa_1,Alfa_2,Betta_1,Betta_2 : real; { Eiyooeoeeaiou a?aie/iuo
oneiaee }
Data : word; { Ia?aiaiiay ?aaeeia aaiaea ia/aeueiuo aeaiiuo }
Function Fun_U (Xm,t:real):real; { Ooieoeey U (oi/iia ?aoaiea) }
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; { Ooieoeey F }
begin
if SelectFunction=1 then
Fun_F:=2*t-v*6*Xm+a*3*SQR(Xm)-b*(SQR(Xm)*Xm+SQR(t));
if SelectFunction=2 then
Fun_F:=3*SQR(Xm)*SQR(t)+10*Xm+4*SQR(t)*t*Xm-v*2*SQR(t)*t+
a*(2*Xm*SQR(t)*t+10*t+SQR(SQR(t)))-b*(SQR(Xm)*SQR(t)*t+10*Xm*t+Xm*SQR(SQ
R(t)));
if SelectFunction=3 then
Fun_F:=SQR(Xm)*COS(Xm*t)+4*SQR(Xm)*SIN(t)-v*(2*COS(Xm*t)*t-
Xm*SIN(Xm*t)*SQR(t)-8*COS(t))+a*(SIN(Xm*t)+Xm*t*COS(Xm*t)-8*COS(t)*Xm)-
b*(Xm*SIN(Xm*t)-4*SQR(Xm)*COS(t));
if SelectFunction=4 then
Fun_F:=EXP(Xm)-v*(t*EXP(Xm))+a*(t*EXP(Xm))-b*(t*EXP(Xm));
if SelectFunction=5 then
Fun_F:=EXP(t)-v*(-SIN(Xm))+a*(COS(Xm))-b*(SIN(Xm)+EXP(t));
end;
Function Betta_Zero (time:real): real; { Ooieoeey Betta[0] aeey
i?iaiiee }
begin
If SelectFunction=1 then
Betta_Zero:=(h/(Betta_1*h-Alfa_1))*(Alfa_1*3*SQR(x0)+
Betta_1*(SQR(x0)*x0+SQR(time)));
If SelectFunction=2 then
Betta_Zero:=(h/(Betta_1*h-Alfa_1))*(Alfa_1*(2*x0*SQR(time)*time+
10*time+SQR(SQR(time)))+Betta_1*(SQR(x0)*SQR(time)*time+10*x0*time+SQR(S
QR(time))*x0));
If SelectFunction=3 then
Betta_Zero:=(h/(Betta_1*h-Alfa_1))*(Alfa_1*(SIN(x0*time)+
x0*time*COS(x0*time)-8*x0*COS(time))+Betta_1*(x0*SIN(x0*time)-4*SQR(x0)*
COS(time)));
If SelectFunction=4 then
Betta_Zero:=(h/(Betta_1*h-Alfa_1))*(Alfa_1*(time*EXP(x0))+
Betta_1*(time*EXP(x0)));
If SelectFunction=5 then
Betta_Zero:=(h/(Betta_1*h-Alfa_1))*(Alfa_1*(COS(x0))+
Betta_1*(SIN(x0)+EXP(time)));
end;
Function U_End (time,Alf,Bet:real): real; { Ooieoeey Um aeey i?iaiiee }
begin
If SelectFunction=1 then
U_End:=(Alfa_2*h*3*SQR(X)+Betta_2*h*(SQR(X)*X+SQR(time))
+ Bet*Alfa_2)/(Alfa_2-Alf*Alfa_2+h*Betta_2);
If SelectFunction=2 then
U_End:=(Alfa_2*h*(2*X*SQR(time)*time+10*time+SQR(SQR(time)))+
Betta_2*h*(SQR(X)*SQR(time)*time+10*X*time+SQR(SQR(time))*X)
+Bet*Alfa_2)/(Alfa_2-Alf*Alfa_2+h*Betta_2);
If SelectFunction=3 then
U_End:=(Alfa_2*h*(SIN(X*time)+X*time*COS(X*time)-8*X*COS(time))+
Betta_2*h*(X*SIN(X*time)-4*SQR(X)*COS(time))+Bet*Alfa_2)/(Alfa_2-Alf*Alf
a_2+h*Betta_2);
If SelectFunction=4 then
U_End:=(Alfa_2*h*(time*EXP(X))+Betta_2*h*(time*EXP(X))+Bet*Alfa_2)/
(Alfa_2-Alf*Alfa_2+h*Betta_2);
If SelectFunction=5 then
U_End:=(Alfa_2*h*(COS(X))+Betta_2*h*(SIN(X)+EXP(time))+Bet*Alfa_2)/
(Alfa_2-Alf*Alfa_2+h*Betta_2);
end;
Procedure PrintArray; { I?ioeaaeo?a ia/aoe ianneaa U }
begin
WriteLn; For m:=0 to Kol_voX do begin Write(U_m[m]:15:4);
Write(U,U_m[m]:15:4); end;
WriteLn; WriteLn(U);
end;
{ Iniiaiay i?ia?aiia }
Begin
Assign(U,’u’); { Oaee aeey caiene cia/aiee ooieoeee }
Rewrite(U); { Ioe?uoea oaeea aeey caiene }
TextBackGround(0); { Auai? ooieoeee aeey ?aaiou }
ClrScr; TextColor(10); GoToXY(20,8); Write(‘Aaaaeeoa iiia?
aua?aiiie ooieoeee (1-5):’);
Metka1: ch:=ReadKey;
If ch=’1′ then SelectFunction:=1
else If ch=’2′ then SelectFunction:=2
else If ch=’3′ then SelectFunction:=3
else If ch=’4′ then SelectFunction:=4
else If ch=’5′ then SelectFunction:=5
else
begin
Sound(400); Delay(100);
NoSound; GoTo Metka1;
end;
GoToXY(59,8);TextColor(12);WriteLn(SelectFunction); TextColor(11);
GoToXY(11,12);
Write(‘Au aoaeaoa ?aaioaoue ni noaiaea?oiuie ia?aiao?aie (oeeo?a
~1~)’);
GoToXY(22,13); Write(‘eee aaaaeaoa naie aeaiiua (oeeo?a ~2~) ?’);
Metka2: ch:=ReadKey;
If ch=’1′ then Data:=1
else If ch=’2′ then Data:=2
else
begin
Sound(400); Delay(100); NoSound; GoTo Metka2;
end;
TextBackGround(9); TextColor(10); ClrScr;
{ Aaiae ia/aeueiuo aeaiiuo }
WriteLn; WriteLn(‘——————————– Aaiae aeaiiuo
———————————¬’);
For k:=1 do 21 do WriteLn(‘¦
¦’);
WriteLn(‘L————————————————————–
—————-‘);
TextColor(15); Window(3,3,77,23); Write(‘ Aaaaeeoa iaeanoue ?ann/aoa
ii X io: ‘);
If Data=1 then
begin
x0:=0; Write(x0:1:0); WriteLn;
end
else ReadLn(x0);
Write(‘ aei: ‘);
If Data=1 then
begin
X:=1; Write(X:1:0); WriteLn;
end
else ReadLn(X);
WriteLn; Write(‘ Aaaaeeoa eiee/anoai ?acaeaiee ii iai?aaeaieth X: ‘);
If Data=1 then begin Kol_voX:=30; Write(Kol_voX:2); WriteLn; end
else ReadLn(Kol_voX);
WriteLn;WriteLn; Write(‘ Aaaaeeoa iaeanoue ?ann/aoa ii a?aiaie io: ‘);
If Data=1 then begin t0:=0; Write(t0:1:0); WriteLn; end else
ReadLn(t0);
Write(‘ aei: ‘);
If Data=1 then begin T:=1; Write(T:1:0); WriteLn; end else
ReadLn(T);
WriteLn; Write(‘ Aaaaeeoa eiee/anoai ?acaeaiee ii a?aiaie: ‘);
If Data=1 then begin Kol_voT:=30; Write(Kol_voT:2); WriteLn; end
else ReadLn(Kol_voT);
WriteLn;WriteLn; WriteLn(‘ Aaaaeeoa eiyooeoeeaiou’); Write(‘ a=’);
If Data=1 then begin a:=1; Write(a:1:0); WriteLn; end else
ReadLn(a);
Write(‘ b=’);
If Data=1 then begin b:=1; Write(b:1:0); WriteLn; end else
ReadLn(b);
Write(‘ v=’);
If Data=1 then begin v:=0.001; Write(v:1:3); WriteLn; end else
ReadLn(v);
Write(‘ Alfa-1=’);
If Data=1 then begin Alfa_1:=1; Write(Alfa_1:1:0); WriteLn; end
else ReadLn(Alfa_1);
Write(‘ Betta-1=’);
If Data=1 then begin Betta_1:=1; Write(Betta_1:1:0); WriteLn; end
else ReadLn(Betta_1);
Write(‘ Alfa-2=’);
If Data=1 then begin Alfa_2:=1; Write(Alfa_2:1:0); WriteLn; end
else ReadLn(Alfa_2);
Write(‘ Betta-2=’);
If Data=1 then begin Betta_2:=1; Write(Betta_2:1:0);
WriteLn;TextColor(14);
Write(‘ Iaaeieoa ethaoth eeaaeoo’); ReadKey; end
else ReadLn(Betta_2);
{ Eioa?oaen ye?aia i?e auaea/a ?acoeueoaoa }
TextBackGround(3); TextColor(1); Window(1,1,80,25); ClrScr; WriteLn;
WriteLn(‘a===================== ?acoeueoao
==========================¬’);
For k:=1 to 21 do
WriteLn(‘¦
¦’);
WriteLn(‘==============================================================
=====-‘);
TextColor(0); TextBackGround(7); GoToXY(2,23);
WriteLn(‘ Aeey i?iaeieaeaiey iaaeieoa ethaoth eeaaeoo’);
TextBackGround(3); Window(3,4,77,22);
TextColor(15); ClrScr;
{ Au/aneaiea oaaa naoee }
tau:=(T-t0)/Kol_voT; h:=(X-x0)/Kol_voX;
{ Aaiae aeaiiuo i?e time=t0 }
For m:=0 to Kol_voX do
begin
Xm:=x0+h*m; U_m[m]:=Fun_U(Xm,t0);
end;
TextColor(14); WriteLn(‘A?aiy ?aaii ‘,time:3:3); TextColor(15);
WriteLn(U,’A?aiy ?aaii ‘,time:3:3);
PrintArray;
{ Ia/aei ?ann/aoa }
time:=t0;
Repeat
time:=time+tau;
WriteLn; TextColor(14); WriteLn(‘A?aiy ?aaii ‘,time:3:3);
TextColor(15);
WriteLn(U,’A?aiy ?aaii ‘,time:3:3);
{ 1 yoai. ?aoaaony iaoiaeii neaey?iie i?iaiiee }
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 yoai, /anoue ia?aay. ?aoaaony iaoiaeii neaey?iie i?iaiiee }
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])/
(a_progonka*Alfa[m-1]+b_progonka);
end;
_U_1_2[Kol_voX]:=U_End(time,Alfa[Kol_voX-1],Betta[Kol_voX-1]);
For m:=Kol_voX-1 downto 1 do _U_1_2[m]:=Alfa[m]*_U_1_2[m+1]+Betta[m];
_U_1_2[0]:=Alfa[0]*_U_1_2[1]+Betta[0];
{ 2 yoai, /anoue aoi?ay. ?aoaaony iaoiaeii neaey?iie i?iaiiee }
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,a,b,v)+2*_U_1_2[m]/tau-a_progonka*Betta
[m-1])/
(a_progonka*Alfa[m-1]+b_progonka);
end;
_U_1[Kol_voX]:=U_End(time,Alfa[Kol_voX-1],Betta[Kol_voX-1]);
For m:=Kol_voX-1 downto 1 do _U_1[m]:=Alfa[m]*_U_1[m+1]+Betta[m];
_U_1[0]:=Alfa[0]*_U_1[1]+Betta[0];
{ 3 yoai. Ieii/aoaeueiia cia/aiea }
For m:=0 to Kol_voX do
U_m[m]:=2*_U_1[m]-U_[m];
PrintArray; { Auaiae ?acoeueoaoa ia ye?ai e aai caienue a oaee }
For m:=0 to Kol_voX do { ?ann/ao oi/iiai cia/aiey ooieoeee }
begin z[m]:=Fun_U(x0+m*h,time); end;
{ Auaiae ioeaee ?an/aoa ia ye?ai e a oaee }
Error:=0;
For m:=0 to Kol_voX do
begin
a:=Abs(U_m[m]-z[m]); If ErrorT;
Close(U); { Cae?uoea oaeea n ?acoeueoaoaie }
End.
Ieienoa?noai iauaai e i?ioanneiiaeueiiai
ia?aciaaiey ?O
I?aiao?aneee ainoaea?noaaiiue oieaa?neoao
Einoeooo yia?aaoeee e eioi?iaoeee
eaoaae?a: Eioi?iaoeee
? A N * A O I I – A ? A O E * A N E I A
C A Ae A I E A
« *eneaiiia ?aoaiea iiaeaeueiiai o?aaiaiey
aeenneiaoeee, eiiaaeoeee e eeiaoeee »
Auiieiee : nooaeaio a?. 97 EAeIA
Aieeia A. A
.
Distributed by BRS Corporation
http://www.osu.ru/~BRS
E-mail: brs-99@mail.ru
I?iaa?ee : Iaoooiaa O. I.
I?aiao?a – 1999
( 4 )
( 17 )
Нашли опечатку? Выделите и нажмите CTRL+Enter
