найти значение полинома четвертого порядка

[trung]

здравствуйте
здравствуйте
помогите дополнять программу
текст программы:
Program Lab9;
Uses Crt;
Const n=8;
m=4;
x:array [1..n] of real=(1.3,1.5,2.4,3.8,5,6.8,7.1,8.2 );
y:array [1..n] of real=(0.3,2.2,6.4,5.9,1,-0.2,-1.3,0);
Var Mat:array [1..m+1,1..m+1] of real;
a,b:array [1..m+1] of real;
P,Nevaz:array [1..n] of real;
Sum,Sumnev:real;
i,j,k,s:integer;
Function Power(x:real;y:integer):real;
Var St:real;
i:integer;
begin
St:=1;
For i:=1 to y do
St:=St*x;
Power:=St;
end;

Begin
ClrScr;
For i:=1 to m+1 do
For j:=1 to m+1 do
For k:=1 to n do
Mat[i,j]:=Mat[i,j]+Power(x[k],i+j-2);

For i:=1 to m+1 do
For k:=1 to n do
b[i]:=b[i]+y[k]*Power(x[k],i-1);

s:=1;
Repeat
For j:=s+1 to m+1 do
begin
Mat[s,j]:=Mat[s,j]/Mat[s,s];
end;
b[s]:=b[s]/Mat[s,s];
i:=s+1;

Repeat
For j:=s+1 to m+1 do
begin
Mat[i,j]:=Mat[i,j]-Mat[i,s]*Mat[s,j];
end;

b[i]:=b[i]-b[s]*Mat[i,s];
i:=i+1;
Until i>m+1;

s:=s+1;
Until s>m;

writeln ('Otvet:');
a[m+1]:=b[m+1]/Mat[m+1,m+1];
Repeat
s:=s-1;
Sum:=0;
For k:=s+1 to m+1 do
Sum:=Sum+Mat[s,k]*a[k];
a[s]:=b[s]-Sum;
Until s<2;


For i:=1 to n do
begin
Sum:=0;
For j:=1 to m+1 do
Sum:=Sum+a[j]*Power(x[i],j-1);

P[i]:=Sum;
Nevaz[i]:=P[i]-y[i];
Writeln('P',i,'=',P[i]:6:3,' Nevyazka',i,'=',Nevaz[i]:6:3);
end;
Sumnev:=0;
For i:=1 to n do begin
Sumnev:=Sumnev+Nevaz[i];
end;
writeln('Summarnaya nevyazka=',Sumnev:6:14);
Readln;
end.



помогите дописать эта выражение в проге

P(x)=a*x+a*x^2+a*x^3+a*x^4