double population = 6.8e+7; const int Ndata = 153; double data_time[Ndata] = {45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192, 193, 194, 195, 196, 197}; double data_infected[Ndata] = {1, 1, 1, 3, 3, 7, 13, 17, 24, 38, 51, 62, 62, 116, 150, 202, 240, 274, 402, 554, 709, 927, 1170, 1187, 1280, 1326, 1353, 1380, 1462, 1505, 1585, 1655, 1686, 1749, 1845, 1934, 2003, 2028, 2173, 2272, 2415, 2506, 2605, 2783, 3034, 3158, 3300, 3465, 3635, 3953, 4220, 4361, 4546, 4793, 4996, 5350, 5647, 5951, 6336, 6783, 7220, 7572, 7808, 8232, 8895, 9420, 10015, 10652, 11350, 12074, 12739, 13524, 14355, 15515, 16433, 17200, 18003, 19137, 20125, 21343, 22583, 23615, 24264, 25937, 27403, 29240, 30967, 32683, 34357, 35812, 37525, 40792, 43434, 45973, 48285, 50879, 52991, 55421, 58568, 61927, 65736, 70038, 73533, 76334, 80412, 83890, 87715, 92681, 97302, 101590, 106108, 111796, 118375, 124590, 131800, 138134, 144264, 151209, 159333, 168061, 177124, 187977, 196750, 205721, 215855, 224665, 238339, 250687, 264184, 276242, 287796, 298292, 311049, 324221, 337594, 350879, 364328, 373628, 381798, 394948, 408052, 421996, 434200, 445433, 452529, 459761, 471123, 482169, 493183, 503290, 511485, 516862, 521318}; double data_active[Ndata]={1, 1, 1, 3, 3, 7, 13, 17, 24, 38, 51, 62, 62, 116, 150, 202, 240, 274, 402, 550, 697, 915, 1139, 1156, 1249, 1295, 1322, 1330, 1412, 1410, 1490, 1560, 1591, 1654, 1750, 1839, 1593, 1618, 1763, 1862, 2005, 2096, 1702, 1880, 2131, 2255, 2245, 2410, 2580, 2480, 2747, 2888, 3073, 3320, 2923, 3277, 3574, 3569, 3787, 4234, 4474, 4826, 4655, 5079, 5742, 5437, 5842, 6295, 6993, 7329, 7063, 7441, 7877, 8509, 9135, 9240, 9053, 10187, 10021, 11239, 11483, 11698, 11523, 12486, 13033, 14147, 14851, 15874, 17066, 17499, 17843, 19481, 20346, 21715, 23921, 24780, 23985, 23916, 25316, 26921, 28886, 31507, 33666, 34271, 36081, 38970, 39890, 42355, 45694, 48146, 51063, 54922, 58401, 60479, 64706, 69209, 73650, 77666, 83308, 86062, 90826, 96750, 103435, 107873, 113556, 117823, 125278, 132455, 136469, 141368, 149555, 152013, 150356, 158630, 159411, 168649, 173269, 178763, 173654, 165773, 171792, 176225, 171146, 180356, 177604, 172448, 173156, 172568, 167012, 160829, 164258, 158825, 157567}; double data_recovered[Ndata] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 12, 12, 31, 31, 31, 31, 31, 50, 50, 95, 95, 95, 95, 95, 95, 95, 410, 410, 410, 410, 410, 410, 903, 903, 903, 903, 1055, 1055, 1055, 1473, 1473, 1473, 1473, 1473, 2073, 2073, 2073, 2382, 2549, 2549, 2746, 2746, 3153, 3153, 3153, 3983, 4173, 4357, 4357, 4745, 5676, 6083, 6478, 7006, 7298, 7960, 8950, 8950, 10104, 10104, 11100, 11917, 12741, 13451, 14370, 15093, 16116, 16809, 17291, 18313, 19682, 21311, 23088, 24258, 24364, 26099, 29006, 31505, 33252, 35006, 36850, 38531, 39867, 42063, 44331, 44920, 47825, 50326, 51608, 53444, 55045, 56874, 59974, 64111, 67094, 68925, 70614, 73543, 76025, 81999, 86298, 91227, 93315, 97848, 102299, 106842, 113061, 118232, 127715, 134874, 138241, 146279, 160693, 165591, 178183, 182230, 191059, 194865, 208144, 229175, 236260, 245771, 263054, 265077, 274925, 287313, 297967, 309601, 326171, 342461, 347227, 358037, 363751}; double data_dead[Ndata] ={0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 2, 3, 5, 5, 5, 9, 9, 11, 12, 13, 18, 18, 24, 25, 25, 27, 27, 34, 48, 50, 52, 54, 58, 58, 65, 75, 79, 86, 87, 90, 93, 103, 103, 116, 123, 131, 138, 148, 153, 161, 178, 186, 194, 206, 206, 219, 238, 247, 261, 264, 286, 312, 339, 369, 397, 407, 429, 481, 524, 552, 577, 611, 643, 683, 705, 755, 792, 848, 908, 952, 998, 1080, 1162, 1210, 1284, 1354, 1423, 1480, 1568, 1625, 1674, 1737, 1831, 1877, 1930, 1991, 2102, 2205, 2292, 2340, 2413, 2456, 2529, 2657, 2749, 2844, 2952, 3026, 3199, 3310, 3502, 3602, 3720, 3860, 3971, 4079, 4172, 4346, 4453, 4669, 4804, 4948, 5033, 5173, 5368, 5940, 6093, 6343, 6655, 6769, 7067, 7257, 7497, 7812, 8005, 8153, 8366, 8539, 8884}; double ntests[Ndata]; /* double ntests[15]= [421, 2346, 21380, 17109, 17525, 24850, 41947, 75091, 88302, 132775, 144283, 146210, 123414, 119391, 114412]; */