function [y] = plugflow_function(R)

y = (log((1+R*(1-0.95))/(R*(1-0.95)))) - ((R+1)/(R*(1+R*(1-0.95))));