x = 0:0.001:1.01;
h = -0.01:0.001:0.01;
x0 = 1;
pochodna = [];

for i = 1:length(h)
pochodna = [pochodna, (exp(x0+h(i)) - exp(x0))/h(i);];
end

t = x0+h(1):0.001:x0+h(length(h));
plot(t,pochodna,'g');
hold on;
plot(x,exp(x),'b');
axis([0.95 1.05 2.7 2.74]);