A novel numerical technique is presented to calculate the T-matrix for a single particle through the use of the volume integral equation for electromagnetic scattering. It is based on the method called Coupled Dipole Approximation (CDA), see O. J. F. Martin et al.¹. The basic procedure includes the parallel use of the Lippmann-Schwinger and the Dyson equations to iteratively solve for the T-matrix and the Green’s function dyadic respectively. The boundary conditions of the particle are thus automatically satisfied. The method can be used for the evaluation of the optical properties (e.g. Müller matrix) of anisotropic, inhomogeneous and asymmetric particles, both in far and near field, giving as output the T-matrix, which depends only on the scatterer itself and is independent from the polarization and direction of the incoming field. Estimation of the accuracy of the method is provided through comparison with the analytical spherical case (Mie theory) as well as non-spherical cubic ice particles.