A tensor method is developed in this paper to deal with the propagation of non-Gaussian and non-symmetric optical beams. Three tensors are defined to describe any profile beams. Their transformation rules through non-symmetric systems are derived by using two-dimension Huygens integral. The beam quality factor of arbitrary beam is discussed. As an example, the transformation of high-order Hermite-Gaussian beam through a cylindrical lens system is given.