order-increasing forward layer recursion
[A,B] = frwrec(r)
number of layers
reflection coefficients = [r(1),...,r(M+1)]
(M+1)x(M+1) matrices whose columns are the reflection polynomials
frwrec(r) order-increasing forward layer recursion - from r to A,B
add sample pgm