Flow over bicycle

PART1: %Program to claculate the drag force clear allclose allclc % Inputs% Drag coefficientc_d=0.8 % Area m^2A=1 % Density kg/m^3rho=1.2 % velocity m/sv= [1:50]; drag_force=0.5*v.^2*rho*c_d*A; plot(v,drag_force)xlabel('Velocity')ylabel('Drag Force') Part 2: % drag force for various shapes and their drag coefficient clear allclose allclc % Inputs% Area m^2A=1 % Density kg/m^3rho=1.2 % velocity m/sv= [1:50] %streamline bodyc1=0.04 %half spherec2=0.43 %spherec3=0.6 %conec4=0.5 %long cylinderc5=0.82 %cubec6=1.05 drag_force=0.5*v.^2*rho*c3*A; plot(c3,drag_force)xlabel('cd')ylabel('Drag Force')grid on