function lund_sig,n,wang=wang

; Purpose: return significant level (95%) for the Lund changepoint
; Category: statistics
; See-also: lund; lund_plot
; Inputs: n - length of series

if (keyword_set(wang)) then begin
  ns=[10    , 20    , 30    , 40    , 50    , 60    , 70    , 80    , 90    , 100   , 150   , 200   , 250   , 300   , 350   , 400   , 450   , 500   , 550   , 600   , 800   , 1000  , 2500]
  fs=[15.559, 11.948, 11.396, 11.148, 11.068, 11.072, 11.059, 11.072, 11.080, 11.085, 11.127, 11.208, 11.310, 11.406, 11.392, 11.416, 11.474, 11.537, 11.552, 11.549, 11.673, 11.749, 12.064]
endif else begin
  ns=[10   , 25  , 50  , 75  , 100 , 200 , 300 , 400 , 500 , 750 , 1000, 2500, 5000]
  fs=[11.56, 7.37, 6.92, 6.88, 6.91, 7.01, 7.11, 7.18, 7.14, 7.37, 7.42, 7.65, 7.85]
endelse

if (n lt ns(0) or n gt max(ns)) then message,/cont,'Warning: outside range for n'

return,interpol(fs,ns,n)

end