Problem:
Find the laplacian of the function $f(x,y,z)$ where $f(x,y,z) = \sinh \left(x^2 z\right)$
from sympy.vector import *
from sympy import *

C = CoordSys3D('C')
x, y, z = C.x, C.y, C.z
f = sinh(x**2*z)
print(laplacian(f))
