f:sqrt((x-a)*(b-x))/x;

assume(a<b);

assume(a>0);

integrate(f,x,a,b);