Coronaviruses are a group of RNA (ribonucleic acid) viruses with the capacity for rapid mutation and recombination. Coronaviruses are known to cause respiratory or intestinal infections in humans and animals. In this paper, a biologically compatible set of nonlinear fractional differential equations governing the outbreak of the novel coronavirus is suggested based on a model previously proposed in the literature. Then, this set is numerically solved utilizing two new methods employing sine–cosine and Bernoulli wavelets and their operational matrices. Moreover, the convergence of the solution is experimentally studied. Furthermore, the accuracy of the solution is proved via comparing the results with those obtained in previous research for the primary model. Furthermore, the computational costs are compared by measuring the CPU running time. Finally, the effects of the fractional orders on the outbreak of the COVID-19 are investigated.