close all
clear all
clc

vel = [-0.15:0.01:0.15];
rho = 1027;
Cd = 1.5;
A = 0.1018;
v_nom = 0.1;


for i = 1:length(vel)
    quadratic_drag(i) = -0.5*rho*Cd*A*vel(i)*abs(vel(i));
    linear_drag(i) = -0.5*rho*Cd*A*vel(i)*v_nom;
end

figure(1)
plot(vel, quadratic_drag)
hold on
plot(vel, linear_drag)
xlabel("meters/s")
grid on