close all
clear all

figure('Position', [55 100 673 844],'PaperOrientation','landscape','PaperPosition',[0 0 11 8.5])

Amp = [0:5:40];  %inches
EnergyLossPerCycle = [0 27 110 243 436 689 1006 1391 1853];  %N-m

H = 2*Amp/39.4;  %Meters
Period = 2*pi*sqrt(9.8*H/0.21)/(0.877*9.8);  %Seconds from PM Spectrum

PowerLoss = EnergyLossPerCycle./Period;
PowerLoss(1) = 0;


subplot(2,1,1)
plot(Amp,EnergyLossPerCycle)
xlabel('stroke amplitude (inches)');
ylabel('Joules lost to heat per cycle');


subplot(2,1,2)
plot(Amp,PowerLoss);
xlabel('stroke amplitude (inches)');
ylabel('Average power lost to heat (Watts)');




print('-dpdf','MaxLoss');
