Drag force on cycle

Objective: % calculate the drag forceclear allclose allclc % inputsc_d = 0.8; %drag coefficientA = 0.1; %frontal arearho = 1.2; %densityv = [0:50]; %velocitydrag_force = 0.5*c_d*A*rho*v.^2; %plotfigure(1);plot(v, drag_force); %various drag coefficient problemc_d = linspace(0,1,51);drag_force = 0.5.*c_d.*(A*rho).*v.^2figure(2);plot(c_d,…