Book excerpt: Bipolar neutrosophic matrices (BNM) are obtained by bipolar neutrosophic sets. Each bipolar neutrosophic number represents an element of the matrix. The matrices are representable multi-dimensional arrays (3D arrays). The arrays have nested list data type. Some operations, especially the composition is a challenging algorithm in terms of coding because there are so many nested lists to manipulate. This paper presents a Python tool for bipolar neutrosophic matrices. The advantage of this work, is that the proposed Python tool can be used also for fuzzy matrices, bipolar fuzzy matrices, intuitionistic fuzzy matrices, bipolar intuitionistic fuzzy matrices and single valued neutrosophic matrices.