Presented is the family of super-Gaussian spatial optical solitons in non-local nonlinear media modelled by 1+1 dimensional non-local nonlinear Schrdinger equation (NNLSE), in the strongly non-local case, an approximate analytical solution is obtained for an arbitrary response function by a variational approach. The solution with a sine beams width shows that the shift of phase, the spatial chirp and the critical input power are proportional to the order of the super-Gaussian optical beams, the period of the beam width is related not only to the material, but also to the orders of both the beams and the phase. The probability of quasi-square wave spatial solitons is also involved.