ender: TObject) ;: integer; .Clear; .Clear; .Clear; .Step:=StrToFloat (LabeledEdit1.Text); i:=1 to Meth.N do.x [i]:=StrToFloat (StringGrid1.Cells [i, 1]) ;. y [i]:=StrToFloat (StringGrid1.Cells [i, 2]) ;; i:=0 to Meth.N - 1 do.x1[i]:=StrToFloat(StringGrid1.Cells[i+1,1]);.y1[i]:=StrToFloat(StringGrid1.Cells[i+1,2]);;RadioGroup1.ItemIndex =0 then.Graf; RadioGroup1.ItemIndex=1 then.Newton ;;
procedure TForm1.Button3Click (Sender: TObject); RadioGroup1.ItemIndex=0 then.Text:=FloatToStr(Meth.PointLag(StrToFloat(Edit2.Text))).Text:=FloatToStr(Meth.PointNew(StrToFloat(Edit2.Text)));
end;.
Код програми Паскаль по інтерполяційної формулою Лагранжа:
unit Unit2;
interface, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms ,, Grids, StdCtrls, Math;=array [1..10] of real;=array [0..9] of real;=array [1..10] of array [1..10] of real;=class: integer;, y, a: Tvector;, x1: Tvector1;: TMatr;: real ;: string;
//function Lagrange(q:real):real;Graf;PointLag(point:real):real;Razn();RealD(xk:real;x0:real):real;Factorial(Number:integer): integer; Newton (); PointNew (point: real): real ;;
implementationUnit1;
{function Methods.Lagrange (q: real): real;, j: integer;, s: real;:=0; i:=1 to n do:=1; j:=1 to N doj lt; gt; i then:=s * (qx [j])/(x [i] -x [j]);:=L + y [i] * s ;;. Memo1.Lines.Add ( x= + FloatToStr (q) + y= + FloatToStr (L));
//Form1.Series1.AddXY (q, l);:=l; }
procedure Methods.Graf () ;: array [1..10] of array [1..10] of real ;: array [1..10] of real ;: array [1..10] of real;, j, k: byte ;, step, count: real;, My: real; i:=1 to N do [i]:=y [i]; j:=1 to N do [ j, i]:=Power (x [j], i - 1);
end;
{Прямий хід - виключення змінних}
for k:=1 to n - 1 doi:=k + 1 to n do [i, k]:=- MassivA [i, k]/MassivA [k, k]; j: =k + 1 to n do[i,j]:=MassivA[i,j]+MassivA[i,k]*MassivA[k,j];;[i]:=MassivB[i]+MassivA[i,k]*MassivB[k];;[n]:=MassivB[n]/MassivA[n,n];
{Зворотний хід - знаходження коренів} k:=n - 1 downto 1 do
begin:=MassivB [k]; j:=k + 1 to n do:=h-MassivX [j] * MassivA [k, j]; [k]:=h/MassivA [ k, k] ;; i:=1 to n do [i]:=MassivX [i];:= raquo ;;
//побудова многочленаi:=0 to N - 1 do ((i=N - 1) or (i=0)) then (i=N - 1) then (a [Ni] lt ; 0) then:=Polinom + FloatToStr (a [Ni]):=Polinom + + + FloatToStr (a [Ni]);:=Polinom + FloatToStr (a [Ni]) + x ^ + IntToStr (Ni - 1); (a [Ni] lt; 0) then:=Polinom + FloatToStr (a [Ni]) + x ^ + IntToStr (Ni - 1):=Polinom + + + FloatToStr (a [Ni]) + x ^ + IntToStr (Ni - 1) ;;;. Memo2.Lines.Add (Polinom);
//побудова графіка:=x [1]; step lt;=x [N] do:=step;:=0; i:=1 to N do:=My + a [i] * Power (Mx, i - 1) ;;. Series1.AddXY (Mx, My);:=step + 0.001 ;;
//ущільнення таблиці:=x [1];:=StrToFloat (Form1.LabeledEdit1.Text); step lt;=x [N] do:=step;:=0; i:=1 to N do:=My + a [i] * Power (Mx, i - 1) ;;. Memo1.Lines.Add ( x= + Format (%. 3f , [Mx]) + y= + Format (%. 4f , [My]));:=step + count ;;;
function Methods.PointLag (point: real): real;: byte; i:=1 to N do:=result + a [i] * Power (point, i - 1) ;;
//кінцівок разностіMethods.Razn () ;, j: byte; j:=1 to N - 1 do [j, 1]:=y [j + 1] -y [j]; i:=2 to N doj:=1 to N-i + 1 do [j, i]:=mas [j + 1, i - 1] -mas [j, i - 1] ;;
//висновок DMethods.RealD (xk: real; x0: real): real;:=(xk-x0)/(x [2] -x [1]) ;;
function Methods.Factorial (Number: integer): integer;: integer;:=1; i:=1 to Number do:=Result * Number ;;
procedure Methods.Newton () ;, j: integer ;, Ndop, d: real ;, xk, point: real;:=x [1];:=StrToFloat (Form1.Edit5.Text );:=StrToFloat (Form1.Edit4.Text);
//d:=RealD (xk, xn) ;; (point lt;=xk) do:=y [1]; i:=1 to N - 1 do:=RealD (point, xn);:= d * mas [1, i]/Factorial (i); (i gt; 1) thenj:=1 to i - 1 do:=Ndop * (dj) ;;:=New + Ndop ;;. Series1.AddXY (point , New) ;. Memo1.Lines.Add ( x= + Format (%. 3f , [point]) + y= + Format (%. 4f , [New ]));:=point + Step ;;. Memo2.Lines.Add (Polinom) ;;
function Methods.PointNew (point: real): real ;, Ndop, d: real ;, j: byte;:=y [1]; i:=1 to N - 1 do:=RealD (point, x [1]);:=d * mas [1, i]/Factorial (i); (i gt; 1) thenj:=1 to i - 1 do:=Ndop * (dj) ;;:=New + Ndop; ;;
end.
