I need help writing this code in Matlab, 
$$\frac{dx}{dt} = x-2y+1$$ $$\frac{dy}{dt} = 3y-3x$$
Using MATLAB, construct a phase portrait using the quiver command. Then, approximate the solutions for the initial conditions Y(0) = (-1,3), Y(0) = (-1,-2), and Y(0) = (1,1). Your solutions should be plotted from t = 0 to t = 10 on the same axes as your phase portrait. The graph should include a legend so that the solutions are clearly visible. Make comments about both the phase portrait and your solutions.
Also, if there is any version of pplane8 that works with MATLAB 2017B that would also be very appreciated.