drag force vs velocity

close all 
clear all
clc

%calculation of drag force using f=0.5 * d*a*v.^2*c_d

%drag coefficient
c_d=0.6;

%frontal area in m^2
a=0.6;

%velocity in m/s
velocity=[1:100];


%density in kg/m^3
d=2.3;


%drag force
drag_force=0.5*d*a*velocity.^2*c_d

%plot
plot(velocity,drag_force,\'linewidth\',2,\'color\',\'b\')
xlabel(\'velocity\')
ylabel(\'drag_force\')
grid on