 1.375 1.625 1.875 2.125 2.375 2.625 2.875 3.125 3.375 3.625 3.875 4.125 4.375 4.750 5.250 5.750
 2.716 2.850 2.757 2.536 2.184 1.752 1.506 1.176 0.973 0.860 0.819 0.782 0.822 0.920 0.922 0.916
 0.339 0.211 0.249 0.497 0.611 0.725 0.686 0.627 0.557 0.546 0.512 0.490 0.495 0.552 0.626 0.675


Fits are for the log10(density)
openw, 3, txt_name
printf, 3, Lset, form='(16f6.3)'
printf, 3, denset, form='(16f6.3)'
printf, 3, stdset, form='(16f6.3)'
close,3

