This paper investigates the propagation of scalar and vector Hermite–Gaussian (HG) solitons in strong nonlocal media with exponential-decay response. The evolution equations for the parameters of the single HG beam are obtained by variational approach and the analytical results are confirmed by numerical simulation in section 2. Both the analytical and numerical solutions show that the critical power is increase with the increase of the order. Section 3 numerically studied the vector HG soliton and found that, the total critical power and the initial powers of the components should satisfy a complex relation. Because of the mutual attraction between the components, the stability of the quasi high-order vector HG soliton is better than the corresponding single HG beam during propagation.