Flow over bicycle using matlab

clc
clear all
close all

%Inputs
c_d = 0.8

%Area
A = 0.1

%Density
rho = 1.2

% Velocity (m/s)
V = [1,2,3,4,5];

% Drag force calculation
drag_force = 0.5*rho*A*V.^2*c_d

%plot graph for velocity vs drag force
plot(V , drag_force)
xlabel('Velocity')
ylabel('drag force')
grid on

clc
clear all
close all

%Inputs
c_d = [0.47,0.42,0.50,0.80,0.82];

%Area
A = 0.1

%Density
rho = 1.2

% Velocity (m/s)
V = 1

% Drag force calculation
drag_force = 0.5*rho*A*V.^2*c_d

%plot graph for velocity vs drag force
plot(c_d , drag_force)
xlabel('coefficient of drag')
ylabel('drag force')
grid on