Flow over Bicycle

Incorporate a function to plot Drag force vs Velocity.(25 points)

%PROGRAM TO DETERMINE DRAG FORCE ON CYCLIST
clear all
close all
clc
%INPUTS
c_d=0.8;
A=0.1;
rho=1.2;
v=linspace(1,20,10); %linear values b/w 1 to 100
drag_force=0.5*rho*A.*v.^2*c_d
df(v,drag_force); %call function to plot v vs df

%function for plotting drag force Vs velocity
function[]=df(v,drag_force)
plot(v,drag_force,\'linewidth\',3,\'color\',\'g\')
xlabel(\'VELOCITY\')
ylabel(\'DRAGFORCE\')
title(\'PROGRAM TO DETERMINE DRAG FORCE ON CYCLIST\')
grid on
legend(\'Velocity VS Dragforce\')
end

Document the drag coefficient for various shapes and try to plot drag force as a function of drag coefficient. (25 points)

%PROGRAM TO DETERMINE DRAG FORCE ON CYCLIST
clear all
close all
clc
%INPUTS
%C_d for shapes(flatplate,prism,bullet,sphere,airfoil)respectively
c_d=[1.28,1.14,0.295,0.5,0.045];
A=0.1;
rho=1.2;
v=10; %for velocity 10
drag_force=0.5*c_d.*rho*A*v^2
plot(c_d,drag_force,\'linewidth\',3,\'color\',\'g\')
xlabel(\'Drag Coefficient\')
ylabel(\'DRAGFORCE\')
title(\'PROGRAM TO DETERMINE DRAG FORCE FOR DIFFERENT SHAPES\')
grid on
legend(\'Drag Coefficient VS Dragforce\')