flow over bicycle both part written seperately

% part 1

%coefficient of drag
cd=0.6;
%velocity vector
vel=1:100;
%density
den=1.2;
%frontal area
fa=1.5;
%drag force
df=0.5*cd*den*fa*vel.^2;
%plot drag force vs velocity
plot(vel,df)
xlabel(\'VELOCITY->\')
ylabel(\'DRAG FORCE->\')

%part 2

%frontal area of different shapes
fa_square=1*1;
fa_rectangle=1*1.2;
fa_circle=pi*1^2/4;
fa_triangle=0.5*1*1;
fa=[fa_square,fa_rectangle,fa_circle,fa_triangle];
%density,vel,coefficient of drag
dens=1.2;
vel=1;
cd=0.8;
%drag force
df=0.5*fa*dens*vel*vel*cd;
bar(df);
xlabel(\'Same Drag Coeff of square,rectangle,circle,triangle\');
ylabel(\'Drag Force\');