%Ricky Palomares
%Dr. Lin
%ARO305-03
%
%HW3-P3 TAFT
clc
clear all
%Note: n is normalized x, that is x/c = n
%Note: Tau is normalized zmax, that is, Tau = zmax/c
tau = 0.04; %max camber normalized to chord

%unit contstants
V = 1; %freestream velcoity

%alpha sweep for CP plot
%part 6 of handout
alpha = (-12:0.5:15);
xCP = (alpha + 4 * tau)./(4*alpha + 8 * tau);
figure(1)
plot(alpha, xCP);
xlabel('Alpha, Degrees'); ylabel('xCP, percent Chord');
title('xCp verus Alpha')
%MAKE TABLE

%Chordwise lifting pressure alpha sweep
%part 9 of handout
alphalift = [-10, (-2*tau)*180/pi, 0, 1, 3, 5]*pi()/180; %aoa to radians
nsweep = linspace(0, 1, 50); %50 steps to get to x/c = 1

figure(4)
hold on
markers = ['-+'; '-.'; '-o';'-x'; '-*'; '->'];
for i  = 1:length(alphalift)
    liftinggam(i, :) = 2*(2 .* V .*(alphalift(i).*(( 1 - nsweep)./(nsweep))...
        .^0.5 + 8 .* tau .* (nsweep .* ( 1 - nsweep)).^0.5));
    plot(nsweep, liftinggam(i, :), markers(i,:))
end

xlabel('X/c'); ylabel('Delta Cp');
title('Chordwise Lifting Pressure Coefficent');
legend('AoA = -10 deg', 'AoA = -2Tau deg', 'AoA = 0 deg', 'AoA = 1 deg',...
    'AoA = 3 deg', 'AoA = 5 deg');
axis([0, 1, -2, 2]);