def METHOD_NAME(self):
    # Initial concentrations
    self.solver.run(initials=self.y0)
    print(self.solver.vol)
    assert self.solver.vol is None
    self.solver.vol = 1e-20
    assert self.solver.vol == 1e-20