int du/(sqrt(1/R-1/u)))