.png)
Try this:
t = linspace(0, 2, 25);
h = t(2)-t(1);
s = sin(2*pi*t/max(t));
ds = gradient(s,h);
dt = gradient(t,h);
figure
quiver(t,s, dt,ds)
grid
producing:
%20-%202019%2006%2013.png)
This calculates the derivative from the data (here sin(t)). You may be able to recover your derivative directly from your ODE, so you may not need to calculate the gradient of your solved vector. I will leave that for you to determine.
