An algorithm based on со variant (intrinsic) paeudodifferential calculus for computing DeWitt-Seeley-Gilkey coefficient* of elliptic operators on compact curved manifold with torsion is considered. To implement this algorithm a program package for algebraic and tensor manipulations has been written in С language. Some computer results are presented.