function skydiver_graphical
% Function to plot falling speed at t = 4 sec for various
% values of skydiver mass m

npts = 101;
m = linspace(75,175,npts)';
v = ones(npts,1);
g = 9.81;
c = 0.25;
t = 4;

for i = 1:npts
    v(i,1) = skydiver_func1(m(i,1),g,c,t);
end
vdesired = 36*ones(npts,1);

plot(m,v,'r-',m,vdesired,'b-')