//example 7.6
//page 272
clc; funcprot(0);
//initialisation of variable
A=50;//area
Cd=0.3;
U=11.11;//speed
rho=1.177;
Fd=1/2*Cd*rho*U^2*A;
disp(Fd,"drag force (N)");
clear