plot x = 2* (exp(x*y)-1.0)/(exp(y)-1.0)