Метод Ньютона Delphi7

Описание программы

Осуществляется решение системы нелинейных алгебраических уравнений (СНАУ) методом Ньютона. С клавиатуры задается начальное приближение и точность решения.

Фрагмент кода из программы:

function Func(x:real):real;
function D_Func( x:real):real;
function Iter_Func(x,nju:real):real;
  implementation

{$R *.dfm}


//-----------------------------------------------------------------------------
function Func(x:real):real; //функция!
begin
   Func:=x*sqr(x)-3*x-2*exp(-x);
end;

//-----------------------------------------------------------------------------
function D_Func(x:real):real;//производная функции!
begin
  D_Func:=3*sqr(x)-3*x+2*exp(-x);
end;

//----------------------------------------------------------------------------
function Iter_Func(x,nju:real):real;
begin
  Iter_Func:=x-nju*Func(x);
end;


//------------------------------------------------------------------------------
// Метод по приведению к результату х2
function NetownsMetod(Interval:TInterval; eps:real; var Result2:real):integer;
var
  x1,x2,Verification:real;
  CountIter:integer;
begin
  x1:=Interval.x1;
  x2:=x1-Func(x1)/D_Func(x1);
  CountIter:=0;
    while abs(x1-x2)>eps do
    begin
      inc(CountIter);
      x1:=x2;
      x2:=x1-Func(x1)/D_Func(x1);
    end;

  Result2:=x2;
  NetownsMetod:=CountIter;
end;

Содержание архива: