get_Bt_Cpp function

Calculate the impact matrix $B_t$ for all $t$ for models with a non-Gaussian conditional distribution with mutually independent shocks.