
clear;
MAXW=input('Max Vertical velocity?? in m/s ');
dz=8;                                      % delta z
nstep=6000/dz;  
z=1000:dz:7000;
W(1)=0;
for i=2:nstep+1;
W(i)=-1*(MAXW-.5)*(z(i)/1000-1)*(z(i)/1000-7)/9+.5;
if (W(i)<=0) 
W(i)=0.5;
end
dt(i)=dz/W(i);                                % delta t is now an array
end

figure
plot(W,z)
axis([0 MAXW 0 7000])
xlabel('vertical velocity (m/s)')
ylabel('height (m)')

# plot for munich conference
clear;

[W25,z25,dt25]=parabolicw(2,5);
[W55,z55,dt55]=parabolicw(5,5);
[W27,z27,dt27]=parabolicw(2,7);
[W57,z57,dt57]=parabolicw(5,7);



figure
plot(W25,z25,'bx',W55,z55,'g+',W27,z27,'r*',W57,z57,'md')
axis([0 5.5 0 7000])
xlabel('vertical velocity (m/s)')
ylabel('height (m)')
