solving ordinary differential equations the newextended differential equations with the quantity become vsi psv psi vsi ir ( ( ( ( we shall refer to this model as the sirv model the new equation for poses no difficulties when it comes to the numerical method in forward euler scheme we simply add an update nc pts the program needs to store tin an additional array vand the plotting command must be extended with more arguments to plot versus as well the complete code is found in the file sirv py using : and : as values for the vaccine efficiency parameterthe effect of vaccination is seen in fig (other parameters are as in fig fig the effect of vaccinationp : (leftand : (right