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

[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

[TS*]

f(x)=0.35*x-23.04*x^0.33+46.63*x^0.21

[Lapp]

f(x)=0.35*x-23.04*x^0.33+46.63*x^0.21

TS*, эта функция, если говорить строго, не имеет корней на области определения. Если расширить область определения значением х=0 по непрерывности, то есть корень х=0. Для этого нужно ввести в f условный оператор типа такого:
if x<1e-10 then f:=0 else f:=...

Если тебя интересуют подробности поведения этой функции, то задай вопрос в Математике. Я бы на твоем месте проверил, насколько она соответствует реальным данным (с графика). Кстати, программу интерполяции ты выложил, а вот сами данные - нет. А без них ничего не скажешь..

2 Volvo:
прога на той ссылке несколько странная.. На отрицательные величины она почему-то не ругается, а выдает числа (фантасмагорические). Идея замены возведения в степень на умножение сама по себе заслуживает внимания (и, кстати, допускает отрицательные значения аргумента при целых показателях) - может, в некоторых (явно не в этом) случаях это может оказаться полезным. Но именно эта функция все равно требует переделки..