Нелинейное уравнение методом хорд (подсказать)

[Swex]

Решить нелинейное уравнение методом хорд.

Вот мой код:

PHP код:

unit Unit1;

interface

uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, ComCtrls, ExtCtrls, TeeProcs, TeEngine, Chart, Math, Series;

type
  TForm1 = class(TForm)
    PageControl1: TPageControl;
    TabSheet1: TTabSheet;
    Edit1: TEdit;
    Edit2: TEdit;
    Edit3: TEdit;
    Edit4: TEdit;
    Memo1: TMemo;
    Button1: TButton;
    Button2: TButton;
    Chart1: TChart;
    Series1: TLineSeries;
    Edit5: TEdit;
    procedure Button1Click(Sender: TObject);
    procedure Button2Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;

var
  Form1: TForm1;

implementation

{$R *.dfm}

function f(x,a:extended):extended;
begin
f := a*power(cos(abs(x)),3)+a*sin(a*x)+a*power((sin(x)/cos(x)),4)+a*power(ln(x),2)+a;
end; 

function pp(x,a:extended):extended;
begin
pp := power((-a),3)*sin(a*x);
end; 


procedure TForm1.Button1Click(Sender: TObject);
var j:extended;
var x,x0,e,c,x1,x2,a:extended;
begin
e:= strtofloat(Edit3.Text);
X1 := strtofloat(Edit1.Text);
X2 := strtofloat(Edit2.Text);
a := strtofloat(Edit5.Text);
x := x1;
c := x2;
if f(c,a)*pp(a,c)<0 then begin c := x1; x := x2 end;
repeat
  x0 := x;
  x := (c*f(a,x0)-x0*f(a,c))/(f(a,x0)-f(a,c));
until abs(f(x,a))/0.00001<=e;
j := x;
Memo1.Lines.Add('Корень для метода хорд: '+floattostrf(j,ffgeneral,5,3));
end;
end. 


Но вылазит ошибка!!! Подскажите что исправить...

[Cannibal]

Где Ошибка? Какая ошибка?

[Swex]

Ошибка с таким содержанием: invalid floating point operation. В строчке, где описывается функция f. Я так понял, что она ругается на ln, который не модет быть 0 и меньше...но что можно исправить?

[Вадим Мошев]

Сообщение от Swex

Ошибка с таким содержанием: invalid floating point operation. В строчке, где описывается функция f. Я так понял, что она ругается на ln, который не модет быть 0 и меньше...но что можно исправить?

Нет, если такая функция задана, то ничего менять не надо, просто надо не допускать, чтобы функция принимала такие значения, при которых она не определена.

[Swex]

А как этого добиться, чтобы не допускала? промежуток я брал от 1 и дальше...