Vectors forms and partial derivatives?

Hella rusty folks! I'm trying to see if I can code this as a unsolved math function. Or would it be easier if I solved it on paper and then applied the written out formula?

For a scalar function ρ ( x , y , z ) = z^4 cos ( 3 π x ) and a vector v = ( 3 z , − 2 x , 4 y ) , derive the equation for ∇^2 ⋅ ρ. In your written work, write out the formula fully. Also, include a matlab code that allows you to solve for ∇^2 ⋅ ρ, given any values of [x,y,z]. Here enter the answer for [x=3,y=2,z=1]. Keep 3 sig figs.

Answer 1

[Anas Deiranieh's]

Ok, I figured it out.

laplacian works like the ∇^2 (del/nabla) operator and using syms this was do-able.