Mathematical tools called differential and integral operators are used to model real world problems in all fields of science as they are able to replicate some behaviors observed in real world like fading memory, long-range dependency, power law, random walk and many others. Very recently the world has faced a serious challenge since the breakout of corona-virus started in Wuhan, China. The deathly disease has killed about 1720000 and infected more than 2 millions humans around the globe since December 2019 to 21 of April 2020. In this paper, we analyzed a mathematical model for the spread of COVID-19, we first start with stability analysis, present the optimal control for the system. The model was extended to the concept of non-local operators for each case, we presented the positiveness of the system solutions. We presented numerical solutions are presented for different scenarios.