не работает прога, составить структуированную программу табулирования значений функций Опции

[Fallen59]

задание:
Составить структурированную программу табулирования значений функции S(x) и Y(x) в указанном диапазоне изменения аргумента. Сумму вычислить из k членов ряда. Шаг изменения аргумента принять равным Δx=|xk-xn|/9.
S(x)=1+(2*x/2)+...+((n*n+1)*x^n)/(2^n(*n!))
Y(x)=((x*x/4)+(x/2+1))*e^(x/2)
xn=0.1
xk=1
k=30.



не работает прога, вроде все правидьно но результаты S(x) и Y(x) не совпадают помогите кто нить плиз((

Program z3n20;
uses crt;
const
  n=30;
  xn=0.1;
  xk=1;
var
  dx,x,xi,s,y,f: real;
  i:byte;
begin
 clrscr;
 dx:=abs(xk-xn)/9;
 x:=xn;
 writeln('ЙНННЛННННННЛНННННН»');
 writeln('є x є   s  є   y  є');
 writeln('МНННОННННННОНННННН№');
 while x<=xk do
   begin
   y:=(sqr(x)/4+x/(2+1))*exp(x/2);
   s:=1;
   xi:=x;
   for i:=1 to n do
     begin
   f:=((sqr(i)+1)*x)/((sqr(i-1)+1)*2*i);
   xi:=xi*f;
   s:=s+xi;
   end;
    writeln('є',x:1:1,'є',s:2:4,'є',y:2:4,'є');
     writeln('МНННОННННННОНННННН№');
               x:=x+dx;
              end;
      writeln('ИНННКННННННКННННННј');
 readln;
end.