close all
clear all

m = 1000;  %kg
a = .76/2;  %radius (m)
b = 1.7/2;  %m
g = 9.8;  %m/s^2
rho = 1025; %kg/m^3

Vol = 4*pi*a*b*b/3; %m^3

W = rho*g*Vol-m*g;
u = rho*Vol;  %Added mass (from Newman Pg 147)


t = 0:.05:2.0;
n = 1;
for V=-1:-1:-3
y(n,:) = (.5*W/(m+u))*t.^2+V.*t;
vel(n,:) = (W/(m+u)).*t+V;
n = n+1;
end

plot(t,y)

legend('Vel = -1 m/s','Vel = -2 m/s','Vel = -3 m/s','Location','southwest');


foo = axis;

foo(4) = 0;
axis(foo);

xlabel('time (seconds after anchor impact)');
ylabel('meters');
foo = sprintf('Landing Transient of Termination after Anchor Touchdown \n Inertial Considerations only, Drag neglected \n m = %d kg, %2.1f m X %2.1f m ellipsioid',m,2*b,2*a);
title(foo);


print -depsc2 'LandingTransient'