Fix RMS uncertainty calculation

src/data/timing/foucalc.sl

@@ -191,7 +191,7 @@ define foucalc()
       
       variable sigrms = sqrt( sum(signormpsd[[idxmin:idxmax-1]] * df) );
       variable noirms = sqrt( sum(noinormpsd[[idxmin:idxmax-1]] * df) );
-      variable errrms = sqrt((2.*sigrms^2*noirms^2+noirms^4)/(2.*numseg*length(idx)*sigrms^2));
+      variable errrms = sqrt((4.*sigrms^2*noirms^2+2.*noirms^4)/(numseg*length(idx)*sigrms^2));  % Vaughan+03 Eq. 11
       set_struct_field(rmsquant, "sigrms"+istr, sigrms);
       set_struct_field(rmsquant, "errrms"+istr, errrms);
       set_struct_field(rmsquant, "noirms"+istr, noirms);

src/data/timing/lag_energy_spectrum.sl

@@ -388,7 +388,7 @@ define lag_energy_spectrum(lc_list, dimseg) {
       % copy-paste the calculation here!)
       variable sigrms = sqrt( sum(fouquant.signormpsd2[[fouquant_f.idxmin:fouquant_f.idxmax-1]] * fouquant_f.df) );
       variable noirms = sqrt( sum(fouquant.noinormpsd2[[fouquant_f.idxmin:fouquant_f.idxmax-1]] * fouquant_f.df) );
-      variable errrms = sqrt((4.*sigrms^2*noirms^2+2.*noirms^4)/(numseg*fouquant_f.n_freqs*sigrms^2));
+      variable errrms = sqrt((4.*sigrms^2*noirms^2+2.*noirms^4)/(numseg*fouquant_f.n_freqs*sigrms^2));  % Vaughan+03 Eq. 11
 
       dat.rms[ii,jj] = sigrms;
       dat.rms_err[ii,jj] = errrms;