лабы / лаб 2 / unit1

unit Unit1;

{$mode objfpc}{$H+}

interface

uses
Classes, SysUtils, FileUtil, Forms, Controls, Graphics, Dialogs;

type
TForm1 = class(TForm)
private
{ private declarations }
public
{ public declarations }
end;

const m =5;
var
Form1: TForm1;
a:array [1..m,1..m] of double =((4,-1,0,0,0),
(1,4,-1,0,0),
(0,1,4,-1,0),
(0,0,1,4,-1),
(0,0,0,1,4));
b:array [1..m] of double = (2,6,10,14,24);
i,j,k,n:integer;d:double; x:array[1..m] of double;

implementation

{$R *.lfm}

begin

writeln('Enter n:');
readln(n);
for i:=0 to n-1 do
begin
for j:=0 to n-1 do
read(a[i,j]);
end;

for j:=0 to m-1 do
read(a[i,j]);

for i:=1 to n-1 do

begin

for k:=i+1 to n do

begin

d:=a[k,i]/a[i,i];

for j:=i+1 to n do

a[k,j]:=a[k,j]-a[i,j]*d;

b[k]:=b[k]-b[i]*d;

end;

end;

{................О Б Р А Т Ы Й Х О Д...........}

for i:=n downto 1 do

begin

for j:=i+1 to n do

b[i]:=b[i]-a[i,j]*x[j];

x[i]:=b[i]/a[i,i];

end;

writeln(' О Т В Е Т: ');

end.