The entropy function has its maximum at x=1/2
Yet when one uses solve:
x=var('x')
assume(x>0)
solve(derivative(-x*log(x) -(1-x)*log(1-x),x)==0,x)
sage returns something that is correct, but curiously non-specific
[log(x) == log(-x + 1)]
It almost seems like sage can't figure out that x=1/2 solves the above even after the assumption on x.
