xmin=2;
xmax=11;
nbPts=10;
T = (xmax-xmin);
h = T/(nbPts-1);

for i=1:nbPts

    x(i)=xmin+(T/(nbPts-1))*(i-1);
    y(i) = 2*x(i).^3 ;

end

der = diff(y)/h;

figure(1)
hold on;
plot(x,y);
plot(x(:,1:end-1),der);  %Vue que la matrice der est plus petite que celle de x, on prend l'avant dernière valeur de x
end