flow over bicycle

%Calculating the drag force against a cyclist

clear all
close all
clc

%Inputs
%Drag coefficient
c_d = 0.7
%Front area 
A=0.2;
%Velocity
v=[1:50]

%Drag force
D_F= 0.5*rho*A*v.^2*c_d

%plotting
plot(v,D_F)
xlabel('velocity')
ylabel('Drag Force')



%Calculating the drag force for different shapes

clear all
close all
clc

%Inputs
%Drag coefficient
c_d = [0:1];
%Front area 
A = 0.3; 
%Density
rho = 1.5;
%Velocity 
v = 20; 

%Calculating the drag force
D_F = 0.5*c_d.*A*rho*v^2;

%Plotting
plot(c_d, D_F)

xlabel('Drag coefficient')
ylabel('Drag Force')