diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py
index d498b66..962516a 100644
--- a/src/seismic_hazard_forecasting.py
+++ b/src/seismic_hazard_forecasting.py
@@ -1,4 +1,9 @@
 # -*- coding: utf-8 -*-
+def datenum_to_datetime(datenum):
+    # MATLAB epoch starts at year 0, Python's datetime starts at year 1
+    matlab_epoch = datetime(1, 1, 1)
+    days = datenum - 367  # Subtracting the offset for MATLAB's epoch BUT THIS COULD BE PROBLEMATIC
+    return matlab_epoch + timedelta(days=days)
 
 def calc_bins2(dates, unit, bin_dur, dates_calc, rate_forecast, rate_unc_high, rate_unc_low, multiplicator, filename = "activity_rate.png", **kwargs):