f(sin) в pascal

[A_Orlov]

Пожалуйста, помогите решить задачу на Pascal.
Дано вещественное число X и целое число N (> 0). Найти значение выражения:X – X3/(3!) + X5/(5!) – … + (–1)N•X2•N+1/((2•N+1)!) (N! = 1•2•…•N). Полученное число является приближенным значением функции sin в точке X.

[Филантроп]

uses crt;
var x,n,s,t,ch,zn:real;
    i:integer;
begin
clrscr;
write('x=');
readln(x);
write('n=');// типа 0,0001
readln(n);
s:=x;//начальная сумма
ch:=x;//числитель первой дроби
zn:=1;//знаменатель первой дроби
t:=ch/zn;//значение первого члена ряда
i:=1;//номер
while i<=n do
  begin
   ch:=ch*(-1)*x*x;//домножаем числитель на х в квадрате и меняем знак
   zn:=zn*2*i*(2*i+1);//домножаем знаменатель на 2 следующих номера
   t:=ch/zn;//вычисляем очередной член ряда
   s:=s+t;//считаем очередную сумму
   inc(i);//увеличиваем номер
 end;
writeln('Y=        ',s:0:6);//сумма
writeln('sin(',x:0:2,')=',sin(x):0:6);//аналитическое значение
readln
end.