Calculating the Drag force over a bicycle

%calculating drag force over a bicycle.

clear all

close all

clc

%CASE-1

%Inputs

cd=0.8;

a=0.2;

rho=1.2;

v=[1:30];

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

figure(1);

plot(v,drag_force);

xlabel(\'v\');

ylabel(\'drag_force\')

%CASE-2

%Velocity is constant as,

V=25;

a=0.2;

rho=1.2;

%Different values of C_d are,

c_d=[0.47,0.42,0.50,1.05,0.80,1.15];

Drag_force=rho*a*V^2*c_d.*0.5;

figure(2);

plot(c_d,Drag_force);

xlabel(\'c_d\')

ylabel(\'Drag_force\')

end