Pascal, значения функции Y(x) и ее разложения в ряд S(x)

[Cors_wp]

Доброго времени суток. Имеется задание:

"Вывести на экран таблицу значений функции Y(x) и ее разложения в ряд S(x) для x изменяющихся от xn до xk с шагом h=(xn-xk)/10. Близость значений S(x) и Y(x) во всем диапазоне значений x указывает на правильность вычисления S(x) и Y(x)."



Xn,Xk,S(x),n,Y(x).

Помогите исправить нижеприведенный код для моего условия. В математике совсем не "бум-бум"...

const n=14;
      k=10;
      a=0.1;
      b=1;
 
var  s,t,x,h,y:real;
     i:integer;
begin
h:=(b-a)/k;
writeln('Табулирование функции Y=(x^2/4+x/2+1)*e^x/2, представленной ее разложением');
writeln(' в ряд  на интервале [0,1;1] при n=14 и с шагом ',h:0:2);
writeln('---------------------------');
writeln('|  x  |   S(x)  |   Y(x)  |');
writeln('---------------------------');
x:=a;
while x<=b+h/2 do
 begin
  t:=1;
  s:=1;
  for i:=1 to n do
   begin
    t:=t*(x/2)/i;
    s:=s+t*(i*i+1);
   end;
  y:=(x*x/4+x/2+1)*exp(x/2);
  writeln('|',x:5:2,'|',s:8:5,' |',y:8:5,' |');
  x:=x+h;
 end;
writeln('---------------------------');
end.

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

[Cors_wp]

тему можно закрывать