Flow over Bicycle coding

% entering value of variables and constants
cd = 0.8; 
a = 1;
rho = 1.2;
v = [1 2 3 4 5];
%calculating drag force
drag_force = 0.5*cd*a*rho*v.^2;
plot(v,drag_force);

% entering value of variables and constants
v = 10; 
a = 1;
rho = 1.2;
cd = [0.1 0.2 0.3 0.4 0.5];
%calculating drag force
drag_force = 0.5.*cd*rho*a*v^2;
%plotting drag force vs drag coefficient
plot(cd,drag_force);