Three numerical approaches to find mutually unbiased bases using Bell inequalities