clear; close all
%
%  Skydiver problem
%
m = 90;     %% kg
g = 9.81;   %% m/sec^2
c = 30;     %% kg/sec
power = 1.1;
v0 = 0;
t = [0 15]; %% 15 sec duration
colors = 'rgbkm';
for power = [1 1.1 1.2 1.5 2]
    etc = [m g c power];
    [t,v] = ode45(@skyfun,t,v0,[],etc);
    plot(t,v,colors(1));
    colors(1) = [];
    hold on
end
grid
legend('p=1','1.1','1.2','1.5','2','Location','NorthWest')
xlabel('Tim, sec'); ylabel('Velo, m/sec');