Problem:

Find the gradient of the function $f(x,y,z)$ where $f(x,y,z) = \sqrt{z \left(x^5+y\right)}$

from sympy.vector import *
from sympy import *

C = CoordSys3D('C')
x, y, z = C.x, C.y, C.z
f = sqrt(z*(x**5+y))
print(gradient(f))
