PROGRAM TO CALCULATE DRAG FORCE

% A PROGRAM TO CALCULATE DRAG FORCE ON BICYCLE

clear all

close all

clc

%Inputs

%Drag coeffecient for different geometries

c_d=1inspace(0.04,1.15,20);

%Area m/s^2

a=0.125

%Density of air kg/m^3

rho=1.225

%Velocity m/s assuming speed varying from 1-60km/hr

v=1inspace(0.2,16.67,20);

drag_force=rho*a*v.^2.*c_d*0.5;

%plotting drag force against velocity

plot(v,drag_force)

xlabel(Velocity)

ylabel(Drag force)

close all

%plotting drag force against drag coeffecient

plot(c_d,drag_force)

xlabel(Drag cooeffecient)

ylabel(Drag force)