Flow over Bicycle and variation of Drag co-efficient with Drag forces

%program to calculate the drag force over a bicycle
clc
clear all
close all

%density
rho=1.225;

%velocity
v=[1:60];

%Area for a cyclist(Hoods) (sourced from internet)
A=0.32;

%Drag co efficient for a cyclist(Hoods)(sourced from internet)
Cd=0.88;

Drag_force=0.5*rho*A*v.^2*Cd;
plot(Drag_force,v)
xlabel(\'Velocity\')
ylabel(\'Drag Force\')
grid on;

%program to document various Cd values and plot against drag force
clc
close all
clear all

%density
rho=1.225;

%velocity
v=20;

%Area
A=0.02;

%Drag co efficients
Cd=[1.28 1.14 0.295 0.07 0.045]
%Cd of flat plate=1.28
%Cd of Prism=1.14
%Cd of bullet=0.295
%Cd of sphere=0.07
%Cd of Airfoil=0.045

%Drag force
drag_force=0.5*Cd.*v^2*A*rho

plot(drag_force,Cd)
xlabel(\'dragforce\')
ylabel(\'Drag coefficient\')

grid on;