// Puissance (1+x) Version 1

temps=0;
tic()

function z=factorielle(n)
  z=1;
  for i=1:n
    z=z*i;
  end
endfunction

function alfa=coeffDL(n)
   for i=1:1:n
    alfa=alfa*(alfa-i);
  end
endfunction

function y=elementSerie(n,x,alfa)
  y=alfa*(x^n)/factorielle(n);
endfunction

y=elementSerie(2,0.1,1.2)

function S=sommeSerie(x,alfa,eps)
  S=1;
  n=1;
  en=elementSerie(n,x,alfa);//initialisations
  while(abs(en)>eps) // tant que l'élément courent est plus grand que l'erreur admise
    S=S+en;
    n=n+1;
    en=elementSerie(n,x,alfa) ;
  end
endfunction
 
//tests

x=0.2;
alfa=1.3;
eps=0.0001;
S=sommeSerie(x,alfa,eps)
erreur=S-((1+x)^(alfa))
temps=toc() 