%ME 482: Aerodynamics
%Student: Ian Chase
%Collaborators: Rachel Cao, Ariana Munro
%Date: 10/23/2014
clear all; close all; clc;

% Calculates d, Vt, Cp & Lambda for 4,8,16 & 32 paneled circles
[d4,Vt4,Cp4,L4] = AeroHW8(4);
[d8,Vt8,Cp8,L8] = AeroHW8(8);
[d16,Vt16,Cp16,L16] = AeroHW8(16);
[d32,Vt32,Cp32,L32] = AeroHW8(32);

%% Output
% displays the lambda values for  4,8,16 & 32 paneled circles
disp('Lambda of 4 panel circle')
disp(L4);
disp('Lambda of 8 panel circle')
disp(L8);
disp('lambda of 16 panel circle')
disp(L16);
disp('lambda of 32 panel circle')
disp(L32);

% Plots Tangent Velocity against degrees
subplot(2,1,1);
plot(d4*(180/pi),Vt4);
hold all
plot(d8*(180/pi),Vt8);
hold all
plot(d16*(180/pi),Vt16);
hold all
plot(d32*(180/pi),Vt32);
hold all

% Plots Coefficient of pressure against degrees
subplot(2,1,2);
plot(d4*(180/pi),Cp4);
hold all
plot(d8*(180/pi),Cp8);
hold all
plot(d16*(180/pi),Cp16);
hold all
plot(d32*(180/pi),Cp32);
hold all