function [y] = plugflow_function(R)

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

