o:=x0 + j * h;:=r * (xz);:=Result + r;;; TForm1.PolyGauss2: string;:=floattostr (y [0]); i:=1 to n do:=Delta (i, y, - (succ (i) div 2)) / (fact (i) * IntPower (h, i)); z=0 then continue ; (z> 0) then Result:=Result + «+»; z <> 1 then Result:=Result + floattostr (z); j:=- (i div 2) to (i div 2) -1 + i mod 2 do:=x0 + j * h; z=0Result:=Result + «x» if z> 0 then Result:=Result+«(x-»+floattostr(z)+«)»Result:=Result+«(x»+floattostr(z)+«)»;;;Tform1.PolyG2(x:extended):extended;:=(y[0]);i:=1 to n do:=Delta (i, y, - (succ (i) div 2)) / (fact (i) * IntPower (h, i)); z=0 then continue;:=z; j:=-(i div 2) to (i div 2) - 1 + i mod 2 do:=x0 + j * h;:=r * (xz);:=Result + r;;; TForm1.PolyStirling: string;:= floattostr (y [0]); i:=1 to n do:=(Delta (i, y, - (i div 2)) + Delta (i, y, - (succ (i) div 2))) / (2 * fact (i) * IntPower (h, i)); not (odd (i)) then z:=z / 2; z=0 then continue; (z> 0) then Result:=Result + « + »; z <> 1 then Result:=Result + floattostr (z); j:=- (i div 2) +1- i mod 2 to (i div 2) - 1 + i mod 2 do:= x0 + j * h; z=0Result:=Result + «x» if z> 0 then Result:=Result+«(x-»+floattostr(z)+«)»Result:=Result+«(x»+floattostr(z)+«)»;not(odd(i)) then:=- 2 * x0; z=0 then Result:=result+«2x»z<0Result:=Result+«(2x»+floattostr(z)+«)»Result:=Result+«(2x+»+floattostr(z)+«)»;;;;Tform1.PolyS(x:extended):extended;:=(y[0]);i:=1 to n do:=(Delta (i, y, - (i div 2)) + Delta (i, y, - (succ (i) div 2))) / (2 * fact (i) * IntPower (h, i)); not (odd (i)) then z:=z / 2; z=0 then continue;:=z; j:=- (i div 2) +1- i mod 2 to (i div 2)- 1 + i mod 2 do:=x0 + j * h;:=r * (xz); not (odd (i)) then r:=r * (2 * (x-x0));:=Result + r;;; TForm1.Button1Click (Sender: TObject);:=strtofloat (EditXo.text);:=(Strtofloat (EditXn.Text)-x0) / n; i:=0 to n do [i]:=strtofloat (setka.Cells [i +1,1]);:=PolyNuton1;:=PolyNuton2; i:=0 to n do [i-(n div 2)]:=y [i];:=x0 + (n div 2) * h;:=PolyGauss1;:=PolyGauss2;:=PolyStirling; Polynoms do [1,0]:=PN1; Cells [1,1]:=PN2; [1,2]:=PG1; Cells [1,3]:=PG2; [1,4]:=PS;;;; TForm1.FormCreate (Sender: TObject);:=4; Polynoms do [1]:=1593;
Cells [0,0]:=«I??????? ??????? »; [0,1]: =« II??????? ??????? »; [0,2]: =« I??????? ?????? »; [0,3]: =« II??????? ?????? »; [0,4]: =«??????? ????????? »;
end; setka do [0,0]:=« x »; [0,1]: =« y »;;; TForm1.EditXnExit (Sender: TObject); StrTofloat (EditXn.Text)
begin;; .5. Результати численних експериментів.
Добре обумовлені матриці (Введення користувачем)
№ Розмірність матріциВектор невязки (r) Норма матриці | | r | |-равномерное13r1=0,00001 r2=0 r3=0 | | r | |=10124r1=0 r2=0,00002 r3= 0,00003 r4=0,00003 | | r | |=153310r1=0,00002...
