flow over bicycle both part written seperately

Objective: % part 1 %coefficient of dragcd=0.6;%velocity vectorvel=1:100;%densityden=1.2;%frontal areafa=1.5;%drag forcedf=0.5*cd*den*fa*vel.^2;%plot drag force vs velocityplot(vel,df)xlabel(\'VELOCITY->\')ylabel(\'DRAG FORCE->\') %part 2 %frontal area of different shapesfa_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…