Необходимо исправить код программы, желательно перенести его в процедуру!

[leshka_good]

Есть вот такая часть кода программы, Sum считает скалярные произведения, в итоге получается const. Идет зависимость от степени многочлена, необходимо реализовать для n-й степени!

For j:=0 to n1 do begin     //mass P(t) & X^n1
For i:=0 to n do begin
                 P[i]:=StrToFloat(StringGrid1.Cells[n1+1,i+1]);
                 X[i,j]:=StrToFloat(StringGrid1.Cells[j,i+1]);
                 end;
                  end;


For j:=0 to n1 do begin
                  StringGrid1.Cells[j+2+n1,0]:='||H'+IntToStr(j)+'||';
                  StringGrid1.Cells[j+3+n1+n1,0]:='e'+IntToStr(j);
                  StringGrid1.Cells[j+4+n1+n1+n1,0]:='H'+IntToStr(j);

For i:= 0 to n-1 do            
              begin
              m:=0;
            Repeat
              For k:= 0 to j do begin
              For l:= 0 to n-1 do begin
                                  Sum[m,j]:=Sum[m,j]+X[l,k]*P[l]*h*E[l,k-1];                              
                                  end;
                                  Su[m,j]:=Sum[m,j];
                                  Sum[m,j]:=0;
                                  end;

              StringGrid1.Cells[j+4+n1+n1+n1,i+1]:=FloatToStr(X[i,j]-Su[m,j]*E[i,j-1]);
              StringGrid1.Cells[j+2+n1,i+1]:=FloatToStr(Power(X[i,j]-Su[m,j]*E[i,j-1],2)*P[i]*h);
              Hi[i,j]:=X[i,j]-Su[m,j]*E[i,j-1];
              C[j]:=C[j]+Power(X[i,j]-Su[m,j]*E[i,j-1],2)*P[i]*h;
              StringGrid1.Cells[j+2+n1,n+1]:=FloatToStr(C[j]);
              StringGrid1.Cells[j+2+n1,n+2]:=FloatToStr(sqrt(C[j]));
              m:=m+1;
            Until m<>n1;
            end;
For i:= 0 to n-1 do            //ei
              begin
              E[i,j]:=Hi[i,j]/sqrt(C[j]);
              StringGrid1.Cells[j+3+n1+n1,i+1]:=FloatToStr(E[i,j]);
              end;

end;