Is this really correct for varying x?  
I'm not able to follow the algebra right now, but since trend and exog vary over time, my impression is that this should be a lag polynomial  (1-ar(L)) (trend_t + X_t*beta)  for conditional prediction, instead of (1-ar(1)) * (...)