Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
15
Добавлен:
17.04.2013
Размер:
8 Кб
Скачать



 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.








Соседние файлы в папке curstau