Flow over bicycle

------------------------------% program to calculate drag force on bicycle
clear all
close all
clc

%inputs

c_d=0.8 % drag coefficient
v=[1 20 25 40] % velocity
A=0.1 % frontal area
rho=1.2 % density
drag_force=c_d*rho*A*v.^2*0.5
figure(1)
plot(v,drag_force,\'color\',\'r\')
xlabel(\'Velocity\')
ylabel(\'Drag force\')
grid on

-----------------------------------------% drag coefficient c_d for various shapes
v1=30
c_d=[0.47 0.42 0.50 1.05 0.80 0.82 1.15 0.04 0.09] % for various shapes
drag_force=rho*A*0.5*c_d.*v1^2
plot(c_d,drag_force,\'color\',\'b\')
xlabel(\'c_d\')
ylabel(\'drag force\')
grid on