Drag Force on Bicycle

%Drag Force
clear all
close all
clc

%inputs

dragF= 0; % N
ForceDrag =0; % N
a= 0.8; %m^2
cd= 0.5
c_d= [0.47 0.04 1.05 0.80 1.15]; % sphere = 0.47 ,streamlined_body = 0.04, cube= 1.05, angled cube= 0.8, short cylinder = 1.15
rho = 1.125; %kg/(m^3)
v1= 2; %m/s
%to input range of values
v = [1:5]; %m/s

%formula

dragF = a*cd*v.*v.*rho*0.5
plot(v, dragF,\'color\',\'b\'); % plotting
axis ([0 8 0 10]) % setting limits for axis
xlabel(\'velocity\')
ylabel(\'Drag Force\')
title({\'Drag Force vs Velocity\'}); % title of the graph
grid on

figure;
ForceDrag = a*c_d.*v1*v1*rho*0.5
plot(c_d,ForceDrag, \'color\',\'r\');
axis ([0 1.5 0 5])
xlabel(\' c_d\')
ylabel(\'Drag Force\')
title({\'Drag Force vs drag coefficient(C_d)\'});
grid on