Multispectral bioluminescence tomography (BLT) attracts increasinglymore attention in the area of optical molecular imaging. In this paper, we analyzethe properties of the solutions to the regularized and discretized multispectralBLT problems. First, we show the solution existence, uniqueness, and itscontinuous dependence on the data. Then, we introduce stable numerical schemes andderive error estimates for numerical solutions. We report some numerical resultsto illustrate the performance of the numerical methods on the quality of multispectralBLT reconstruction.