We now return to model (6)- (9) and illustrate some of the properties discussed in the previous sections. The parameter values that we use for numerical simulations are in Table 1 [41] . We consider an imperfect flu vaccine for which the waning rate is about 0.15. The relationship between b 2 and b 3 is not very obvious; consequently, we can either assume equality or that the former is slightly greater than the latter. Transmission dynamics of infectious diseases with and without media coverage have already been carried out in previous studies, but these models do not account for the vaccination coverage. Therefore, we illustrate some numerical results for the model with optimal control when media coverage has (i) a beneficial effect (see Figures 2 and 3) (ii) and adverse effect on the vaccination rate (see Figures 4 and 5 ).