Matlab Program

I would like to create a program that does the same simulation per attached published paper. I would not use the same bessel functions the researcher used. Therefore, I would like to calculate Q (calculated in appendix A) for a general pre defined function let's say B1 and B2. In other words, Q will be a subroutine that will call B1 and B2 instead of J and Y (as in the attached paper). General steps of programming -Creating a subroutine that calculates Q[j] for B1 and B2.. for any layer -Creating program that calculates the M matrix M= Q[j+]/Q[j-] -Creating another program that calculates T. The program will be compared by setting B1=J and B2=Y to get the same results of the paper.