// Test de binom et binomrec

Timp=zeros(10,10)
Trec=zeros(10,10)

for i=1:10
    for j=1:i
        Timp(i,j)=binom(j-1,i-1)
        Trec(i,j)=binomrec(j-1,i-1)
    end
end

disp(Timp)
disp(Trec)

// Comparaison de la vitesse d'execution de binom et binomrec
timer();binom(10,20);disp(timer())
timer();binomrec(10,20);disp(timer())
// J'obtiens 0.000098 pour binom et 1.596935 pour binomrec

T=[0:0.001:1]

scf(0)
clf()
plot2d(T',[f(T)',bernstein(5,f,T)',bernstein(10,f,T)',bernstein(50,f,T)'],[3:6])


// Pour tester la fonction bernsteinOP on peut par exemple faire disp([bernsteinOP(5,f,T);bernstein(5,f,T)])

// Comparaison de la vitesse d'execution de bernstein et bernsteinOP
timer();bernstein(50,f,T);disp(timer())
timer();bernsteinOP(50,f,T);disp(timer())

// J'obtiens 1.461051 pour bernstein et 0.530149 pour bernsteinOP ce qui montre que bernsteinOP est plus performant
