Hi everyone,

I have below integral. The result must be real. I would like to simplify and get to real part but I couldnt.

> fs := (1-1/sqrt((Uo+U*sin(x))^2+a^2))*(Uo+U*sin(x)); fss := `assuming`([2*(int(fs*sin(x), x = 0 .. Pi))/Pi], [U::real, U::positive, a::real, a::positive, Uo::real, Uo::positive]);

I used evalc but it didn't give good results as well.

thanks