Написать программу разложения (x+y) в степени n по формуле бинома Ньютона

program EvgeniiSp;
var
 koef: real;
 k, n, a, b, c: integer; 
 
 function fact(x: integer): integer;
  var
   i, rez: integer;
  begin
   rez:=1;
   for i:=1 to x do rez:=rez*i;
   fact:=rez; 
  end;
  
begin
 k:=0;
 write('Введите значение n: ');
 readln(n);
 while n >= k do
  begin
   // вычисление биноминальных коэффициентов
   a:=fact(n);
   b:=fact(k);
   c:=fact(n-k);
   koef:=a/(b*c);
   write(koef,'*a^',n-k,'*b^',k); // ^ - означает возведение в степень
   // вычисление показателей степеней
   k:=k+1;
   if n >= k then write('+');
  end;
end.

Program nuton;
Uses crt;
Var a,l,j,n,p,k,g,b:integer;
Begin
Clrscr;
Write('Введите степень суммы a и b: ');Readln(n);
Write('(a+b) в степени ',n,' = ');
l:=0; j:=0; a:=1; k:=0; p:=0; b:=n; g:=0;
Repeat
Repeat
write('a');
l:=l+1;
if l<>n then write('*');
Until l=n;
if a>=2 then
Repeat
k:=k+1;
write('*b');
Until k=p;
l:=j; n:=n-1;
write('+');
a:=a+1;
p:=p+1;
k:=0;
Until n=0;
Repeat
Write('b');
g:=g+1;
if g<>b then
Write('*');
Until g=b;
End.