diff --git a/matlab/imcforecast.m b/matlab/imcforecast.m
index e71d2971624fe73c69a9298e87ca588d2c654ce5..bc7c0a903bfcfacd4303f4c08ea4c7835df07944 100644
--- a/matlab/imcforecast.m
+++ b/matlab/imcforecast.m
@@ -180,7 +180,9 @@ else
     cL = size(constrained_paths,2);
 end
 
-constrained_paths = bsxfun(@minus,constrained_paths,trend(idx,1:cL));
+for t=1:cL
+    constrained_paths(:,t) = constrained_paths(:,t)-trend(idx,t);
+end
 
 randn('state',0);