clc;
clear all;
close all;
% intervalo de tiempo
t=0:0.001:10;
% función a derivar
f=sin(t);
% derivada de la función sen(t) con respecto al tiempo.
% es decir se obtiene $\dot{f}(t_k)=\frac{d f(t_k)}{dt_k}$
dfdt = diff(f)./diff(t);
% el vector df tiene dimensión n-1
t1=0:0.001:9.999;% nueva base de tiempo con dimensión n-1 
%compara la derivada aproximada con el método analítico
plot(t1,dfdt,t1,cos(t1))