Bipolar neutrosophic matrices (BNM) are obtained by bipolar neutrosophic sets. Eachbipolar neutrosophic number represents an element of the matrix. The matrices are representablemulti-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 manynested lists to manipulate. This paper presents a Python tool for bipolar neutrosophic matrices. Theadvantage of this work, is that the proposed Python tool can be used also for fuzzy matrices, bipolarfuzzy matrices, intuitionistic fuzzy matrices, bipolar intuitionistic fuzzy matrices and single valuedneutrosophic matrices.