// Sinus Version 1
tic()
function z=factorielle(n)
  z=1;
  for i=1:n
    z=z*i;
  end
endfunction

function y=elementSerie(n,x)
  m=modulo(n,2)
  if(m==0)
      y=0;
  else
      y=((x^n)*(-1)^((n-1)/2))/factorielle(n);
  end 
endfunction

y=elementSerie(1,1)

function S=sommeSerie(x,eps)
  // x est le point pour lequel on veut approcher la fonc
  // eps est la tolérance
  S=0;
  n=1;
  en=elementSerie(n,x);//initialisations
  while(abs(en)>eps) // tant que l'élément courent est plus grand que l'erreur admise
    S=S+en;
    n=n+2;
    en=elementSerie(n,x) ;
  end
endfunction
 
//tests

x=0.21466;
eps=0.0001;
S=sommeSerie(x,eps)
erreur=S-sin(x)
temps=toc() 