Computes a preliminary design of a two stage to orbit vehicle.

It computes the delta-V required given an initial velocity (Mach 6 at 30 km) ------------------------------------------------------------------------- See also LambertTOF, AtmDens2, Plot2D, Mag -------------------------------------------------------------------------

%--------------------------------------------------------------------------
%   Copyright 1995-1998 Princeton Satellite Systems, Inc.
%   All rights reserved.
%--------------------------------------------------------------------------

altitude = 30;
g        = AtmDens2( altitude );

v        = 6*340.3*.8671;

r1       = [6378.165;0;0] + [30;0;0];
r2       = [-6758;0;0];

dT       = linspace(.5,10)*3600;
v1       = zeros(3,length(dT));

for k = 1:length(dT)
  vT  = LambertTOF(r1,r2,dT(k),'direct',3.98600436e5,1e-8);
  v1(:,k) = vT(:,1);
end

vMag = Mag(v1) - Mag(v)/1000;

Plot2D( dT, vMag, 'Time (hr)', 'dV1 (km/s)', 'Second Stage' )


%--------------------------------------
% PSS internal file version information
%--------------------------------------