A Multiregion Discrete-Time Epidemic Model of <i>Mycobacterium tuberculosis</i> Infections: An Optimal Control Approach

<div>Region of Oriental in Morocco. This region is divided into eight subregions (or domains)<span class="nowrap"><svg height="14.8173pt" id="M7" style="vertical-align:-5.52897pt" version="1.1" viewbox="-0.0498162 -9.28833 47.9173 14.8173" width="47.9173pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M300 -147C201 -63 143 98 143 270S200 602 300 686L282 710C136 610 70 450 70 271V270C70 89 136 -72 282 -170L300 -147Z"></path></g><g transform="matrix(.013,0,0,-0.013,4.498,0)"><path d="M719 146L691 147C666 77 657 72 612 72H528C572 125 688 224 688 364C688 546 558 665 384 665C212 665 77 541 77 361C77 221 188 124 239 72H150C110 72 97 80 73 147L46 146L63 0H307L314 18C284 71 178 203 178 360C178 529 256 631 383 631C509 631 588 523 588 363C588 197 485 67 453 18L458 0H701L719 146Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,14.339,3.132)"><path d="M400 606C400 634 383 656 353 656C316 656 294 620 294 593C294 564 317 545 343 545C375 545 400 573 400 606ZM366 351C379 413 381 451 356 451C323 451 251 408 183 341L199 313C223 335 267 365 277 365C285 365 284 354 277 312C245 132 222 27 193 -100C182 -148 160 -188 131 -188C113 -188 90 -178 75 -170C64 -164 55 -168 48 -175C38 -185 24 -203 24 -222S48 -257 71 -257C89 -257 131 -241 186 -192C243 -141 286 -46 310 74L366 351Z"></path></g><g transform="matrix(.013,0,0,-0.013,18.765,0)"><path d="M275 270C275 450 212 609 64 710L45 686C145 604 203 442 203 270S147 -63 45 -147L64 -170C213 -68 275 89 275 270Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,23.263,3.132)"><path d="M389 0V32C297 38 291 46 291 118V635C234 613 175 595 109 583V556L161 554C203 552 207 547 207 497V118C207 46 201 38 110 32V0H389Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,27.694,3.132)"><path d="M552 105V164L122 352L552 540V599L59 377V326L552 105ZM552 -75V-16L59 205V146L552 -75Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,33.275,3.132)"><path d="M400 606C400 634 383 656 353 656C316 656 294 620 294 593C294 564 317 545 343 545C375 545 400 573 400 606ZM366 351C379 413 381 451 356 451C323 451 251 408 183 341L199 313C223 335 267 365 277 365C285 365 284 354 277 312C245 132 222 27 193 -100C182 -148 160 -188 131 -188C113 -188 90 -178 75 -170C64 -164 55 -168 48 -175C38 -185 24 -203 24 -222S48 -257 71 -257C89 -257 131 -241 186 -192C243 -141 286 -46 310 74L366 351Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,37.189,3.132)"><path d="M552 105V164L122 352L552 540V599L59 377V326L552 105ZM552 -75V-16L59 205V146L552 -75Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,42.77,3.132)"><path d="M253 635C142 635 71 555 71 471C71 401 113 353 178 316C144 294 106 267 91 252C70 231 46 202 46 157C46 51 132 -12 241 -12C328 -12 442 52 442 168C442 257 379 305 311 343C354 373 380 398 388 407C406 428 417 458 417 487C417 568 348 635 253 635ZM240 600C288 600 340 563 340 481C340 423 313 387 279 360C207 394 148 427 148 499C148 549 182 600 240 600ZM252 23C188 23 129 73 129 163C129 217 162 266 211 298C289 260 360 217 360 144C360 69 312 23 252 23Z"></path></g></svg>:</span> Oujda Angad <span class="nowrap"><svg height="12.4698pt" id="M8" style="vertical-align:-3.18147pt" version="1.1" viewbox="-0.0498162 -9.28833 23.9537 12.4698" width="23.9537pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M300 -147C201 -63 143 98 143 270S200 602 300 686L282 710C136 610 70 450 70 271V270C70 89 136 -72 282 -170L300 -147Z"></path></g><g transform="matrix(.013,0,0,-0.013,4.498,0)"><path d="M719 146L691 147C666 77 657 72 612 72H528C572 125 688 224 688 364C688 546 558 665 384 665C212 665 77 541 77 361C77 221 188 124 239 72H150C110 72 97 80 73 147L46 146L63 0H307L314 18C284 71 178 203 178 360C178 529 256 631 383 631C509 631 588 523 588 363C588 197 485 67 453 18L458 0H701L719 146Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,14.339,3.132)"><path d="M389 0V32C297 38 291 46 291 118V635C234 613 175 595 109 583V556L161 554C203 552 207 547 207 497V118C207 46 201 38 110 32V0H389Z"></path></g><g transform="matrix(.013,0,0,-0.013,19.285,0)"><path d="M275 270C275 450 212 609 64 710L45 686C145 604 203 442 203 270S147 -63 45 -147L64 -170C213 -68 275 89 275 270Z"></path></g></svg>,</span> Berkane <span class="nowrap"><svg height="12.4698pt" id="M9" style="vertical-align:-3.18147pt" version="1.1" viewbox="-0.0498162 -9.28833 23.9537 12.4698" width="23.9537pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M300 -147C201 -63 143 98 143 270S200 602 300 686L282 710C136 610 70 450 70 271V270C70 89 136 -72 282 -170L300 -147Z"></path></g><g transform="matrix(.013,0,0,-0.013,4.498,0)"><path d="M719 146L691 147C666 77 657 72 612 72H528C572 125 688 224 688 364C688 546 558 665 384 665C212 665 77 541 77 361C77 221 188 124 239 72H150C110 72 97 80 73 147L46 146L63 0H307L314 18C284 71 178 203 178 360C178 529 256 631 383 631C509 631 588 523 588 363C588 197 485 67 453 18L458 0H701L719 146Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,14.339,3.132)"><path d="M414 144C384 79 371 75 317 75H135L276 221C367 316 408 376 408 465C408 570 327 635 237 635C179 635 131 609 100 575L42 494L67 471C94 510 138 565 205 565C277 565 321 517 321 435C321 348 258 270 195 195C146 137 88 81 33 26V0H411C423 44 433 88 446 135L414 144Z"></path></g><g transform="matrix(.013,0,0,-0.013,19.285,0)"><path d="M275 270C275 450 212 609 64 710L45 686C145 604 203 442 203 270S147 -63 45 -147L64 -170C213 -68 275 89 275 270Z"></path></g></svg>,</span> Nador <span class="nowrap"><svg height="12.5794pt" id="M10" style="vertical-align:-3.29107pt" version="1.1" viewbox="-0.0498162 -9.28833 23.9537 12.5794" width="23.9537pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M300 -147C201 -63 143 98 143 270S200 602 300 686L282 710C136 610 70 450 70 271V270C70 89 136 -72 282 -170L300 -147Z"></path></g><g transform="matrix(.013,0,0,-0.013,4.498,0)"><path d="M719 146L691 147C666 77 657 72 612 72H528C572 125 688 224 688 364C688 546 558 665 384 665C212 665 77 541 77 361C77 221 188 124 239 72H150C110 72 97 80 73 147L46 146L63 0H307L314 18C284 71 178 203 178 360C178 529 256 631 383 631C509 631 588 523 588 363C588 197 485 67 453 18L458 0H701L719 146Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,14.339,3.132)"><path d="M290 377C321 398 342 415 358 430C378 450 389 473 389 502C389 578 329 635 238 635H237C184 635 137 610 109 578L64 515L88 493C112 529 154 573 208 573S303 542 303 482C303 409 233 370 141 341L149 308C165 313 190 319 215 319C272 319 341 283 341 193C342 98 292 43 222 43C163 43 122 72 96 94C88 101 79 100 70 94C61 87 47 73 46 60C44 47 48 37 62 23C76 10 118 -12 165 -12C238 -12 430 62 430 223C430 297 379 359 290 375V377Z"></path></g><g transform="matrix(.013,0,0,-0.013,19.285,0)"><path d="M275 270C275 450 212 609 64 710L45 686C145 604 203 442 203 270S147 -63 45 -147L64 -170C213 -68 275 89 275 270Z"></path></g></svg>,</span> Driouch <span class="nowrap"><svg height="12.4698pt" id="M11" style="vertical-align:-3.18147pt" version="1.1" viewbox="-0.0498162 -9.28833 23.9537 12.4698" width="23.9537pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M300 -147C201 -63 143 98 143 270S200 602 300 686L282 710C136 610 70 450 70 271V270C70 89 136 -72 282 -170L300 -147Z"></path></g><g transform="matrix(.013,0,0,-0.013,4.498,0)"><path d="M719 146L691 147C666 77 657 72 612 72H528C572 125 688 224 688 364C688 546 558 665 384 665C212 665 77 541 77 361C77 221 188 124 239 72H150C110 72 97 80 73 147L46 146L63 0H307L314 18C284 71 178 203 178 360C178 529 256 631 383 631C509 631 588 523 588 363C588 197 485 67 453 18L458 0H701L719 146Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,14.339,3.132)"><path d="M462 177V227H365V632H320C217 496 116 350 21 208V177H284V109C284 43 280 38 190 31V0H451V31C369 38 365 43 365 108V177H462ZM284 227H88C155 336 218 430 282 519H284V227Z"></path></g><g transform="matrix(.013,0,0,-0.013,19.285,0)"><path d="M275 270C275 450 212 609 64 710L45 686C145 604 203 442 203 270S147 -63 45 -147L64 -170C213 -68 275 89 275 270Z"></path></g></svg>,</span> Guercif <span class="nowrap"><svg height="12.5794pt" id="M12" style="vertical-align:-3.29107pt" version="1.1" viewbox="-0.0498162 -9.28833 23.9537 12.5794" width="23.9537pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M300 -147C201 -63 143 98 143 270S200 602 300 686L282 710C136 610 70 450 70 271V270C70 89 136 -72 282 -170L300 -147Z"></path></g><g transform="matrix(.013,0,0,-0.013,4.498,0)"><path d="M719 146L691 147C666 77 657 72 612 72H528C572 125 688 224 688 364C688 546 558 665 384 665C212 665 77 541 77 361C77 221 188 124 239 72H150C110 72 97 80 73 147L46 146L63 0H307L314 18C284 71 178 203 178 360C178 529 256 631 383 631C509 631 588 523 588 363C588 197 485 67 453 18L458 0H701L719 146Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,14.339,3.132)"><path d="M158 548H390L417 615L410 623H122L83 318C105 326 143 337 185 337C296 337 350 275 350 188C350 116 308 42 225 42C164 42 122 74 100 93C90 101 82 99 72 92C60 82 51 68 50 59C48 46 52 38 66 24C82 9 125 -12 172 -12C225 -11 292 15 346 59C408 108 437 166 437 226C437 309 371 397 242 397C214 397 170 382 133 369L158 548Z"></path></g><g transform="matrix(.013,0,0,-0.013,19.285,0)"><path d="M275 270C275 450 212 609 64 710L45 686C145 604 203 442 203 270S147 -63 45 -147L64 -170C213 -68 275 89 275 270Z"></path></g></svg>,</span> Taourirt <span class="nowrap"><svg height="12.5794pt" id="M13" style="vertical-align:-3.29107pt" version="1.1" viewbox="-0.0498162 -9.28833 23.9537 12.5794" width="23.9537pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M300 -147C201 -63 143 98 143 270S200 602 300 686L282 710C136 610 70 450 70 271V270C70 89 136 -72 282 -170L300 -147Z"></path></g><g transform="matrix(.013,0,0,-0.013,4.498,0)"><path d="M719 146L691 147C666 77 657 72 612 72H528C572 125 688 224 688 364C688 546 558 665 384 665C212 665 77 541 77 361C77 221 188 124 239 72H150C110 72 97 80 73 147L46 146L63 0H307L314 18C284 71 178 203 178 360C178 529 256 631 383 631C509 631 588 523 588 363C588 197 485 67 453 18L458 0H701L719 146Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,14.339,3.132)"><path d="M141 347C171 480 264 541 327 570C364 587 403 599 436 605L429 641C389 634 339 622 300 609C191 570 37 458 37 239C37 86 126 -12 245 -12C367 -12 454 89 454 210C454 314 380 397 272 397C252 397 230 390 208 380L141 347ZM231 338C323 338 366 257 366 174C366 108 340 27 263 27C180 27 129 123 129 244C129 268 130 292 135 310C158 323 193 338 231 338Z"></path></g><g transform="matrix(.013,0,0,-0.013,19.285,0)"><path d="M275 270C275 450 212 609 64 710L45 686C145 604 203 442 203 270S147 -63 45 -147L64 -170C213 -68 275 89 275 270Z"></path></g></svg>,</span> Jerada <span class="nowrap"><svg height="12.5794pt" id="M14" style="vertical-align:-3.29107pt" version="1.1" viewbox="-0.0498162 -9.28833 23.9537 12.5794" width="23.9537pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M300 -147C201 -63 143 98 143 270S200 602 300 686L282 710C136 610 70 450 70 271V270C70 89 136 -72 282 -170L300 -147Z"></path></g><g transform="matrix(.013,0,0,-0.013,4.498,0)"><path d="M719 146L691 147C666 77 657 72 612 72H528C572 125 688 224 688 364C688 546 558 665 384 665C212 665 77 541 77 361C77 221 188 124 239 72H150C110 72 97 80 73 147L46 146L63 0H307L314 18C284 71 178 203 178 360C178 529 256 631 383 631C509 631 588 523 588 363C588 197 485 67 453 18L458 0H701L719 146Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,14.339,3.132)"><path d="M453 623H65C60 580 56 530 46 472H80C103 538 110 548 175 548H389C310 380 195 168 90 0L98 -12L175 -2C273 206 366 408 462 610L453 623Z"></path></g><g transform="matrix(.013,0,0,-0.013,19.285,0)"><path d="M275 270C275 450 212 609 64 710L45 686C145 604 203 442 203 270S147 -63 45 -147L64 -170C213 -68 275 89 275 270Z"></path></g></svg>,</span> and Figuig <span class="nowrap"><svg height="12.5794pt" id="M15" style="vertical-align:-3.29107pt" version="1.1" viewbox="-0.0498162 -9.28833 23.9537 12.5794" width="23.9537pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M300 -147C201 -63 143 98 143 270S200 602 300 686L282 710C136 610 70 450 70 271V270C70 89 136 -72 282 -170L300 -147Z"></path></g><g transform="matrix(.013,0,0,-0.013,4.498,0)"><path d="M719 146L691 147C666 77 657 72 612 72H528C572 125 688 224 688 364C688 546 558 665 384 665C212 665 77 541 77 361C77 221 188 124 239 72H150C110 72 97 80 73 147L46 146L63 0H307L314 18C284 71 178 203 178 360C178 529 256 631 383 631C509 631 588 523 588 363C588 197 485 67 453 18L458 0H701L719 146Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,14.339,3.132)"><path d="M253 635C142 635 71 555 71 471C71 401 113 353 178 316C144 294 106 267 91 252C70 231 46 202 46 157C46 51 132 -12 241 -12C328 -12 442 52 442 168C442 257 379 305 311 343C354 373 380 398 388 407C406 428 417 458 417 487C417 568 348 635 253 635ZM240 600C288 600 340 563 340 481C340 423 313 387 279 360C207 394 148 427 148 499C148 549 182 600 240 600ZM252 23C188 23 129 73 129 163C129 217 162 266 211 298C289 260 360 217 360 144C360 69 312 23 252 23Z"></path></g><g transform="matrix(.013,0,0,-0.013,19.285,0)"><path d="M275 270C275 450 212 609 64 710L45 686C145 604 203 442 203 270S147 -63 45 -147L64 -170C213 -68 275 89 275 270Z"></path></g></svg>.</span></div>

Discrete Dynamics in Nature and Society

fig1

Figure 1

Figure 1: A Multiregion Discrete-Time Epidemic Model of <i>Mycobacterium tuberculosis</i> Infections: An Optimal Control Approach 

Discrete Dynamics in Nature and Society

A Multiregion Discrete-Time Epidemic Model of Mycobacterium tuberculosis Infections: An Optimal Control Approach

Figure 1