diff --git a/cara/montecarlo.py b/cara/montecarlo.py index d45acc4a..3e0c208e 100644 --- a/cara/montecarlo.py +++ b/cara/montecarlo.py @@ -344,8 +344,9 @@ class MCExposureModel: exposure = np.zeros(self.concentration_model.infected.samples) for start, stop in self.exposed.presence.boundaries(): - concentrations = np.asarray([self.concentration_model.concentration(t) for t in np.linspace(start, stop)]) - integrals = np.trapz(concentrations, axis=0) + times = np.arange(start, stop, 0.001) + concentrations = np.asarray([self.concentration_model.concentration(t) for t in times]) + integrals = np.trapz(concentrations, times, axis=0) exposure += integrals return exposure * self.repeats