This paper proposes a probabilistic approach for the solution of elastic buckling of columns, involving uncertainties, using stochastic element free Galerkin method. In the present work, modulus of elasticity is modeled as a homogeneous random field. Karhunen-Loeve expansion and shape function method are used to represent random field and their effectiveness is compared in modeling the same in a computationally viable manner. Both Gaussian and non-Gaussian field are considered for the present study. The stochastic eigenvalue problem is solved for first and second moment characteristics of buckling load, using perturbation analysis. Numerical examples of columns with different boundary conditions are solved. Monte Carlo simulation is used as a validation tool. The obtained results are found in good agreement with those obtained by Monte Carlo simulation. © 2017 Elsevier Ltd