clc clear all disp('name') disp('section') disp('date') disp('This program finds the lift of an airfoil') % p=density (kg/m^3) p =1.247; % a=area (m^2) a =2.4510; % cl=lift coef (unitless) cl =0.8100; % v=velocity (m/sec) v =input('enter the velocity (m/s): '); % lift formula is below lift=(1/2)*p*a*cl*(v^2); fprintf('velocity is %g m/s\n',v) fprintf('the lift force is %g N\n',lift)