С точностью до [LATEX]\varepsilon = {10}^{-5}[/LATEX] вычислить приближенное значения S

uses crt;
var
i,n:integer;
x,e, S,y, Pogreshnost:real;
factorial: longint;
 
Begin
// Считываем ввод переменных
WRITELN('Введите x');
READ(x);
WRITELN('Введите N');
READ(n);
 
e:= 0.00001 ;
 
// находим y
 y:= power(((power(x,2)/4+x/2+1)),x)*power(e,2);
 
repeat
write('x = ', x);
 
// находим факториал числа
factorial := 1;
 
    for i:=2 to n do
        factorial := factorial * i;
// находим S
S:= S + (power(n,2)+1)/factorial*power((x/2),n);
// Pohubka
Pogreshnost := Abs((S-y)/y)*100;
// Выводим ответ на экран
WRITE('; S = ', S);
WRITE('; y = ', y);
WRITE('; Pogreshnost = ', Pogreshnost);
WRITELN();
 
n:=n+1;
 
until S<y;
 
End.

uses crt;
var e,a,b,x,h,s,y,p,t:real;
    n:integer;
begin
e:=0.00001;
write('a=');
read(a);
repeat
write('b>a  b=');
read(b);
until b>a;
repeat
write('h<',b-a:0:1,' h=');
read(h);
until h<b-a;
writeln('-------------------------------------');
writeln('|  x  |    S    |    y    |    p    |');
writeln('-------------------------------------');
x:=a;
while x<b+h/2 do
 begin
  t:=1;
  s:=t;
  n:=0;
  while abs(t)*(n*n+1)>e do
   begin
    n:=n+1;
    t:=t*(x/2)/n;
    s:=s+t*(n*n+1);
   end;
  y:=(x*x/4+x/2+1)*exp(x/2);
  p:=abs((s-y)/y*100);
  writeln('|',x:5:2,'|',s:9:5,'|',y:9:5,'|',p:9:5,'|');
  x:=x+h
 end;
writeln('-------------------------------------');
end.