Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

3курс / WEGSTEIN

.PAS
Скачиваний:
14
Добавлен:
11.04.2015
Размер:
1.22 Кб
Скачать
Uses CRT;
const
n=5;{j}
h=0.2;
epsilon=0.0001;
type
koef=array[0..n] of real;
fun=function(n:integer;p:koef;t:real):real;
var
P:koef;
k: integer;
t,z: real;
ft:text;
zag:string;

function Polynom(n:integer;p:koef;t:real):real;
var
z: real;
k: integer;
begin
{ }
z:=P[0];
for k:=1 to n do
z:=z*t+P[k];
Polynom:=z;
end;
function phi(x:real):real;
begin
phi:=x+Polynom(n,p,x);
end;
procedure root(x0,x1:real;var x:real);
var
y1,y2,x2:real;
delta:real;
begin
y1:=phi(x0);y2:=phi(x1);

repeat
x2:=y2-(y2-y1)*(y2-x1)/((y2-y1)-(x1-x0));
writeln(x1,x2);
x0:=x1;x1:=x2;
y1:=y2;y2:=phi(x2);
delta:=abs(Polynom(n,P,x2));
until delta<epsilon;
x:=x2;
end;




begin
clrscr;
assign(ft,'c:\Root\in1.txt');
reset(ft);
readln(ft,zag);
writeln(zag);
writeln('koef ');
for k:=0 to n do
begin
readln(ft,P[k]);
writeln(P[k]);
end;
close(ft);
root(1,2,t);
assign(ft,'c:\Root\out1.txt');
rewrite(ft);
writeln(ft,t);
close(ft);

end.
Соседние файлы в папке 3курс