Application and Optimal Control for an HBV Model with Vaccination and Treatment

<div>Contour plots of the basic reproductive number <svg height="12.3054pt" id="M202" style="vertical-align:-3.29108pt" version="1.1" viewbox="-0.0498162 -9.01432 17.3641 12.3054" width="17.3641pt" 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="M852 669C820 657 792 643 768 628C705 667 627 687 547 687C310 687 154 548 154 412C154 332 212 282 296 282C434 282 515 426 523 557H503C474 402 410 312 298 312C241 312 195 352 195 408C195 535 341 657 546 657C613 657 683 644 738 609C640 538 574 431 486 271C374 71 314 16 189 16C118 16 69 60 69 95H71C78 84 93 83 101 83C129 83 144 106 144 131C144 160 121 179 93 179C60 179 34 151 34 108C34 42 100 -15 207 -15C364 -15 491 73 603 304C614 297 639 293 650 293C681 293 683 254 683 245C683 181 599 150 599 70C599 21 645 -12 692 -12C747 -12 804 30 868 103L851 119C789 52 741 18 706 18C685 18 674 29 674 52C674 92 750 157 750 226C750 257 727 287 696 310V312C812 315 876 394 876 473C876 525 849 577 793 612C812 627 832 639 859 650L852 669ZM814 486C814 417 765 345 686 337C669 348 649 354 626 354C673 451 708 533 766 588C799 561 814 525 814 486Z" id="g198-19"></path></g><g transform="matrix(.0091,0,0,-0.0091,12.272,3.132)"><path d="M245 635C92 635 37 457 37 312C37 149 91 -12 244 -12C395 -12 449 166 449 312C449 469 395 635 245 635ZM243 598C332 598 358 454 358 312C358 173 334 26 245 26C158 26 128 174 128 313S152 598 243 598Z" id="g50-49"></path></g></svg> versus the parameter pairs (<span class="nowrap"><svg height="6.34998pt" id="M203" style="vertical-align:-0.2063899pt" version="1.1" viewbox="-0.0498162 -6.14359 7.47218 6.34998" width="7.47218pt" 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="M548 455L523 469C505 448 490 439 461 439C416 439 374 447 316 447C144 447 23 310 23 161C23 47 89 -12 179 -12C312 -12 429 106 429 244C429 303 419 344 368 386L371 388C404 383 440 378 474 378C507 378 530 411 548 455ZM350 274C350 191 308 25 198 25C144 25 108 75 108 157C108 264 170 395 276 395C296 395 320 384 333 360C347 334 350 316 350 274Z" id="g113-240"></path></g></svg>,</span><svg height="6.1673pt" id="M204" style="vertical-align:-0.2063904pt" version="1.1" viewbox="-0.0498162 -5.96091 8.42472 6.1673" width="8.42472pt" 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="M615 290C615 392 563 448 529 448C511 448 491 441 484 434L483 425C518 400 546 349 546 273C546 158 504 46 417 46C381 46 356 77 356 142C356 184 380 293 392 334L387 340L307 323C305 288 296 228 287 188C265 94 221 52 180 52C144 52 112 89 112 171C112 278 164 372 241 430L226 454C106 394 23 288 23 155C23 40 79 -12 139 -12C184 -12 245 32 285 85C292 31 322 -12 376 -12C423 -12 470 18 501 45C570 104 615 184 615 290Z" id="g113-246"></path></g></svg>) (left) and (<svg height="12.4698pt" id="M205" style="vertical-align:-3.18147pt" version="1.1" viewbox="-0.0498162 -9.28833 21.7223 12.4698" width="21.7223pt" 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="M475 507C475 612 440 712 326 712C139 712 23 420 23 215C23 96 58 -12 180 -12C369 -12 475 293 475 507ZM391 522C391 486 387 448 379 394H126C155 538 222 677 310 677C386 677 391 571 391 522ZM373 346C344 193 283 22 189 22C126 22 106 114 106 196C106 243 111 293 118 346H373Z" id="g113-230"></path></g><g transform="matrix(.0091,0,0,-0.0091,6.149,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" id="g50-51"></path></g><g transform="matrix(.013,0,0,-0.013,11.095,0)"><path d="M95 130C70 130 46 113 46 88C46 72 54 64 59 64C93 55 121 33 121 -3C121 -41 93 -68 44 -88L55 -117C117 -98 186 -56 186 22C186 91 131 130 95 130Z" id="g113-45"></path></g><g transform="matrix(.013,0,0,-0.013,16.239,0)"><path d="M387 375C387 402 357 448 257 448C172 448 82 404 82 326C82 289 108 255 156 241V239C85 223 23 181 23 116C23 39 89 -12 182 -12C265 -12 336 31 378 91L361 114C320 73 269 47 216 47C157 47 115 82 115 137C115 191 160 219 218 219C243 219 262 218 272 217L304 259L302 266C295 265 281 264 255 264C195 264 163 294 163 335C163 377 200 416 249 416C293 416 321 389 329 342C331 332 335 329 341 329C355 329 387 352 387 375Z" id="g113-247"></path></g></svg>) (right).</div>

Discrete Dynamics in Nature and Society

fig3

Figure 3

Figure 3: Application and Optimal Control for an HBV Model with Vaccination and Treatment 