MATLAB: How to add direction line into the phase plane plot

direction fieldphase planetrajectory

Yes, how to add direction line into the phase plane plot? I have the codes below which it work successful but I want to add some direction arrow into my phase plane.
function phase()
IC = [1 1;1 2;1 3;1 4;1 5];
hold on
for ii = 1:length(IC(:,1))
    [~,X] = ode45(@EOM,[-5 5],IC(ii,:));
    u = X(:,1);
    w = X(:,2);
    plot(u,w,'r')
end
xlabel('u')
ylabel('w')
grid
end
function dX = EOM(t, y)
dX = zeros(2,1);
u  = y(1);
w  = y(2);
A  = 1;
B  = 1;
dX = [(w*u^2 - B*u);...
      (A - w - w*u^2)];
end
And the plot below is from the codes above:
And I want to add the direction line, which something like this:
%

Best Answer

Alex, use quiver. In
 quiver(x,y,u,v)
the x, y are the grid points, u and v you get from the dX vector from the differential equations function.
Related Question