<s>[INST] Prove $\sqrt{3a + b^3} + \sqrt{3b + c^3} + \sqrt{3c + a^3} \ge 6$ 
If $a,b,c$ are non-negative numbers and $a+b+c=3$, prove that:
$$\sqrt{3a + b^3} + \sqrt{3b + c^3} + \sqrt{3c + a^3} \ge 6.$$


Here's what I've tried:
Using Cauchy-Schawrz I proved that:
$$(3a + b^3)(3 + 1) \ge (3\sqrt{a} + \sqrt{b^3})^2$$
$$\sqrt{(3a + b^3)(4)} \ge 3\sqrt{a} + \sqrt{b^3}$$
$$\sqrt{(3a + b^3)} \ge \frac{3\sqrt{a} + \sqrt{b^3}}{2}$$
Also I get:
$$\sqrt{(3b + c^3)} \ge \frac{3\sqrt{b} + \sqrt{c^3}}{2}$$
$$\sqrt{(3c + a^3)} \ge \frac{3\sqrt{c} + \sqrt{a^3}}{2}$$
If I add add 3 inequalities I get:
$$\sqrt{(3a + b^3)} + \sqrt{(3b + c^3)} + \sqrt{(3c + a^3)} \ge \frac{3\sqrt{a} + \sqrt{b^3}}{2} + \frac{3\sqrt{b} + \sqrt{c^3}}{2} + \frac{3\sqrt{c} + \sqrt{a^3}}{2}$$
Now i need to prove that:
$$\frac{3\sqrt{a} + \sqrt{a^3}}{2} + \frac{3\sqrt{b} + \sqrt{b^3}}{2} + \frac{3\sqrt{c} + \sqrt{c^3}}{2} \ge 6 = 2(a+b+c)$$
It's enough now to prove that:
$$\frac{3\sqrt{a} + \sqrt{a^3}}{2} \ge b+c = 3-a$$
$$\frac{3\sqrt{b} + \sqrt{b^3}}{2} \ge a+c = 3-b$$
$$\frac{3\sqrt{c} + \sqrt{c^3}}{2} \ge b+a = 3-c$$
All three inequalities are of the form:
$$\frac{3\sqrt{x} + \sqrt{x^3}}{2} \ge 3-x$$
$$3\sqrt{x} + \sqrt{x^3} \ge 6-2x$$
$$(3\sqrt{x} + \sqrt{x^3})^2 \ge (6-2x)^2$$
$$9x + x^3 + 6x^2 \ge 36 - 24x + 4x^2$$
$$x^3 + 2x^2 + 33x - 36 \ge 0$$
$$(x-1)(x^2 + 3x + 33) \ge 0$$
Case 1:
$$(x-1) \ge 0 \ \ \ \ \text{ for any }\ x \geq 1$$
$$(x^2 + 3x + 33) \ge 0 \ \ \ \ \text{ for any x in R} $$
Case 2:
$$0 \ge (x-1) \ \ \ \ \text{ for any }\ 1 \geq x$$
$$0 \ge (x^2 + 3x + 33) \ \ \ \ \text{there are no solutions in R} $$
This proves that for $$x \geq 1$$ $$(x-1)(x^2 + 3x + 33) \ge 0$$ is true and so it is
$$\sqrt{3a + b^3} + \sqrt{3b + c^3} + \sqrt{3c + a^3} \ge 6$$, but a, b, c can be every non-negative number. I proved it's true for $$a,b,c \geq 1$$, but i can't for $$a,b,c \geq 0$$
 [/INST] I'm not sure if "proof with Mathematica" counts, but here it is:
In order to find the local extrema of function $f(a,b,c) = \sqrt{3a+b^3}+\sqrt{3b+c^3}+\sqrt{3c+a^3}$ under the constraint $a+b+c=3$, we differentiate the corresponding lagrangian and arrive at a system of equations
$$\left\{\begin{aligned}& (a^2c^2-a^2+1)^2(3a+b^3)=(b^2a^2-b^2+1)^2(3b+c^3)=(c^2b^2-c^2+1)^2(3c+a^3) \\& a+b+c=3 \\& c \ge a,b \ge 0 \end{aligned}\right.$$
Obviously, $a=b=c=1$ is one solution; with $f(a,b,c)=6$. In order to find the other solutions I plug this system into Mathematica. After putting the pieces together, Mathematica simplifies this to a 48-degree equation, with only 1 additional solution satisfying the constraints: $a \approx 0.335288, b \approx 0.744425, c \approx 1.920288$. The precision does not matter much, since at this point the value of the function is $f(a,b,c) \approx 6.65093$, which means it's a local maximum.
</s>