@jakubi  Dear friend,

I think this is a good answer. I'm trying to calculate the stress tensor divergence (momentum equation) in curvilinear coordinates, however the results that I get in Maple, using your procedure, are a little different, comparing with some literature. I can`t understand why. I `m calculating the connection using Christoffel from Differential Geometry package.

This is the tensor

This my results:

This is the result from literature:

I really appreciate your help