Ошибка в решении уравнения методом деления отрезка по полам(быстрее всего в записи уравнения), не могу найти корень, , 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.
Значения a и b какие вводишь? X не должен быть отрицательным ни на каком этапе...
Нулевым, кстати, тоже...
a=-100; b=100...
Ну, вот тебе и ответ...
Ln(-100) чему равен?
... понял, а можна как-то поменять выражения (Exp(а * Ln(x))?...
Ответ должен быть:
x=99.9969482
?
Если да, то покажу тебе, КАК решать...
Нет... Не пойдет... У этой функции ВООБЩЕ есть решения? Кстати, саму функцию приведи...
f(x)=0.35*x-23.04*x^0.33+46.63*x^0.21
Функция получена с помощью апроксимации на програмке <Numeri>,ввожу точки полученые из графика, а потом в функции пользователя нахожу коэфициенты.
Прикрепленные файлы
numeri.zip ( 263.58 килобайт )
Кол-во скачиваний: 221
Вот так придется заменять функцию возведения в степень (найдено у нас на форуме... http://forum.pascal.net.ru/index.php?s=&showtopic=987&view=findpost&p=9689), и заданную тебе функцию:
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;
Спасиба!