program to plot velocity vs dragforce and drag coefficiant vs dragforce

% A program to calculate the dragforce on a cyclist with respect to various speed

clear all

close all

clc

% input

%drag coefficiant

c_d =.8

% area = m^2

a=.1

% density = kg/m^3

rho =1.2

% velocity

v = [1:50];

drag_force = rho*a*v^2*c_d*.5;

plot (v, drag_force)

xlabel(\'velocity\');

ylabel(\'drag_force\');

% Program to plot drag coefficiant to dragforce for various shapes

clear all

close all

clc

% input

% drag coefficiant of flat plate

c_d = 1.28

% drag coefficiant of sphere

c_d = 1.14

% drag coefficiant of bullet

c_d = .295

% drag coefficiant of airfoil

c_d = .045

% frontal area = m^2 ( same for all shape )

a= .3

% density = kg/m^3

rho = 1.2

% velocity = m/s

v = 40

c_d = [ 1.28 1.14 .295 .5 .045]

drag_force = rho*a*v^2*c_d*.5

plot (c_d, drag_force)

xlabel (\'drag coefficiant\');

ylabel(\'drag_force\');