Метод Ньютона Опции

[Merhaba]

Добрый Вечер!!! Помогите Пожалуйста устранить ошибку!

program MetodNuytona;
const eps = 1e-6;
var x,x0: real;
Function f(x: real): real;
begin
f:=ln(exp(ln(2)*x)-5)/ln(3)+ln(exp(ln(2)*x)-7)/ln(3)-1;
end;
Function f1(x: real): real;
begin
f1:=exp(ln(2)*x)*ln(2)/((exp(ln(2)*x)-5)*ln(3)+exp(ln(2)*x)*ln(2))/((exp(ln(2)*x)-5)*ln(3));
end;
Begin
writeln;
writeln('Решение нелинейного уравнения f(x)=0');
writeln;
write('x0 = ');
read(x0);
writeln;
x:=x0;
repeat
x:=x-f(x)/f1(x);
until abs(f(x)/f1(x))<eps;
writeln('Метод Ньютона**:');
writeln;
writeln('x=',x:8:5);
End.

Free Pascal почему-то работает только с x0=3 (исключительно с ним)..
Выражение:

log_3(2^x-5)+log_3(2^x-7)=1

Где может быть ошибка?

Сообщение отредактировано: Merhaba - 14.12.2011 1:10

[IUnknown]

Значит, проверяй, правильно ли вычисляется производная в твоей программе. Похоже, что нет...

[Merhaba]

Значит, проверяй, правильно ли вычисляется производная в твоей программе. Похоже, что нет...

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