drag force on a cycle

Objective: clcclear allclose all %calculating drag force on a cycle %inputs%coefficient of drag frocec_d=0.8; %frontal areaA=0.1; %densityrho=1.2; %velocityv=[1:50]; drag_coeff= 0.5*rho*v.^2*A*c_d plot(v,drag_coeff)