close all
clear all

d = 30*.0254;
l = 33*.0254;

nu = 1e-6;

Amplitudes = 1:5;
Periods = 3:20;


for i = 1:length(Amplitudes)
    A = Amplitudes(i);
    leg_str{i} = ['Amp = ' num2str(A) ' m'];
   for j = 1:length(Periods);
       T = Periods(j);
       w = 2*pi/T;
        
       U = A*w/sqrt(2);
       
   ReynoldsNumber(i,j) = U*d/nu;    
       
   end
end
       
figure
plot(Periods,ReynoldsNumber);
xlabel('Wave Period (s)');
ylabel('Reynolds Number');
title(['ESP diameter = '  num2str(d/.0254) ' in']);
legend(leg_str);