
Добавил:
Upload
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Новая папка_2 / Program511
.pas PROGRAM NEW;
{$N+,E+}
uses CRT;
var
eps0:real; m1:real; m2:real; x0:real;
var
x:real; eps:real; k:integer;
function f(x:real):real;
begin
f:= -(x*x*x+3)/4
end;
function f1(x:real):real;
begin
f1:=-(3*x*x)/4
end;
procedure Newton (eps:real; var x:real; var k:integer);
var y:real;
begin
k:=0;
repeat
y:=x;
x:=y-f(y)/f1(y);
k:=k+1;
until abs(y-x)<eps
end;
begin
Writeln(' eps0 ');
Readln(eps0);
Writeln(' m1 ');
Readln(m1);
Writeln(' m2 ');
Readln(m2);
Writeln(' x0 ');
Readln(x0);
x:=x0;
eps:=eps0;
{eps:=sqrt(2*m1*eps0/m2); }
Newton(eps,x,k);
WriteLn('количество итераций = ', k:3);
WriteLn('приближённое значение корня = ', x:12:8);
end.
{$N+,E+}
uses CRT;
var
eps0:real; m1:real; m2:real; x0:real;
var
x:real; eps:real; k:integer;
function f(x:real):real;
begin
f:= -(x*x*x+3)/4
end;
function f1(x:real):real;
begin
f1:=-(3*x*x)/4
end;
procedure Newton (eps:real; var x:real; var k:integer);
var y:real;
begin
k:=0;
repeat
y:=x;
x:=y-f(y)/f1(y);
k:=k+1;
until abs(y-x)<eps
end;
begin
Writeln(' eps0 ');
Readln(eps0);
Writeln(' m1 ');
Readln(m1);
Writeln(' m2 ');
Readln(m2);
Writeln(' x0 ');
Readln(x0);
x:=x0;
eps:=eps0;
{eps:=sqrt(2*m1*eps0/m2); }
Newton(eps,x,k);
WriteLn('количество итераций = ', k:3);
WriteLn('приближённое значение корня = ', x:12:8);
end.
Соседние файлы в папке Новая папка_2