load Orbit91.txt
figure

time = Orbit91(:,1);
N65 = Orbit91(:,6);
S65 = Orbit91(:,7);
N15 = Orbit91(:,8);
S15 = Orbit91(:,9);

[f,P,Pred_norm,Pwhite_norm,P95,P95chi2]=myspectrum(time,N65);
subplot(2,2,1);
loglog(f,P,'k'); hold on;
loglog(f,Pred_norm,'r');
loglog(f,Pwhite_norm,'b');
title('Summer N 65')

[f,P,Pred_norm,Pwhite_norm,P95,P95chi2]=myspectrum(time,S65);
subplot(2,2,2);
loglog(f,P,'k'); hold on;
loglog(f,Pred_norm,'r');
loglog(f,Pwhite_norm,'b');
title('Summer S 65');

[f,P,Pred_norm,Pwhite_norm,P95,P95chi2]=myspectrum(time,N15);
subplot(2,2,3);
loglog(f,P,'k'); hold on;
loglog(f,Pred_norm,'r');
loglog(f,Pwhite_norm,'b');
title('Summer N 15');

[f,P,Pred_norm,Pwhite_norm,P95,P95chi2]=myspectrum(time,S15);
subplot(2,2,4);
loglog(f,P,'k'); hold on;
loglog(f,Pred_norm,'r');
loglog(f,Pwhite_norm,'b');
title('Summer S 15');





