Trace[z, r] = T[z, r] + alphay * (T[z-1, r]-2*T[z, r]+T[z+1, r]) + alphax * (T[z, r-1]-2*T[z, r]+T[z, r+1]) if int(np.around(xm/dx)) <= r <= int(np.around(xM/dx ...