итерационные циклы

program lab2;
uses crt;
const eps = 0.00001;
pi = 3.14159265358979;
var i,m:byte;
s,x,xn,xk,hx,max,n,a,b,c,k1,k2,y,sy:real; {d=summa riada}
ch:char;
begin
clrscr;
writeln ('enter star and final means of xn,xk,cross of change func');
readln(xn,xk,hx);
writeln(' xn= ',xn:5:2,' xk= ', xk:5:2, ' hx= ', hx:5:2);
writeln('г===========T============T==========T==========T==========¬');
 
m:=round((xk-xn)/hx)+1;
x:=xn;
k1:=0.5*sqrt(5+sqrt(5));
k2:=(sqrt(5)-1)/2;
 
for i:=1 to m do
a:=1; b:=x; c:=-1; s:=-1*(k1-x*k2); n:=0;
while abs(s)>eps do
 
begin
n:=n+1;
a:=a*x*x/(2*n+1)*(2*n+2);
b:=b*x*x/(2*n+2)*(2*n+3);
c:=-c;
s:=c*(a*k1-b*k2);
end;
s:=s/2;
x:=x+hx;
y:=cos(x+pi/10);
if abs(y)>abs(s) then
max:=y
else
max:=s;
if (abs(y)>1E-5) or (abs(s)>1E-5) then
sy:=100*abs((y-s)/max)
else
sy:=0;
writeln('¦ x ¦ n ¦ s ¦ y ¦ sy,% ¦');
writeln('L===========¦============¦==========¦==========¦==========-');
writeln('¦ ',x:6:2,' ¦ ',n:3:0,' ¦ ',s:9:6,'¦', y:9:6,' ¦',sy:10:1,'¦');
 
writeln('¦============¦===========¦==========¦===========¦=========¦');
 
repeat
ch:=readkey;
until ord(ch)=27;
end.