fastcpd_ar function

Find change points efficiently in AR($p$) models