Hi Jitse Niesen, I know what you're saying, I was thinking along the same lines.
It's easy to show that the unconstrained function has no stationary points, but when you add the constraint (a+b+c=3) then it’s a lot more messy.
If you substitute c=3-a-b to make it a function of two variables you get at least two stationary points. One is exactly at (a,b)=(1,1) and there's another somewhere near (a,b)=(0.76,0.74). The first one (1,1) turns out to be a local minima and I'm pretty sure that the second one is a saddle point.
They were pretty messy derivatives and I used maple to help me with those results. The point I was making before was just that if there was some other easy argument (say some sort of symmetry argument or whatever) to show that the maximum occurred when one of the variables was zero then the whole thing could be reduced pretty easily to a one dimensional problem.
BTW. With c=0 and b=3-a then resulting 1d problem gives the following cubic:
[tex]5\,a^3 -13\,a^2 + 15\,a - 9 = 0[/tex].
Numerically this has the (approx) solution a=1.369, (and b=1.631, c=0) which gives f_max approx 12.765.