A solution for the elastic buckling of flat rectangular plates with centerline boundary conditions subjected to non-uniform in-plane axial compression is presented. The loaded edges are simply supported, the non-loaded edges are free, and the centerline is simply supported with a variable rotational stiffness. The Galerkin method is used to establish an eigenvalue problem and a series solution for plate buckling coefficients is obtained by using combined trigonometric and polynomial functions that satisfy the boundary conditions. It is demonstrated that the formulation approaches the classical solution of a plate with a clamped edge as the variable rotational stiffness is increased. The variation of buckling coefficient with aspect ratio is presented for various stress gradient ratios. The coupling between plate aspect ratio, centerline rotational stiffness, and gradient of applied compressive stress is illustrated and discussed. The solution is applicable to stiffened plates and I-shaped beams that are subjected to biaxial bending or combined flexure and torsion, and is important to estimate the reduction in elastic buckling capacity due to stress gradient. © 2005 Elsevier Ltd. All rights reserved.