Drag force of a cyclist

drag force vs velocity

clear all
close all
clc

cd=0.8
a=0.1
density=1.2
velocity= [1:50];


dragforce=density*a*velocity.*velocity.*cd*0.5;

plot(velocity,dragforce)
xlabel('velocity')
ylabel('dragforce')
grid on

drag force vs drag coefficient

clear all
close all
clc

cd=linspace(0,5,10)
a=0.1
density=1.2
velocity= 1;


dragforce=density*a*velocity*velocity*cd.*0.5;

plot(cd,dragforce)
xlabel('drag co-efficient')
ylabel('dragforce')
grid on