Tính tổng dãy

Cho x là số thực nhập từ bàn phím. Tính
       Mexp(x) = (x^o / 0!) + (X^1 / 1!) + (X^2 / 2!) + … + (x^n / n!) + …
   Quá trình tính toán kết thúc khi số hạng cộng vào đã nhỏ hơn một số dương épsilon cho trước.
                                                                 Giải
Uses crt;
Var
       Mc,Ms,Me,x: Real;
Function Mexp(x:real):Real;
       Const e=0.00001;
       Var
           Sh,S:Real; n:Longint;
     Begin
           S:=1; Sh:=1; n:=0;
           While Abs(Sh)>=e do
               Begin
                  Sh:=(Sh*x)/(n+1); S:=S+Sh; n:=n+1;
              End;
          Mexp:=S;
     End;
BEGIN
       ClrScr;
       Write('Nhap x = '); Readln(x);
       Writeln('=> Mexp(',x:6:2,') = ',Me:10:3,' so sanh voi chuan = ',exp(x):10:3);
       Readln;
END.