def create_dead_recover_reactions_model(
    vaccination_states,
    non_vaccination_state,
    virus_states,
    areas,
):
    """Create death and recover reactions for the model.

    Parameters
    ----------
    species : list of strings
        List containing the names of all species.

    vaccination_states : list of strings
        List containing the names of the vaccinations containing a state for
        non-vaccinated individuals.

    non_vaccination_state : str
        Name of state indicates non-vaccinated individuals.

    virus_states : list of strings
        List containing the names of the virus types.

    areas : list of strings
        List containing the names of the areas.

    Returns
    -------
    dead_recover_reactions: dict
        Dictionary that contains the dead/recover reaction names as keys
        and dicitonaries that contain the reactants, products and formulas
        as values.
    """
    dead_recover_reactions = {}
    for index_areas in areas:
        for index_vaccination in vaccination_states:
            for index_virus in virus_states:
                numb_infected = (
                    f"infectious_{index_areas}_{index_vaccination}_{index_virus}"
                )
                numb_recovered = (
                    f"recovered_{index_areas}_{index_vaccination}_{index_virus}"
                )
                numb_dead = f"dead_{index_areas}_{index_vaccination}_{index_virus}"
                key_recover = f"{numb_infected}_to_{numb_recovered}"
                key_death = f"{numb_infected}_to_{numb_dead}"

                if index_vaccination == non_vaccination_state:
                    omega_term = "0"
                else:
                    omega_term = f"omega_{index_vaccination}_{index_virus}"

                dead_recover_reactions[key_recover] = {
                    "reactants": {f"{numb_infected}": 1},
                    "products": {f"{numb_recovered}": 1},
                    "formula": f" (1 - (1 - {omega_term}) * prob_deceasing) * lambda1 * {numb_infected}",
                }
                dead_recover_reactions[key_death] = {
                    "reactants": {f"{numb_infected}": 1},
                    "products": {f"{numb_dead}": 1},
                    "formula": f"(1 - {omega_term}) * prob_deceasing * lambda1 * {numb_infected}",
                }
    return dead_recover_reactions