Operador Laplaciano en coordenadas curvilíneas

El operador Laplaciano en coordenadas cartesianas es: $latex Delta u = frac{partial^2}{partial x^2} u + frac{partial^2}{partial y^2} u + frac{partial^2}{partial z^2} u = frac{partial^2}{(partial x^i)^2} u$. ¿Qué pasa cuando queremos expresarlo en otro sistema de coordenadas curvilineas tal que $latex x=x(q^1,q^2,q^3) = x(q^i)$, $latex y=y(q^i)$, $latex z = z(q^i)$ cualesquiera? Pues despues de un poco …