Flow Over Bicycle Calculating Drag Force

%1. Program to plot Drag Force Vs Velocity with a constant value of Drag Coefficient

clear all
close all
clc

%Inputs
% Drag Coefficient <dimensionless>
dc = 0.9;

% Range of Velocities <m/s>
vel = [1:20];

% Frontal area <m^2>
area = 12.85;

% Density <kg/m^3>
rho = 2.56;

drag_force = 0.5*rho*dc.*vel.*vel*area;

plot( vel, drag_force, \'color\', \'r\', \'linewidth\', 2);
ylabel(\'Drag Force\');
xlabel(\'Velocity\');
grid on;

--------------------------------------------------------------------------------------------------------

%2. Program with Drag Coefficient of different shapes and plot b/w Drag Force & Drag %Coefficent

clear all
close all
clc

%Inputs
% Drag Coefficient for different shapes like bicycle, bus, plane, train respectively %<dimensionless>
dc = [ 0.6 0.8 0.9 1.2 ];

% Constant Velocity <m/s>
vel = 14.56;

% Frontal area <m^2>
area = 12.85;

% Density <kg/m^3>
rho = 2.56;

drag_force = 0.5*rho.*dc*vel^2*area;

plot( dc, drag_force, \'color\', \'b\', \'linewidth\', 1.5);
ylabel(\'Drag Force\');
xlabel(\'Drag Coefficient\');
grid on;