The authors consider the differential operator \[ D_{q,p}=(-1)^{q- p}z\prod^{p}_{j=1}(\vartheta +\mu_ j)-\prod^{q}_{j=1}(\vartheta +\nu_ j-1), \] where \(\vartheta\) stands for the Euler operator zd/dz; \(\mu_ 1,\mu_ 2,...,\nu_ q\) are given complex parameters while p and q are certain natural numbers such that \(\mu_ i-\mu_ j\) are not entire, \(p\geq 1\) and \(q\geq p+1\). Their main aim is to establish, in each one of certain appropriately designed consecutive sectors, a set of fundamental solutions to \(D_{q,p}u=0\) such that the asymptotic expansions of these solutions, valid near the point \(z=\infty\), are prescribed. In the first section of the paper they study the matrices (named: Stokes Matrices) which permit one to translate the system of solution valid in a sector to the solutions in the next sector over the intersections of these sectors. To this end they transform first the given equation through the Mellin transform and then use the classical inverse transform integral. The second section of the paper is devoted to the investigation of the group properties of the operators \(D_{3,1}\) and \(D_{3,2}\).