@elluminate

пишу программу,решающую уравнение методом ньютона. уравнение вида 3ln^2(x)+ 6ln(x) - 5=0. начальное приближение нашел вручную.равно 2. написал код

unit Unit1;

interface

uses
Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
Dialogs, StdCtrls;

type
TForm1 = class(TForm)
Button1: TButton;
Edit1: TEdit;
Memo1: TMemo;

procedure Button1Click(Sender: TObject);
private
{ Private declarations }
public
{ Public declarations }
end;

var
Form1: TForm1;
i,n,t,result,result1:real;

implementation
function f(i:real):real;
begin
result:=3*sqr(ln(i))+6*ln(i)-5;
end;
function d(i:real):real;
begin
result1:=(6*ln(i)+6)/i;
end;

{$R *.dfm}







procedure TForm1.Button1Click(Sender: TObject);
var eps:real;
begin

eps:=0.0001;
i:=strtofloat(edit1.text);
n:=0;

repeat

i:=i-(f(i)/d(i));
memo1.lines.add('x'+floattostr(n)+'='+floattostr(i));
n:=n+1;
until abs(f(i))<eps;
end;



end.




компилируется.но возвращает ошибку деления на ноль.помогите.

0

@igoriy

Ошибка возникает на этой строке Почему не знаю,вроде деления здесь нет.
Можно сделать так:

0

@elluminate

попробовал.смысла-то нет.она как не работала,так и не работает.только сообщает об оштбке теперь =) это не выход

0

@igoriy

Я тут полазил по просторам нета . Синтаксис функции LN таков что Ln(0) - такое неззяя.А у тебя такое получается - в строке которую я показал.Она делает все действия до нуля и прекрасно заполняет memo.

0

@elluminate

я знаю что нельзя.так как исправить

0

@flatgad

просто когда будет 0 то пропусти операцию и продолжай дальше делай через if

0

@igoriy

ну или убери и все тебе же отрицательные то же не покатят

0

@flatgad

Да точно завтыкал i>0 надо

0

@elluminate

справился =)))) спасибо всем за помощь,правда ни один совет не работает. =D вот рабочий код


unit Unit1;

interface

uses
Windows, Math, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
Dialogs, StdCtrls;

type
TForm1 = class(TForm)
Edit1: TEdit;
Label1: TLabel;
Button1: TButton;
Label2: TLabel;
Edit2: TEdit;
procedure Button1Click(Sender: TObject);
private
{ Private declarations }
public
{ Public declarations }
end;

var
Form1: TForm1;
t,i:real;

implementation

{$R *.dfm}
function f(i:real):real;
begin
result:=3*sqr(ln(i))+6*ln(i)-5;
end;
function d(i:real):real;
begin
result:=(6*ln(i)+6)/i;
end;

procedure TForm1.Button1Click(Sender: TObject);
var eps:real;
begin
eps:=0.0001;
i:=strtofloat(edit1.text);

repeat
t:=f(i)/d(i);
i:=i-t;

until abs(f(i))<eps;
edit2.Text:=floattostr(i);


end;

end.



может кому-то пригодится =)

1

@simon140906

Поясните, пожалуйста, про переменную t. Я так и не увидел где задается её значение и какую роль она играет.

0