First create the function file and name it whatever you would like. I prefer phase-portrait.
% Phase Plot Program
% To use this function, do the following:
% >> phase_portrait(x1, x2, y1, y2, tfinal, 'F', N); for example,
% >> phase_portrait(-5, 5, -5, 5, 10, 'F', 5)function [] =...