vlad / TURBO P / tp7 / BIN / 2_2

Program Metod_Nutona;
uses crt;
var a,b,x,y1,y2,x1,x2,y,l,eps,k:real;
begin clrscr;
writeln('‚ўҐ¤iвм §­ зҐ­­п i­вҐаў «г [a,b]:');
readln(a,b);
writeln('‚ўҐ¤iвм в®з­iбвм:');
readln(eps);
x:=0.3;
y1:=(4*sqr(x1)-4*x1+1)/(x1-sqr(x1)+1);
y2:=(8*x*sqr(x1)-8*sqr(x1)+10*x1-5)/sqr(x1-sqr(x1)+1);
repeat
begin
x2:=x1-(y1/y2);
l:=abs(x2-x1);
y1:=(4*sqr(x2)-4*x2+1)/(x2-sqr(x2)+1);
y2:=(8*x2*sqr(x2)-8*sqr(x2)+10*x2-5)/(sqr(x2-sqr(x2)+1));
x1:=x2;
k:=k+1;
end;
until (l<eps) or (y1=0);
writeln('x=',x1:4:2);
y1:=(1-x1+sqr(x1))/(1+x1-sqr(x1));
writeln('Mi­i¬ «м­Ґ §­ зҐ­­п ҐЄбв६г¬г: y(x)=',y1:6:4,' k= ',k:2:0);
readln;
end.