plot (erf(-ln(2/x-1))+1)/2 where 0<x<2