This is a kinda p-recursive method. So, when polynomials of some order, and many repetitions are needed, one can use this method.
This is one of the most accurate methods, especially when high orders are needed. It even has a better accuracy than the q-recursive. This will appear in our paper:
Title: NUMERICAL-STABILITY QUALITY FACTOR FOR ORTHOGONAL POLYNOMIALS: ZERNIKE
RADIAL POLYNOMIALS CASE STUDY
that will be published soon Applied Numerical Mathematics
Few changes, but the code is mainly similar.