Calculating drag force against a cycle

Objective: clcclear allclose all %Inputsrho = 1.2; a = 0.1; v = [1 : 10]; %drag coefficient for different shapescd1 = {\'aerofoil\', \'sphere\', \'cube\', \'flat plate\',\'bullet\'};cd2 = [0.05 0.47 1.05 1.28 0.295];ct=1for i =1:length(cd2);CD1 = cd1(i); CD2 = cd2(i);b = \'drag force for \'; c = strcat(b,CD1); drag_force = 0.5 *…