Drag force calculation by using different velocity values and drag force coefficient values

% Drag force against bicycle
clear all
close all
clc

% Inputs
% density Kg/m^3
rho = 1.225
% area m^2
a = 0.32
% velocity m/s
v = [1:50]
% drag coefficient
cd = 0.88

% Drag force
F = 0.5*rho*a*v.*v.*cd

plot(v,F,\'color\',\'b\')
xlabel(\'Velocity\',\'color\',\'r\')
ylabel(\'Dragforce\',\'color\',\'r\')
grid on

% Drag Force against a cyclist
clear all
close all
clc

% Inputs
% Sphere cd
scd = 0.47
% half sphere cd
hscd = 0.42
% cone cd
ccd = 0.50
% cube cs
cubecd= 1.05
% angled cube cd
angledcubecd = 0.80
% long cylinder cd
lcd = 0.82
% short cylindr cd
scd = 1.15
% streamlined body
streamcd = 0.04
% streamlined half body cd
shcd = 0.09
% Density of air Kg/m^3
rho = 1.225 
% Area m^2
a = 0.32
% Velocity m/s
v = 20
% Coefficient of drag
Cd = [0.47 0.42 0.50 1.05 0.80 0.82 1.15 0.04 0.09]

% Drag Force
F = 0.5*rho*a*v.*v.*Cd;

plot(Cd,F, \'color\',\'r\')
xlabel(\'Drag force coefficient\')
ylabel(\'Dragforce\')
grid on