This is the code I wrote for graphing the Phase Plot of the differential equation dy/dt = y(y-2)(y-2)(y-4)(y-6), but I am getting errors for this line of the code:plot(y, yPrime, 'Color', 'Red', 'LineWidth', 2, eqPoints, zerosVec, 'k*');
Complete Code:
%%Phase Plot
y = -1 : 0.01 : 7; %define the values of y
yPrime = y.*((y-2).^2).*(y-4).*(y-6); %differential equation
eqPoints = roots([1 -14 68 -136 96 0]); %find the critical points
zerosVec = zeros(length(eqPoints), 1); %create a zeros vector with length equal to the number of roots of the polynomial
plot(y, yPrime, 'Color', 'Red', 'LineWidth', 2, eqPoints, zerosVec, 'k*');xlabel('y'); ylabel('dy/dt');title('Phase Plot');
Best Answer