calculate drag force for flow over bicycle

In this program for flow over bicycle, I take five different values of Drag Coefficient for different shape and different velocity.

I use two \'for\' loop to plot two different graph (1)Drag Force vs Velocity and (2)Drag Force vs Drag Coefficient.

In the the first graph it show that for particular Drag Coefficient and change in Velocity how Drag Force is change.

In the second graph it show that for particular Velocity and change in Drag Coefficient how Drag Force is change.

%A programme to calculate drag for flow over bicycle.
close all
clear all
clc

%inputs
%drag coefficient for various shap sphere, half sphere, cone, long cylinder, stremlined body as below.
c_d = [0.47 0.42 0.50 0.82 0.04];

%resistance or front area m^2
A = 0.3; 

%density of fluid which resist force (air kg/m^3)
rho = 1.225;

%velocity of bicyle m/s
V = 1:15;

%equation for drag force
ct = 1;
for i = 1:length(c_d)
    
  drag_force = c_d(i)*A*V.^2*rho*0.5;
  axis ([-0.5 16 -1 36])
  plot (V, drag_force)
  xlabel(\'velocity\')
  ylabel(\'drag force\')
  
  grid on
  hold on
 ct = ct+1;

end
figure(ct)
for j = 1:length(V)
    r = V(j);
    drag_force = A*r^2.*c_d*rho*0.5;
    
    plot(c_d, drag_force)
    axis ([-0.5 1.5 -1 36])
    xlabel(\'drag coefficient\')
  ylabel(\'drag force\')
    grid on
    hold on
end