%% 13.00 hw7
%% Script to make all the plots
%% (After that, I used ginput() to select points on the graphs, as shown
%% in class, then I computed c_g and c_p manually
clear all;
close all

load half_hz.dat
waves = half_hz;
figure(1)
plot(waves(:,1), waves(:,2), 'r');
hold on
plot(waves(:,1), waves(:,3), 'g');
plot(waves(:,1), waves(:,4), 'k');
grid on
zoom on
title('13.00 wave probe data for 0.5 Hertz');
print -depsc f1.epsc


load one_hz.dat
waves = one_hz;
figure(2)
plot(waves(:,1), waves(:,2), 'r');
hold on
plot(waves(:,1), waves(:,3), 'g');
plot(waves(:,1), waves(:,4), 'k');
grid on
zoom on
title('13.00 wave probe data for 1.0 Hertz');
print -depsc f2.epsc


load one_and_a_half_hz.dat
waves = one_and_a_half_hz;
figure(3)
plot(waves(:,1), waves(:,2), 'r');
hold on
plot(waves(:,1), waves(:,3), 'g');
plot(waves(:,1), waves(:,4), 'k');
grid on
zoom on
title('13.00 wave probe data for 1.5 Hertz');
print -depsc f3.epsc

load two_hz.dat
waves = two_hz;
figure(4)
plot(waves(:,1), waves(:,2), 'r');
hold on
plot(waves(:,1), waves(:,3), 'g');
plot(waves(:,1), waves(:,4), 'k');
grid on
zoom on
title('13.00 wave probe data for 2.0 Hertz');
print -depsc f4.epsc