Метод деления отрезка по полам, проблема в уравнении Опции

[TS*]

Ошибка в решении уравнения методом деления отрезка по полам(быстрее всего в записи уравнения), не могу найти корень, , eps ввожу 0.1, код програмы стандартный:

Uses Crt;
Function f(x: real): real;
 Begin
      f :=0.95*x-23.04*exp(0.33*ln(x))+46.63*exp(0.21*ln(x));
 End;
Var
 x, Eps, a, b, c: real;
 n: Integer;
begin
 ClrScr;
 Writeln('vvedite a i b'); Read(a, b);
 WriteLn('vvedite Eps'); Read(Eps);

 n := 0;
 Repeat
   c := (a + b) / 2;
   If (f(a) * f(с)) < 0 Then b := c
   Else a := c;
   Inc(n)
 Until (b - a) <= Eps;

 x := (a + b) / 2;
 WriteLn('корень равен  x=', x:10:7);
 WriteLn('количество итераций = ',n);
 readln;
end. 

Хотя остальные уравнения без выражения типа "exp(a*ln(x))" решаються без проблем...

Сообщение отредактировано: volvo - 3.11.2006 0:07

[volvo]

Вот так придется заменять функцию возведения в степень (найдено у нас на форуме... Здесь), и заданную тебе функцию:

Function power(a, b: real): real;

  function relerror(observed, actual: real): real;
  begin
    if actual = 0.0 then relerror := abs(observed)
    else relerror := abs(observed / actual - 1);
  end;

var
  i: byte;
  res: real;
begin
  if (trunc(b) = b) and (b > 0) then begin
    res := 1;
    for i := 1 to trunc(b) do res := res * a;
    power := res; exit
  end;

  if a = 0 then
    if b = 0 then power := 1
    else
      if b < 0 then
      else power := 0
  else
    if a < 0 then
      if abs(b) < 1e-10 then power := 1
      else
        if relerror(b,round(b)) < 1e-8 then
          power := (1-2*ord(odd(round(b))))*exp(b*ln(abs(a)))
        else
          if (relerror(1/b, round(1/b)) < 1e-8) and odd(round(1/b))
          then power:=-exp(b*ln(abs(a)))
          else
    else
      power := exp(b*ln(a))
end;

{ Твоя функция F(X) = ... }
function f(x: real): real;
var
  res: real;
const eps = 0.00001;
begin
  res := 0.35*x-23.04*power(x, 0.33)+46.63*power(x, 0.21);
  if abs(res) < eps then res := - eps;
  f := res;
end;

[TS*]

Спасиба!