Problem:

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

from sympy.vector import *
from sympy import *

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