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

 Program SNT;
       uses Crt,GB,ENEX,ALG,SAMP,CANON,FINIT,CLSY,RC;
          label 1;
          var n,k,i,j                    :integer;
              l,tk,nr,dc,dh              :real;
              A,B,C,Ad,Bd,Q,Sco,Sob,
              Rsd,Hx,Acr,Ach,Ao,P,E,Ak   :matr;
              X,Xo,_X,U                  :vect;
    BEGIN
    1:TextBackground(9);
      TextColor(15);
      ClrScr;
      writeln('  €ў ­®ў ….Ђ.  12.05.95            ');
      writeln;
      writeln('                    ‘ Ё ­ в Ґ §          ');
      writeln('            жЁда®ў®© бЁб⥬л гЇа ў«Ґ­Ёп  ');
      writeln;
      writeln('    ‚ў®¤ ¬ ваЁж ­ҐЇаҐалў­®© ¬®¤Ґ«Ё ®ЎкҐЄв ');
      write  ('  ђ §¬Ґа­®бвм ¬ ваЁжл  A:    n=');
      readln(n);
      writeln('  Њ ваЁж  A');
      EN(8,n,n,A);
      writeln('  Њ ваЁж  B');
      EN(9+n,n,1,B);
      writeln('  Њ ваЁж  C');
      EN(10+2*n,1,n,C);
      writeln('   Њ ваЁж  гЇа ў«пҐ¬®бвЁ');
      CONT(n,A,B,Sco);
      EX(n,n,Sco);
      DET(n,Sco,dc);
      writeln('   det Sco =',dc);
      ReadLn;
      writeln('   Њ ваЁж  ­ Ў«о¤ Ґ¬®бвЁ');
      OBS(n,A,C,Sob);
      EX(n,n,Sob);
      DET(n,Sob,dh);
      writeln('   det Sob =',dh);
      ReadLn;

      writeln('    „ЁбЄаҐвЁ§ жЁп ­ҐЇаҐалў­ле га ў­Ґ­Ё© б®бв®п­Ёп');
      write('ЏҐаЁ®¤ ¤ЁбЄаҐв­®бвЁ  tk=');
      readln(tk);
      SMP(tk,n,1,A,B,Ad,Bd);
      writeln;
      writeln('     „ЁбЄаҐв­ п  ¬®¤Ґ«м');
      writeln('   ¬ ваЁж  Ad');
      EX(n,n,Ad);
      writeln('   ¬ ваЁж  Bd');
      EX(n,1,Bd);
      ReadLn;
      writeln('   Њ ваЁж  гЇа ў«пҐ¬®бвЁ');
      CONT(n,Ad,Bd,Sco);
      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);
      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('   ‚Ґб®ў п ¬ ваЁж  Q');
      EN(5,n,n,Q);
      writeln('   ‚Ґб®ў п ¬ ваЁж  L');
      write(' l1=');
      readln(l);
      RICC(n,l,Ad,Bd,Q,P,Rsd);
      writeln('    ђҐиҐ­ЁҐ га ў­Ґ­Ёп ђЁЄЄ вЁ: ¬ ваЁж   P');
      EX(n,n,P);
      readln;
      writeln('    Њ ваЁж  ॣг«пв®а  б®бв®п­Ёп  Rsd ');
      EX(1,n,Rsd);
      readln;
      OBST (n,Ad,C,Hx);
      writeln('    Њ ваЁж  ®Ўа в­®© бўп§Ё ­ Ў«о¤ вҐ«п б®бв®п­Ёп  Hx ');
      writeln('     б Є®­Ґз­л¬ ўаҐ¬Ґ­Ґ¬ ЇҐаҐе®¤­®Ј® Їа®жҐбб         ');
      writeln;
      EX(n,1,Hx);
      readln;
      writeln('               ‡ ¬Є­гв п бЁб⥬                     ');
      writeln;
      writeln('   Њ ваЁж  § ¬Є­гв®© бЁб⥬л б ॣг«пв®а®¬ б®бв®п­Ёп');
      writeln('                  Acr=Ad-Bd*Rsd                     ');
      CLSYR(n,Ad,Bd,Rsd,Acr);
      EX(n,n,Acr);
      readln;
      writeln('   ‘®Ўб⢥­­ п ¬ ваЁж  ­ Ў«о¤ вҐ«п       ');
      writeln('                  Ach=Ad-Hx*C            ');
      CLSYH(n,Ad,C,Hx,Ach);
      EX(n,n,Ach);
      readln;
      readln;
      writeln('   Њ ваЁж  § ¬Є­гв®© бЁб⥬л б ॣг«пв®а®¬ б®бв®п­Ёп');
      writeln('                Ё ­ Ў«о¤ вҐ«Ґ¬                      ');
      writeln('               Ao=Ю Acr   Bd*Rsd Ю                  ');
      writeln('                  Ю  0      Ach  Ю                  ');
      writeln;
      CLS(n,Ad,Bd,C,Rsd,Hx,Ao);
      EX(2*n,2*n,Ao);
      readln;
      ClrScr;
      writeln('          ђҐиҐ­ЁҐ га ў­Ґ­Ё© б®бв®п­Ёп ');
      writeln;
      writeln('     ‚ҐЄв®а ­ з «м­ле гб«®ўЁ©    ');
      ENV(2*n,Xo);
      readln;
      writeln('                  _                 _                 _               ');
      writeln('  k    і x1(k)  і x1(k)  і x2(k)  і x2(k)  і x3(k)  і x3(k)  і u(k)  і');
      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,Rsd,_X,U);
      U[1]:=-U[1];
      nr:=0;
      for i:=1 to 2*n do
      begin
        Xo[i]:=X[i];
        nr:=nr+sqr(X[i]);
      end;
      nr:=sqrt(nr);
      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 nr<0.001;
      readln;
      readln;
    END.