variation of draG force with velocity of an object and draG coefficient

%flow over a bicycle
clear all;
close all;
clc;
%Input
rho = 0.2;
A=2;%A is frontal area
cd=linspace(0.1,0.9,10); % draG coeff when velocity is const

%velocity when draG coeff is const

v = [1:50];
%Calculations

cd_1=0.3%a random value to visualise variation of draG with velocity

%variation of draG with velocity

draG_force_1 = 0.5*rho*v.^2.*cd_1*A;

%variation of draG with cd

v_1=30;
draG_force_2 = 0.5*rho*v_1^2.*cd*A

%plottinG
ct=1;
 hold on;

   
Graph_1(v,draG_force_1);
Graph_2(cd ,draG_force_2);


%for j=1:size(cd)
%plot(cd(i) ,draG_force , \'color\' , \'b\');
%end

Two different functions Graph_1 and Graph_2 are defined to show variation of draG_force with velocity and cd respectively.

They are defined as follows:

function makeGraph = Graph_2(cd,draG_force_2)
graph=figure(2);
plot(cd,draG_force_2 , \'color\' , \'b\');
grid on;
end

function makeGraph = Graph_1(v,draG_force_1)
graph=figure(1);
plot(v,draG_force_1 , \'color\' , \'r\');
grid on;
end

The final plots are as shown: