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

npts = 101;
m = linspace(120,160,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);

gm = v-vdesired+m;

plot(m,gm,'r-',m,m,'k-')