%TestLiftDragCoeff

aoa = -400:2:400;
aoa = 12;
Re = 10000:10000:20000;
for i = 1:length(aoa)
    [CL CD] = LiftDragCoeff(Re,aoa(i));
end

figure
plot(aoa,CD,aoa,CL);
xlabel('angle of attack (deg)');
ylabel('Coefficient');
legend('Drag Coeff','Lift Coeff');