How to use the script

1. STSDAS

- stsdas> contisub em.fits continuum.fits

2. Try many possible μ s

- find μ_min and μ_max

3. Calculate the skewness function s(μ) : "_tempskew.txt"

4. Plot "_tempskew.txt"

5. Choose the final μ_opt