Технический форум - Методы нахождения корня на паскале

Технический форум (http://www.tehnari.ru/)

- Помощь студентам (http://www.tehnari.ru/f41/)

Методы нахождения корня на паскале

Здравствуйте. Помогите, пожалуйста разобраться с программой. Дано уравнение f(x)=1+sin(x)-1.14*exp(-x). Найти наименьший корень уравнения с точностью е=0.00005, используя методы половинного деления, ньютона и простой итерации. В принципе, я проги написала. Работают. Заминка только в конце: нужно вывести х0, х1, х2,...,хn для каждого метода. То есть вывести эти приближенные корни, которая прога считает в процессе и мне выдает окончательный результат. Как их вывести? подскажите, пожалуйста. Очень сильно нужно.

Метод половинного деления:

Код:

uses Crt;

function f(x:real):real;

begin

   f:=1+sin(x)-1.14*exp(-x);

end;

var x,eps,a,b,c,d:real;

begin

clrScr;

writeln('vvedite a  b');

readln(a,b);

writeLn('vvedite tochnost eps');

readln(eps);

repeat

c:=(a+b)/2;

if(f(a)*f(c))<0 then b:=c

else a:=c;

until (b-a)<=eps;

x:=(a+b)/2;

d:=(b-a)/2;

writeLn('koren x=',x:0:7, ' d=' ,d:0:7);

readln

end.

Метод Ньютона:

Код:

uses crt;

function F(x:real):real;

begin

F:=1+sin(x)-1.14*exp(-x);

end;

function F1(x:real):real;

begin

F1:=cos(x)+1.14*exp(-x);

end;

function Newton(x1,e:real):real;

var x2,b:real;

begin

x2:=x1;

repeat

b:=x2;

x2:=b-F(b)/F1(b);

until abs(x2-b)>e;

Newton:=x2;

end;

var a,eps:real;

begin

clrscr;

write('Введите начальное приближение a=');

readln(a);

write('Введите точность eps=');

readln(eps);

write('X=',Newton(a,eps):0:5);

readln

end.

Метод простой итерации:

Код:

var a,e,fa,fas,xa,xp,l,fx,fxp:real;

    n:integer;

function FunX(b:real):real;

begin

     fx:=1+sin(b)-1.14*exp(-b);

     fxp:=cos(b)+1.14*exp(-b);

     end;

begin

     write ('vvedite tochnost_');

     readln (e);

     write ('vvedite nachalnoe priblizenie x0_');

     readln (xp);

     n:=0;

     repeat

           xa:=xp;

           fx:=1+sin(xa)-1.14*exp(-xa);

           fxp:=cos(xa)+1.14*exp(-xa);

            l:=(-1)/(fxp+1);

           xp:=xa+l*fx;

           n:=n+1;

     until abs(xp-xa)<e;

     writeln ('otvet: x=', xp:9:4);

     writeln ('kol-vo provedennix vichislenii=',n);

     readln;

end.

В каждом листинге перед оператором until вставьте оператор вывода вычисленного значения.

т.е. через repeat until? можно пример, пожалуйста)

к примеру, в первом листинге вставьте перед строкой с until следующую команду:

Код:

writeln(c);

а как результат проги скопировать в блокнот допустим?

Сделайте вывод в текстовый файл и потом его распечатайте. Также можно реализовать вывод непосредственно на печатающее устройство прямо из программы, указав в операторе Writeln в качестве параметра устройство печати.

Способы подробно описаны в справочной системе среды разработки :)

еще вопрос, в первой проге он выдает числа в конце с Е, а можно их как-то в привычном виде выдать, как-нибудь примерно до 7 знака

Организация вывода в файл:

в разделе переменных определяем файловую переменную
var ....
F: Text;
....

В начале исполнимого модуля связываем файловую переменную с текстовым файлом, указывая полный путь к файлу и открываем его для перезаписи (или создания, при первом обращении к файлу):
begin
...
assign(F,'путь');
rewrite(F);
...

В операторах вывода результатов вычислений в качестве параметра устройства вывода указываем наш файл:

Writeln(F,'text');
Writeln(F,c);

В конце программы (перед конечным оператором end.) закрываем файл командой close(F);

Форматированный вывод значений можно организовать, используя параметры вывода значений. К примеру, форма записи C:4:8 определяет формат вывода числа с по крайней мере 4-мя знакоместами под целую часть и 8 значащими разрядами дробной части. Это используется для организации вывода значений в несколько выравненных по левому краю столбцов и, естественно, для избавления от экспоненциальной формы числа :).

Writeln(F,c:6:7); как пример к вашему вопросу