Решение нелинейных уравнений различными методами

[romaric]

Есть задача решить нелинейное уравнение на Pascal различными методами:
методом итерации, методом половинного деления, методом ньютона.
Можете подсказать что-нибудь?
e^x+ln x-10 x=0
Отрезок содержащий корень: [3;4]

[Вадим Мошев]

Что, например, подсказать? Пробуйте решать самостоятельно, если будут вопросы, задавайте, вероятно, вам смогут помочь.

[romaric]

Попробовал, методом половинного деления:

uses crt;
function F(x,y:real):real;
begin
F:=(exp(ln(y)*x))+(Ln(x))-(10*x);
end;
var a,b,c,e,y,x:real;
begin
clrscr;
a:=3;{левая граница}
b:=4;{правая граница}
y:=3;
e:=0.000001;{точность}
if F(a,y)*F(b,y)<0 then
  begin
   repeat
   c:=(a+b)/2;
   if F(a,y)*F(c,y)<=0 then b:=c
   else a:=c;
   until b-a<e;
   x:=(a+b)/2;
  end
else writeln('На данном интервале корней нет, либо взят большой интервал.');
writeln('Z(x,y)=',x:0:7);
readln
end.

Результат: Z(x,y)=3.0885139

В условии задачи была оговорка, что шаг следует считать так:
i=o,n
Xi=X0+ih
не могу понять, как это изобразить в pascal

[Вадим Мошев]

Прежде всего, по поводу офорления кода. Почитайте это: http://www.programmersforum.ru/showt...23#post1497723

[romaric]

Что не так? [СODE][/СODE] вписано

[Вадим Мошев]

function F(x,y:real):real;
begin
F:=(exp(ln(y)*x))+(Ln(x))-(10*x);
end;

Зачем так?
Вашу функцию можно записать так:

function F(x:real):real;
begin
F:=e(x)+ln(x)-exp(x*ln(10))
end;

Сообщение от romaric

Что не так? [СODE][/СODE] вписано

А вы прочитайте, как правильно пользоваться этим тегом

[romaric]

Не подумал. А по существу будет что? Или будет только набивание ЧСВ модератору, методом тыканья новичка в мелкие ошибки?

[Вадим Мошев]

В условии задачи была оговорка, что шаг следует считать так:
i=o,n
Xi=X0+ih
не могу понять, как это изобразить в pascal

Зачем вам считать шаг? Согласно этому методу, границы движутся навстречу друг другу.

Остальное вроде всё верно. Что проверить, насколько адекватно считает программа, вычислите значение функции, от того, что было найдено в процессе поиска корня.
Кстати

until b-a<e;

Надо писать:

until abs(b-a)<e;

А по существу будет что? Или будет только набивание ЧСВ модератору, методом тыканья новичка в мелкие ошибки?

А вот это зачем? Я вам не хамил.

[romaric]

Теперь по существу. Благодарю. Беру свои слова про ЧСВ обратно.

Решал методом итерации, код:

uses crt;
function F(x,y:real):real;
begin
F:=(exp(ln(y)*x))+(ln(x))-(10*x);
end;
{var x0,x1,a,b,e:real;}
var x0,x1,e:real;
iteraz:integer;
begin
clrscr;
write('Введите приближённое значение X=');
readln(x1);
write('Введите точность e=');
readln(e);
iteraz:=0;
repeat
iteraz:=iteraz+1;
x0:=x1;
x1:=F(x0);
until (abs(x1-x0)<=e);
writeln('Решение уравнения:');
writeln('Точное значение корня.........0.7376');
writeln('Вычисленное значение корня...',x1:6:5);
writeln('Число итераций.............. ',iteraz);
writeln('Программа закончена, нажмите Enter.');
readln;
end.

Программа не компилируется. Выдает: (18,5) Error: Wrong number of parameters specified for call to "fun"
Не могу врубится где ошибка. Точнее в чем именно.

[Вадим Мошев]

function F(x,y:real):real;

Потому что написанная вами функция ожидает два параметра, а вы передаёте ему один:

x1:=F(x0)