function df = dw2df(dw,w)
    %delta wavelength to delta frequency
    df = 299792458 .* dw ./ w.^2;
end