Optical diffraction tomography (ODT) reconstructs a samples volumetric refractive index (RI) to create high-contrast, quantitative 3D visualizations of biological samples. However, standard implementations of ODT use interferometric systems, and so are sensitive to phase instabilities, complex mechanical design, and coherent noise. Furthermore, their reconstruction framework is typically limited to weakly scattering samples, and thus excludes a whole class of multiple-scattering samples. Here, we implement a new 3D RI microscopy technique that utilizes a computational multi-slice beam propagation method to invert the optical scattering process and reconstruct high-resolution (NA > 1.0) 3D RI distributions of multiple-scattering samples. The method acquires intensity-only measurements from different illumination angles and then solves a nonlinear optimization problem to recover the samples 3D RI distribution. We experimentally demonstrate the reconstruction of samples with varying amounts of multiple-scattering: a 3T3 fibroblast cell, a cluster of C. elegans embryos, and a whole C. elegans worm, with lateral and axial resolutions of ≤ 240 nm and ≤ 900 nm, respectively. The results of this work lays groundwork for future studies into using optical wavelengths to probe 3D RI distributions of highly scattering biological organisms.