Решение задачи Proc 45

#include <stdio.h>
#include <math.h>

float power4(float x, float a ,float e){
    float power4_n=1,rez=0;
    int i=0;
 
    while (fabs(power4_n)>e){
        rez+=power4_n;
        ++i;
        power4_n*=x*(a-i+1)/i;
    }
    return rez;
}
 
int main(void)
{
    float x,a,e;
    printf("X:");
    scanf("%f", &x);
    printf("A:");
    scanf("%f", &a);
    int i;
    for(i=1;i<=6;++i){
        printf("E: ");
        scanf("%f", &e);
        printf("(1+x)^a: %f\n",power4(x,a,e));
    }
    return 0;
}

import random
import math

def Power4(x,a,eps):
    if eps <= 0:
        print("Epsilon should be greater than 0")
    y = a*x
    f = 1 + y
    i = 2
    while abs(y) > eps:
        y *= (a-i+1) * x / i
        i += 1
        f += y
        #print(i,":",y,":",f)
    return f

eps = 0.01
a = 10
x = 0.1
for i in range(0, 6):
    #x = -0.5
    #x = 2
    print("( 1 +",x,") ^",a,"=",Power4(x,a,eps),";",(1+x)**a)
    eps /= 10
									

program Proc45;
 
Function Power4(x,a,e:Real):Real;
var
Znak,i:integer;
xi,ni:Real;
TempRez,TempRez2,Rez:Real;
 
begin
xi:=1;
i:=0;
ni:=1;
Znak:=-1;
Rez:=0;
TempRez:=1;
TempRez2:=1;
While (abs(TempRez) > e) do
 begin
 
  Rez:=Rez+TempRez;
  Writeln(i,' ',xi,' ',ni,' ',TempRez,' ',TempRez2,' ',Rez);
  Readln;
 
  inc(i);
  xi:=xi*x;
  ni:=ni*i;
  TempRez2:=TempRez2*(a-i+1);
  TempRez:=TempRez2*xi/ni;
 
 end;
Power4:=Rez;
 
end;
 
var
  X,A,E:real;
  i: integer;
 
begin
   Write('X :');
   Readln(X);
   Write('A :');
   Readln(A);
   for i:=1 to 6 do
    begin
     Write('E :');
     Readln(E);
     Writeln(Power4(X,A,E));
    end;
end.