% f = e^x derivate_ex = @(i) 1 + 0*i ; % f = 1/(1-x) derivate_geom = @(i) factorial(i) ; x = linspace( -1 , .8 , 1000 ) ; plot( x , 1 ./ (1-x) , 'b' , 'LineWidth' , 1 ) ; hold on ; for i_max = 1 : 2 : 40 y = serie_MacLaurin( x , derivate_geom , i_max ) ; plot( x , y ) ; pause( 1 ) ; end