Anovan() with higher squared factors

Hi, I have to conduct a statistic analysis on a set of data with 2 treatments. There is repeated measurements in the response variable. I managed to do the full linear analysis on the form of : yijk = 𝜇 + 𝜏 i + 𝛽 j + (𝜏𝛽)ij + 𝜖 ijk and I would like to perform a higher degree analysis.

I would like to conduct an anovan test where the individual quadratic effects of the treatments are considered.

Thanks!