Интерполяция табличной функции по методу Ньютона

Содержание
Содержание
Задание на курсовую работу.
Интерполяция по Ньютону
Разработка алгоритма решения задачи.
Организация ввода и вывода данных.
Организация контроля ошибочных действий пользователя.
Список литературы.
Приложение 1. Текст кода программы.
Приложение 2. Руководство пользователя программы.

Интерполяция табличной функции по методу Ньютона

begin
inc(n);
SetLength(Nodes,N);
Nodes[N-1].X := X;
Nodes[N-1].Y := Y;
CalcDeltasFla := false;
end;
function TPolynome.CalculateDelta(var eri,erj:integer):boolean;
var i,j:integer;
dx,dy:real;
begin
SetLength(Deltas,N);
Result := false;
erj := 0;
eri := 0;
if N=0 then exit;
for I := 0 to N - 1 do
Deltas[i] := Nodes[i].Y;
for I := 1 to N - 1 do
begin
for J := N - 1 downto I do
begin
dx := Nodes[j].X-Nodes[j-i].X;
if dx = 0 then
begin
// Фиксируем ошиюку
eri := j;
erj := j-I;
Result :=false;
CalcDeltasFla := false;
exit;
end;
// Считаем разделённую разность i-го порядка, для узла j-i
dy := Deltas[j]-Deltas[j-1];
Deltas[j] := dy / dx;
end;
end;
Result:= true;
CalcDeltasFla:= true;
end;
function TPolynome.CalculateY(const x: real; var IsOk:boolean): real;
var i,j:integer;
P:real;
ty:real;
begin
if (CalcDeltasFla or CalculateDelta(i,j)) then
begin
P:=1;
ty:=deltas[0];
for I := 1 to N - 1 do
begin
P := P*(x-Nodes[i-1].x);
ty := ty + P*deltas[i];
end;
Result := ty;
end
else
begin
IsOk := false;
Result := 0;
exit;
end
end;
procedure TPolynome.Clear;
begin
N := 0;
SetLength(Deltas,N);
SetLength(Nodes,N);
CalcDeltasFla := false;
end;
constructor TPolynome.Create;
begin
Clear;
end;
function TPolynome.GetNode(const i: integer): TNode;
begin
if (i<n) then Result := Nodes[i]
else
begin
Result.X := 0;
Result.Y := 0;
end;
end;
end.
Модуль главной формы main.pas.
unit main;
interface
uses
Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
Dialogs, Grids, StdCtrls, ExtCtrls, NewInt;
type
TfrmMain = class(TForm)
dgrdMain: TStringGrid;
btnSelFile: TRadioButton;
Label1: TLabel;
btnImmed: TRadioButton;
edtCount: TEdit;
Label2: TLabel;
Label3: TLabel;
imgWork: TImage;
edtFilename: TEdit;
Label4: TLabel;
btnFileName: TButton;
OD: TOpenDialog;
btnLoad: TButton;
btnCalculate: TButton;
btnExit: TButton;
btnClear: TButton;
procedure FormCreate(Sender: TObject);
procedure btnFileNameClick(Sender: TObject);
procedure FormDestroy(Sender: TObject);
procedure btnExitClick(Sender: TObject);
procedure btnLoadClick(Sender: TObject);
procedure btnSelFileClick(Sender: TObject);
procedure btnImmedClick(Sender: TObject);
procedure edtCountExit(Sender: TObject);
procedure btnCalculateClick(Sender: TObject);
procedure btnClearClick(Sender: TObject);
private
{ Private declarations }
pln:TPolynome;
public
{ Public declarations }
end;
var
frmMain: TfrmMain;
implementation
{$R *.dfm}
procedure TfrmMain.btnFileNameClick(Sender: TObject);
begin
if Od.Execute then
begin
edtFilename.Text :=Od.FileName;
end;
end;
procedure TfrmMain.btnImmedClick(Sender: TObject);
begin
dgrdMain.Options := dgrdMain.Options + [goEditing];
edtCount.Enabled:= true;
btnFileName.Enabled:= false;
edtFilename.Enabled := false;
btnLoad.Enabled := false;
end;
procedure TfrmMain.btnLoadClick(Sender: TObject);
var f : TextFile;
n:integer;
x,y:real;
I,j: Integer;
FileName : string;
begin