Matlab program to calculate drag force against a cyclist

clear all
close all
clc
%input
%drag coefficient
c_d=0.50
%area m^2
a=0.2
%density kg/m^3
rho=1.2
%velocity m/s
v=[1:20];
drag_force=rho*a*v.^2*c_d*0.5
plot(v,drag_force)
xlabel(\'velocity\')
ylabel(\'drag_force\')
grid on