Flow over a cyclist

%program to calculate drag force on a bicycle
clear all
close all
clc

%inputs
%drag coefficient
c_d = linspace(0.01,1.4,50);

%velocity m/s
v = linspace(1000,50000,50);

%area m^2
a = 2;

%density kg/m^3
rho = 1.225;

%drag force
drag_force = 0.5*c_d.*v.^2*rho*a

figure(1)
plot(v,drag_force)
xlabel(\'velocity\')
ylabel(\'drag force\')
grid on

figure(2)
plot(drag_force,c_d)
xlabel(\'drag force\')
ylabel(\'drag coefficient\')
grid on