Here's my mostly pyhthonic model, actually it's done in sagecell.sagemath.org. Problem is, how do you predict the parameters alpha, probability if infection, and beta, probability of recovery. My simulation us interactive in that you can change the values of alpha and beta. A higher alpha ir lower beta with skew the infected curve so it peaks higher and sooner. So how do we reduce alpha? We use vaccines we don't have. How do we increase beta? We use antivirals we don't have. All we can do now us stay home and stay apart to reduce alpha and flatten the infected curve! Stay Healthy, A. Jorge García https://sagecell.sagemath.org/?z=eJytlEtv2zAMgO_-FUSKwXbitE67YVgXr4cdBt8K-9jmoNhMrMGxPIlO0P36UX4l9R6HbQ1QEXx8FCnSplAnz71a-ZDjTlay2kMaJ3BQOZZgmqwAKgRBukgWcbQKwzegdqxCqFXdlIKkqlzfMSPFvBjCg3XK5W6H3xrTsSyYCo0IsiLUOdZY5VjRBccw6Ci054qyLgRskQSkkEDM-jzNySM_Wra2eTqPP3Y508ZkWJPclgiPIwoetaqVtmIAeXqTU-QGPYNhcQcbWUuba0TG1Q4zwvy3vLjnxQMv6XivKAlm6oj6D5ikxyQ9povjYod23vpwKrCCkyxLiAPuXH-x-gJjSGjix8s0CsNdfjg_B8drNE1J95DOxXK7DiGKPkG63t4I97-1YScyUtobujFkN6o84qBdhwH59-AGU236c92zqzu-d1OBkYchr2fj7Ai1OhIVKh6sccKOUgA2JWqXpw2pULk_8x2KQqetLQqvw7vQsdVZMXznmOgpvP7wYeNIK4SrjaNbIdwMveM7dLFu0J6BC18bQyBg36AxAbAb5ApeVMPPgbnMCJTm7MI0GnnUpXl4uHgL5rX53cAef02rNa-PN5wuPROYJ9o8k_2BPIuaxdGfAveZuPlP4aaT5Chplnxnx8k4vAItqj16q-D9bejfO8B_vLsk0mHxONlytZnL9jjb4-gX5m6cpq5JNFWba1Hbr4HXRS-6jH5rk4NNXtjizqYHm76wJZ1tUjeNddNYN23O7Rz6efXWh7pUdDF7YxNbYwrS8KuVai8NyczunXiBrNFHHB1n1jG2jl_sF6Di3n5u1EHxcJ7EEWevgMkEuNfqRMWE6HbEfsnRdJ-FdN0N082_DShF3aOH7aM7afRd1tw54ztxL0rfSXpRs4NtUGT_ebWyt7M7HE-VMSuTqTIZtr1lLNqgRevl_wDu9Rsw&lang=sage&interacts=eJyLjgUAARUAuQ