Try:
FindMaximum[f[x], {x, 1/3}] (* With starting point 1/3*) NMaximize[{f[x], 0 < x < 1}, x] (*{1.09904*10^6, {x -> 0.263896}}*)
Attachments:
