1 var a,b:array[0..100,0..100] of extended; 2     temp,ans:array[0..100] of extended; 3     i,j,k,n:longint; 4     cnt:extended; 5 begin 6     readln(n); 7     for i:=1 to n+1 do 8         for j:=1 to n do 9             read(b[i,j]);10     for i:=1 to n do11         for j:=1 to n do12         begin13             a[i,j]:=2*(b[i+1,j]-b[i,j]);14             a[i,n+1]:=a[i,n+1]+b[i+1,j]*b[i+1,j]-b[i,j]*b[i,j];15         end;16 //--------------------------------rpCardinal Orz--------------------------------------17     for i:=1 to n-1 do18     begin19         for j:=i+1 to n do20             if a[j,i]>a[i,i] then21             begin22                 temp:=a[i];23                 a[i]:=a[j];24                 a[j]:=temp;25             end;26         for j:=i+1 to n do27         begin28             cnt:=a[j,i]/a[i,i]; a[j,i]:=0;29             for k:=i+1 to n+1 do30                 a[j,k]:=a[i,k]*cnt-a[j,k];31         end;32     end;33     ans[n]:=a[n,n+1]/a[n,n];34     for i:=n-1 downto 1 do35     begin36         for j:=i+1 to n do a[i,n+1]:=a[i,n+1]-ans[j]*a[i,j];37         ans[i]:=a[i,n+1]/a[i,i];38     end;39     for i:=1 to n do write(ans[i]:0:3,‘ ‘);40 end.