derivative y=sqrt((b*(x-b))/x3)