matlab coding of threshold parameter R*

i need to write a series of codes which computes the threshold parameter R*. as it stands i am completely lost at what next to do...do find my current matlab script . rather than doing each beta, how do i write a code that does it recursively.. hope that makes sense..now how do i extend that to find the Mu's and in turn R*. if an answer cant be given, tips would be appreciated....... find attached files for detailed explanation