clear;
close all

v = -.15:.01:.15;
rho = 1026;
Cd = 1.5;
A = pi * .18^2;

Fd = -0.5 * rho * Cd * A * v.*abs(v);

clin = 0.1;
Fdlin = -0.5 * rho * Cd * A * v.*clin;

figure;
plot(v,Fd,'b','linewidth',3)
hold on;
plot(v,Fdlin,'r','linewidth',3)
plot(0.10,(-0.5 * rho * Cd * A * .10 *.10),'go','MarkerFaceColor','g')
axis([-0.15 0.15 -2.5 2.5])
grid on;
legend('Quadratic Drag','Linearized Drag','Linearization Velocity');
xlabel('Velocity (m/s)','fontsize', 15)
ylabel('Drag Force (N)','fontsize', 15)
set(gca,'FontSize',12)
textobj = findobj('type', 'text');
set(textobj, 'fontunits', 'points');
set(textobj, 'fontsize', 15);
