% plot the data from the Innolight Temperature sweep

% Innolight temperature

InnTemp = [0.60
    .59
    .56
    .52
    .65] + 39;

FreqDiff = [412
    463
    564
    737
    306];


% fit the data
p = polyfit(InnTemp, FreqDiff, 1);
p1 = p(1)/1000;


% plot the data
close all
plot(InnTemp, FreqDiff, 'o', InnTemp, p(1)*InnTemp + p(2), '-')
xlabel('Innolight temperature (deg C)')
ylabel('Frequency Difference (MHz)')
legend('Data', 'Linear fit')
title('Innolight - Lightwave frequency difference versus Innolight temperature. LW temp = 37.2602 C')
text(39.62, 600, ['Slope = ', num2str(p1, '%5.3f'), 'GHz/K'],...
 'HorizontalAlignment','center',... 
 'BackgroundColor',[.7 .9 .7])
grid on
orient landscape
print('-dpdf', 'Innolight_temp_sweep.pdf')