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.