Скачать: CL | DM;

Помогите с оптимизацией   Найти похожие ветки 

← →
Dragon-zla ©   (2009-10-06 22:36) [0]

Здравствуйте, помогите пожалуйста с оптимизацией кода
Задача состоит в решении уравнения F(x)=0 на отрезке [a,b] с заданной точностью e методом деления отрезка пополам и выводом найденного решения на экран.
Точность = 1e-0,5 (0.00001)
Отрезок [-0.5;3]
Уравнение (1/(sqrt(5)+sin(0.1*x)+ln(1+x))-x)
Код я написал, он компилируется, запрашивает исходные данные, начинает выполнять код, а вот потом ОЧЕНЬ медленно приближается к заданной точности. Подскажите пожалуйста можно оптимизировать его, для ускорения вычисления или это мне с примером "повезло"?
P.S. пробовал изменить точность до 0.1 - бессмысленно
Собственно сам код:
program ecsAdd1_4;

{$APPTYPE CONSOLE}

uses
 SysUtils,crt32,Math;
const
 epsil = 0.00001;
var
 a,b,c,Fa,Fb,Fc : real;
begin
writeln(oem("Введите начальное значение отрезка"));
readln(a);
writeln(oem("Введите конечное значение отрезка"));
readln(b);
while abs(b-a)>epsil do
 begin
 if a*b<0 then
   c := (a+b)/2
 else
   c := (a-b)/2;
 Fa := (1/(sqrt(5)+sin(0.1*a)+ln(1+a))-a);
 Fb := (1/(sqrt(5)+sin(0.1*b)+ln(1+b))-b);
 Fc := (1/(sqrt(5)+sin(0.1*c)+ln(1+c))-c);
 if Fa*Fc<0 then
 b := c
 else
 a := c;
 end;
Fa := (1/(sqrt(5)+sin(0.1*a)+ln(1+a))-a);
Fb := (1/(sqrt(5)+sin(0.1*b)+ln(1+b))-b);
if Fa*Fb<0 then begin
   if a*b<0 then
     c := (a+b)/2
   else
     c := (a-b)/2;
Fc := (1/(sqrt(5)+sin(0.1*c)+ln(1+c))-c);
writeln(oem("Приблизительный корень уравнения = ") ,c);
writeln(oem("Значение функции в этой точке = ") ,Fc);
readkey;
end
else
writeln(oem("Корень не найден"));
readkey;
end.

---

← →
palva ©   (2009-10-06 23:52) [1]

if a*b<0 then
  c := (a+b)/2
else
  c := (a-b)/2;

Можете объяснить на словах, что это такое?

---

← →
Омлет   (2009-10-07 00:11) [2]

У тебя неправильная реализация.
Первая ссылка в яндексе (даже с кодом): http://eco.sutd.ru/Study/Informat/mpd.html

---

Страницы: 1 вся ветка

Скачать: CL | DM;



Память: 0.45 MB
Время: 0.011 c