Let $$\mu _\alpha$$ be the Lebesgue plane measure on the unit disk with the radial weight $$\frac{\alpha +1}{\pi }(1-|z|^2)^\alpha$$ . Denote by $${\mathcal {A}}^{2}_{n}$$ the space of the n-analytic functions on the unit disk $${\mathbb {D}}$$ , square-integrable with respect to $$\mu _\alpha$$ . Extending results of Ramazanov (1999, 2002), we explain that disk polynomials (studied by Koornwinder in 1975 and Wünsche in 2005) form an orthonormal basis of $${\mathcal {A}}^{2}_{n}$$ . Using this basis, we provide the Fourier decomposition of $${\mathcal {A}}^{2}_{n}$$ into the orthogonal sum of the subspaces associated with different frequencies. This leads to the decomposition of the von Neumann algebra of radial operators, acting in $${\mathcal {A}}^{2}_n$$ , into the direct sum of some matrix algebras. In other words, all radial operators are represented as matrix sequences. In particular, we represent in this form the Toeplitz operators with bounded radial symbols, acting in $${\mathcal {A}}^{2}_n$$ . Moreover, using ideas by Engliš (1996), we show that the set of the Toeplitz operators with bounded generating symbols is not weakly dense in $${\mathcal {B}}({\mathcal {A}}^{2}_n)$$ .