diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py
index fc758ac..8cdd3a1 100644
--- a/src/seismic_hazard_forecasting.py
+++ b/src/seismic_hazard_forecasting.py
@@ -1,12 +1,7 @@
 # -*- coding: utf-8 -*-
-from datetime import timedelta, datetime
 import numpy as np
+from igfash.rate import apply_beast, datenum_to_datetime
 
-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):