Simultaneous multi-slice or multi-band (SMS/MB) imaging allows accelerated coverage in magnetic resonance imaging (MRI). Multiple slices are excited and acquired at the same time, and reconstructed using the redundancies in receiver coil arrays, similar to parallel imaging. SMS/MB reconstruction is currently performed with linear reconstruction techniques. Recently, a nonlinear reconstruction method for parallel imaging, Robust Artificial-neural-networks for k-space Interpolation (RAKI) was proposed and shown to improve upon linear methods. This method uses convolutional neural networks (CNN) trained solely on subject-specific calibration data. In this study, we sought to extend RAKI to SMS/MB imaging reconstruction. CNN training was performed on calibration data acquired prior to SMS/MB imaging, in a manner consistent with the existing linear methods. These CNNs were used to reconstruct a time series of functional MRI (fMRI) data. CNN network parameters were optimized using an extensive search of the parameter space. With these optimal parameters, RAKI substantially improves image quality compared to a commonly used linear reconstruction algorithm, especially for high acceleration rates.