//(1+x)PuissanceAlfa Version 1

function A=combinaison(alfa,n)
  A=alfa
  for i=2:1:n
    A=A*(alfa-i+1)/i
  end
endfunction

function en=elementSerie(x,n)
  A=combinaison(alfa,n)
  en=A*x^n
endfunction

function S=sommeSerie(x,alfa,eps)
  S=1+elementSerie(x,1);
  n=2;
  en=elementSerie(x,n)
  while abs(en)>eps
    S=S+en
    n=n+1
    en=elementSerie(x,n)
  end
endfunction

x=0.34;
alfa=(1/3)

s=sommeSerie(x,alfa,eps)
erreur=s-1.34^(1/3)
1.34^0.33
1.34^(1/3)