How can I implement Newton Raphson Method using as initial guess a 3d matrix (mxmxm) in Matlab?

Hello everyone! I made a code to evaluate temperature field in a cube mxmxm. I developed all equations for each boundary condition and then implement a single scan in each point of my cube (I see my cube as a 3d matrix), using "for" control statement, in i,j,k coordinates. Now, I need converge this initial values, and I think NR method is a good guess. Is there a way to implement NR method in Matlab in order to perform this?