Code for numerical experiments using BAOAB with constraints