Power vs. velocity for a propeller

%--------------------------------------------------------------------------
%   Copyright (c) 2005 Princeton Satellite Systems, Inc.
%   All rights reserved.
%--------------------------------------------------------------------------


v   = linspace(0,20);
rho = 1.225;
a   = pi*4;
t   = 5000;
p   = zeros(1,length(v));
w   = zeros(1,length(v));
for k = 1:length(v)
  [p(k), w(k)] = PropellerPower( rho, a, t, v(k) );
end

yL = {'Power (W)' 'Induced Velocity (m/s)'};
Plot2D( v, [p;w], 'Velocity (m/s)', yL, 'Propeller');

%--------------------------------------