flow over bicycle

clc
clear all
close all
%assuming all values
Velocity=[5 10 15 20];
dragCoef=0.5;
Area=2;
airDensity=3;
DragForce=(airDensity*dragCoef*Area*Velocity.*Velocity)/2
hold on
plot(DragForce,Velocity,\'color\',\'r\')
xlabel(\'Drag force\')
dragCoef=[0.5 0.5 0.5 0.5]
plot(DragForce,dragCoef,\'color\',\'b\')