Matlab program to calculate drag force

1) % Incorporate a function to plot Drag force vs Velocity

function [DF] = drag_force
V=linspace(0,100,20);

%inputs
%drag coefficient
Cd=0.8;

%Area
A=0.5;

%Density
rho=1.225;

DF=rho*A*V.^2*Cd*0.5;

plot(V,DF);
xlabel(\'velocity(m/s)\');
ylabel(\'Drag force(N)\');
grid on;

2a) Drive link for drag coefficient for various shapes document :

https://drive.google.com/open?id=10hBtb2JlsvjPX9KaGLLf60_Z4ZJUZ8RW

2b) % to plot drag force as a function of drag coeffiecient

function [DF] = drag_force
Dc=linspace(0,1,20);

%inputs
%Velocity
V=25;

%Area
A=0.5;

%Density
rho=1.225;

DF=rho*A*V.^2*Dc*0.5;

plot(Dc,DF)
xlabel(\'Drag coefficient\');
ylabel(\'Drag force(N)\');
grid on;