Drag Force over a bicycle and Behaviour of Resistance for diff coefficient

%Program to calculate drag force over a bicycle

clear all
close all
clc

%inputs
%Drag coefficient
c_d = 0.8

%Density
p = 2
%Fontal Area
a = 0.3

%Velocity
v = [1:10];
w = 2

drag = p*a*v.^2*c_d*0.5;

%for same frontal area different drag coefficient for different shapes
% shapes inorder flat plate,prism,sphere,bullet,airofoil
c_dall = [1.28 1.14 0.7 0.295 0.045];

q = p*a*w.^2*c_dall*0.5;

%Plotting
figure(1)
plot(v,drag)
xlabel(\'velocity\')
ylabel(\'drag force\')
grid on

figure(2)
plot(c_dall,q)
xlabel(\'coefficient\')
ylabel(\'drag force\')