Add routine to plot theoretical brettel's threshold