Week 1 - Flow over Bicycle

%Program to calculate the drag force of the bicyclist
close all
clear all
clc

%Inputs
%Drag coefficient

c_d = linspace(0,5,50);

%Area m^2
A = 10.0;

%Density kg/m^3
rho = 1.105;

%Velocity
v = 1:50;

d_f = rho*A*v.^2.*c_d*0.5;

figure(1)
plot(v, d_f)
xlabel('Velocity')
ylabel('Drag Force')
grid on

figure(2)
plot(c_d, d_f)
xlabel('coefficent of drag force')
ylabel('Drag Force')
grid on