An EM algorithm implement (in Matlab) for estimating parameters of linear dynamical systems subjected to exogenous input