Langevin function