Flow over Bicycle and calculation of drag force

% this program is to calculate the drag force of air on a set of given values
clear all
close all
clc

%inputs
c_d = 0.8; %drag co-efficient
A = 0.2; %area in m^2
rho = 1.2; %density in kg/m^3
v = [1:100];
v2 = 1; %in m/s, fixed.
c_d2 = [0.47,0.42,0.82,1.15,0.04,0.09];
%{
an array of drag coefficients for different shapes in order:
1. sphere
2. hemisphere
3. long cylender
4. short cylender
5. streamlined body
6. half streamlined body
%}

%calculations
drag_force = 0.5*rho*v.*v*c_d;
drag_force2 = 0.5*rho*v2^2.*c_d2;

%plotting
plot(v,drag_force,\'color\',\'b\')
xlabel(\'velocity\')
ylabel(\'drag force\')
grid on
hold on

figure;
plot (c_d2,drag_force2,\'color\',\'g\')
xlabel(\'co-efficient of drag\')
ylabel(\'drag force\')
grid on
hold off