optimal_p_out function

Table and plot the SONE values