t = linspace(0, pi);
y = sin(t);
dydt = gradient(y, 1/1000);
ddyddt = gradient(dydt, 1/1000);
clf
hold on
plot(t,y, 'r')
plot(t,dydt, 'g')
plot(t,ddyddt, 'o')