Grafico accellerazione(derivata) in matlab con metodo di runge kutta

[vros90]

Salve a tutti io ho questo problema, devo plottare il grafico dell'accelerazione di un sistema oscillante ad un grado di libertà risoto con runge-kutta. sono riuscito a plottare sia il movimento che la velocità ma non riesco con quello dell'accellerazione chi mi da un manina a sbloccarmi ?

qui sotto vi posto parte del codice

function dydt = vdp2(t,y)
dydt = zeros(2,1);
dydt(1)=y(2);
dydt(2)=(F0*ecc..);

nel main:
plot(t, y(:, 1))
grid
xlabel();
ylabel();
title()
hold on;
plot(t,y(:,2),'green')

dal codice si vedce che stampo velocità e spostamento ma l'accellerazione come si fa ??

grazie a tutti delle risposte !!