лаба №1 крэс-05 / shnyaga

Program Matr;
 uses crt;
 const n=8;
       n1=7;
       x:array [0..n] of real=(-4, -3, -2, -1, 0, 1, 2, 3, 4);
       y:array [0..n] of real=(0,   2,  0, -1, 0, 1, 0,-2, 0);
 var m:array [0..n] of real;
     a:array [1..7,1..8] of real;
     i,j:integer;

 procedure Holess;
  var m,t,i1,k:integer;
      max,p,s:real;
      a1:array [0..n] of real;
   begin
    for m:=1 to n1 do
     begin
      if m=1 then
       for j:=2 to n1+1 do a[1,j]:=a[1,j]/a[1,1];
      if m>=2 then
       begin
        max:=abs(a[m,m]);
        for t:=m to n1 do
         if abs(a[t,m])>=max then
          begin
           max:=abs(a[t,m]);
           i1:=t;
          end;
        for i:=1 to n1+1 do a1[i]:=a[m,i];
        for i:=1 to n1+1 do
         begin
          a[m,i]:=a[i1,i];
          a[i1,i]:=a1[i];
         end;
        for i:=m to n1 do
         begin
          p:=0;
          for k:=1 to m-1 do p:=p+a[i,k]*a[k,m];
          a[i,m]:=a[i,m]-p
         end;
        for j:=m+1 to n1+1 do
         begin
          s:=0;
          for k:=1 to m-1 do s:=s+a[m,k]*a[k,j];
          a[m,j]:=(a[m,j]-s)/a[m,m];
         end;
       end;
     end;
   end;

 procedure Res;
  var i,k:integer;
      s:real;
   begin
    for i:=n1 downto 1 do
     begin
      if n1=1 then x[i]:=a[n1,n1+1];
      s:=0;
      for k:=i+1 to n1 do
       s:=s+a[i,k]*x[k];
      x[i]:=a[i,n1+1]-s;
      writeln('M[',i,']=',x[i]:6:5);
     end;
    writeln(x[4]+4*x[5]+x[6]:6:5);
   end;

  begin
   clrscr;
   m[0]:=(-11*y[0]+18*y[1]-9*y[2]+2*y[3])/6;
   m[n]:=(11*y[n]-18*y[n-1]+9*y[n-2]-2*y[n-3])/6;
   for i:=1 to 7 do
    for j:=1 to 8 do a[i,j]:=0;
   for i:=1 to 7 do a[i,i]:=4;
   for i:=1 to 6 do
    begin
     a[i,i+1]:=1;
     a[i+1,i]:=1;
    end;
   a[1,8]:=3*(y[2]-y[0])-m[0];
   a[7,8]:=3*(y[8]-y[6])-m[8];
   for i:=2 to 6 do a[i,8]:=3*(y[i+1]-y[i-1]);
   for i:=1 to 7 do
    begin
     for j:=1 to 8 do write(a[i,j]:6:5,' | ');
     writeln;
    end;
   Holess;
   Res;

   {

   for i:=1 to 7 do
    begin
     for j:=1 to 8 do write(a[i,j]:6:5,' | ');
     writeln;
    end;}
   writeln(-2.06599+4*(-0.90034):6:5);
   readln;
  end.