Курсовые / curstau / SNT2

Program SNT2;
uses Crt,GB,ENEX,ALG,SAMP,CANON,FINIT,CLSY,RC;
label 1;
var n,k,i,j,w :integer;
l,tk,nr,dc,dh,Um :real;
A,B,C,Ad,Bd,Q,Sco,Sob,
Rx,Hx,Acr,Ach,Ao,P,E,Ak,
AAcr,AAAcr,AAch,AAAch :matr;
X,Xo,_X,U,Xoo,_Xo :vect;


BEGIN
1:TextBackground(9);
TextColor(15);
ClrScr;
{ writeln(' €ў ­®ў ….Ђ. 31.01.96 ');}
writeln;
write(' ” ¬Ё«Ёп бв㤥­в  - ');
readln;
write(' ѓагЇЇ  - ');
readln;
writeln;
writeln(' ‘ Ё ­ в Ґ § ');
writeln(' жЁда®ў®© бЁб⥬л гЇа ў«Ґ­Ёп ');
writeln;
writeln(' ‚ў®¤ ¬ ваЁж ­ҐЇаҐалў­®© ¬®¤Ґ«Ё ®ЎкҐЄв ');
write (' ђ §¬Ґа­®бвм ¬ ваЁжл A: n=');
readln(n);
writeln(' Њ ваЁж  A');
EN(11,n,n,A);
writeln(' Њ ваЁж  B');
EN(12+n,n,1,B);
writeln(' Њ ваЁж  C');
EN(13+2*n,1,n,C);
ClrScr;
writeln;
writeln;
writeln(' „ЁбЄаҐвЁ§ жЁп ­ҐЇаҐалў­ле га ў­Ґ­Ё© б®бв®п­Ёп');
writeln;
write('ЏҐаЁ®¤ ¤ЁбЄаҐв­®бвЁ tk=');
readln(tk);
SMP(tk,n,1,A,B,Ad,Bd);
writeln;
writeln(' „ЁбЄаҐв­ п ¬®¤Ґ«м');
writeln(' ¬ ваЁж  Ad');
writeln;
EX(n,n,Ad);
writeln(' ¬ ваЁж  Bd');
writeln;
EX(n,1,Bd);
ReadLn;
ClrScr;
writeln;
writeln(' Њ ваЁж  гЇа ў«пҐ¬®бвЁ ¤ЁбЄаҐв­®© ¬®¤Ґ«Ё');
CONT(n,Ad,Bd,Sco);
writeln;
EX(n,n,Sco);
DET(n,Sco,dc);
writeln(' det Sco =',dc);
ReadLn;
if abs(dc)<10E-12 then
begin
writeln(' det Sco < 10E-12 - ®ЎкҐЄв ­ҐгЇа ў«пҐ¬ - гв®з­ЁвҐ ¬®¤Ґ«м ®ЎкҐЄв ');
ReadLn;
goto 1;
end else
writeln(' Њ ваЁж  ­ Ў«о¤ Ґ¬®бвЁ ¤ЁбЄаҐв­®© ¬®¤Ґ«Ё');
OBS(n,Ad,C,Sob);
writeln;
EX(n,n,Sob);
DET(n,Sob,dh);
writeln(' det Sob =',dh);
ReadLn;
if abs(dh)<10E-12 then
begin
writeln(' det Sob < 10E-12 - ®ЎкҐЄв ­Ґ­ Ў«о¤ Ґ¬ - гв®з­ЁвҐ ¬®¤Ґ«м ®ЎкҐЄв ');
ReadLn;
goto 1;
end else
ClrScr;
writeln;
writeln(' ђ бзҐв ॣг«пв®а  б®бв®п­Ёп б Є®­Ґз­л¬ ўаҐ¬Ґ­Ґ¬ ');
writeln(' ЇҐаҐе®¤­®Ј® Їа®жҐбб  ');
writeln;
writeln(' Њ ваЁж  ॣг«пв®а  б®бв®п­Ёп Rx ');
REGT(n,Ad,Bd,Rx);
writeln;
EX(1,n,Rx);
readln;
writeln(' Њ ваЁж  ®Ўа в­®© бўп§Ё ­ Ў«о¤ вҐ«п б®бв®п­Ёп Hx ');
writeln;
OBST (n,Ad,C,Hx);
EX(n,1,Hx);
readln;

ClrScr;
writeln(' ‡ ¬Є­гв п бЁб⥬  ');
writeln;
writeln(' Њ ваЁж  § ¬Є­гв®© бЁб⥬л б ॣг«пв®а®¬ б®бв®п­Ёп');
writeln(' Acr=Ad-Bd*Rx ');
CLSYR(n,Ad,Bd,Rx,Acr);
writeln;
EX(n,n,Acr);
readln;
MUL(n,n,n,Acr,Acr,AAcr);
MUL(n,n,n,AAcr,Acr,AAAcr);
writeln;
writeln(' Acr*Acr*Acr ');
writeln;
EX(n,n,AAAcr);
readln;
ClrScr;
writeln;
writeln(' ‘®Ўб⢥­­ п ¬ ваЁж  ­ Ў«о¤ вҐ«п ');
writeln(' Ach=Ad-Hx*C ');
CLSYH(n,Ad,C,Hx,Ach);
writeln;
EX(n,n,Ach);
readln;
MUL(n,n,n,Ach,Ach,AAch);
MUL(n,n,n,AAch,Ach,AAAch);
writeln;
writeln(' Ach*Ach*Ach ');
writeln;
EX(n,n,AAAch);
readln;
ClrScr;
writeln;
writeln(' ‚ҐЄв®а б®бв®п­Ёп § ¬Є­гв®© бЁб⥬л ');
writeln(' c ॣг«пв®а®¬ Ё ­ Ў«о¤ вҐ«Ґ¬ ');
writeln;
writeln(' Э x1(k)Э і ');
writeln(' Э x2(k)Э і - ЇҐаҐ¬Ґ­­лҐ б®бв®п­Ёп ®ЎкҐЄв  ');
writeln(' X(k)=Э x3(k)Э х _ ');
writeln(' Эdx1(k)Э і - dxi(k)=xi(k)-xi(k), ');
writeln(' Эdx2(k)Э і _ ');
writeln(' Эdx3(k)Э і xi(k) - ЇҐаҐ¬.б®бв.­ Ў«о¤ вҐ«п');
readln;
ClrScr;
writeln;
writeln(' Њ ваЁж  § ¬Є­гв®© бЁб⥬л б ॣг«пв®а®¬ б®бв®п­Ёп');
writeln(' Ё ­ Ў«о¤ вҐ«Ґ¬ ');
writeln(' Ao=Ю Acr Bd*Rx Ю ');
writeln(' Ю 0 Ach Ю ');
CLS(n,Ad,Bd,C,Rx,Hx,Ao);
writeln;
EX(2*n,2*n,Ao);
readln;
ClrScr;
writeln;
writeln(' ‚ҐЄв®а ­ з «м­ле гб«®ўЁ© ');
writeln(' xi(0) - ­ з «м­лҐ гб«®ўЁп ¤«п ЇҐаҐ¬.б®бв.®ЎкҐЄв  ');
writeln;
write(' x1(0)=');
readln(Xo[1]);
write(' x2(0)=');
readln(Xo[2]);
write(' x3(0)=');
readln(Xo[3]);
writeln(' _ ');
writeln(' xi(0) - ­ з «м­лҐ гб«®ўЁп ¤«п ЇҐаҐ¬.б®бв.­ Ў«о¤ вҐ«п ');
writeln;
writeln(' _ ');
write(' x1(0)=');
readln(_X[1]);
writeln(' _ ');
write(' x2(0)=');
readln(_X[2]);
writeln(' _ ');
write(' x3(0)=');
readln(_X[3]);
readln;
for i:=1 to n do _Xo[i]:=_X[i];
for i:=1 to n do Xo[n+i]:=Xo[i]-_X[i];
for i:=1 to 2*n do Xoo[i]:=Xo[i];
ClrScr;
writeln(' ђҐиҐ­ЁҐ га ў­Ґ­Ё© б®бв®п­Ёп ');
writeln;
writeln(' _ _ _ ');
writeln(' k і x1(k) і x1(k) і x2(k) і x2(k) і x3(k) і x3(k) і u(k) і');
writeln('---------------------------------------------------------------------');
writeln(' ',0 , ' і',Xo[1]:4,'і',_X[1]:4,'і',Xo[2]:4,'і',_X[2]:4,'і',Xo[3]:4,'і',_X[3]:4,'і',U[1]:4);
writeln;
k:=0;
repeat
MULV(2*n,2*n,Ao,Xo,X);
for i:=1 to n do _X[i]:=X[i]-X[n+i];
MULV(1,n,Rx,_X,U);
U[1]:=-U[1];
if k=0 then Um:=abs(U[1]);
if abs(U[1])>Um then Um:=abs(U[1]);
for i:=1 to 2*n do Xo[i]:=X[i];
k:=k+1;
writeln(' ',k:2, ' і',X[1]:8,'і',_X[1]:8,'і',X[2]:8,'і',_X[2]:8,'і',X[3]:8,'і',_X[3]:8,'і',U[1]:8);
readln;
until k=10;
writeln;
writeln(' іUmaxі =',Um);
readln;


repeat
ClrScr;
writeln;
write(' Ќ ¤® Ё§¬Ґ­пвм ЇҐаЁ®¤ ¤ЁбЄаҐв­®бвЁ tk? 1/0 - ');;
readln(W);
if W=1 then
begin
writeln;
write(' ЏҐаЁ®¤ ¤ЁбЄаҐв­®бвЁ tk=');
readln(tk);
SMP(tk,n,1,A,B,Ad,Bd);
REGT(n,Ad,Bd,Rx);
OBST (n,Ad,C,Hx);
CLSYR(n,Ad,Bd,Rx,Acr);
CLSYH(n,Ad,C,Hx,Ach);
CLS(n,Ad,Bd,C,Rx,Hx,Ao);
for i:=1 to 2*n do Xo[i]:=Xoo[i];
writeln(' _ _ _ ');
writeln(' k і x1(k) і x1(k) і x2(k) і x2(k) і x3(k) і x3(k) і u(k) і');
writeln('---------------------------------------------------------------------');
writeln(' ',0 , ' і',Xo[1]:4,'і',_Xo[1]:4,'і',Xo[2]:4,'і',_Xo[2]:4,'і',Xo[3]:4,'і',_Xo[3]:4,'і',U[1]:4);
writeln;
k:=0;
repeat
MULV(2*n,2*n,Ao,Xoo,X);
for i:=1 to n do _X[i]:=X[i]-X[n+i];
MULV(1,n,Rx,_X,U);
U[1]:=-U[1];
if k=0 then Um:=abs(U[1]);
if abs(U[1])>Um then Um:=abs(U[1]);
for i:=1 to 2*n do Xoo[i]:=X[i];
k:=k+1;
writeln(' ',k:2, ' і',X[1]:8,'і',_X[1]:8,'і',X[2]:8,'і',_X[2]:8,'і',X[3]:8,'і',_X[3]:8,'і',U[1]:8);
readln;
until k=10;
writeln;
writeln(' іUmaxі =',Um);
readln;
end;
until W=0;
END.