Meromorphic solutions of non-linear differential equations of the form f n + P (z, f ) = h are investigated, where n ≥ 2 is an integer, h is a meromorphic function, and P (z, f ) is differential polynomial in f and its derivatives with small functions as its coefficients. In the existing literature this equation has been studied in the case when h has the particular form h(z) = p 1 (z)e α1(z) + p 2 (z)e α2 (z) , where p 1 , p 2 are small functions of f and α 1 , α 2 are entire functions. In such a case the order of h is either a positive integer or equal to infinity. In this article it is assumed that h is a meromorphic solution of the linear differential equation h ′′ + r 1 (z)h ′ + r 0 (z)h = r 2 (z) with rational coefficients r 0 , r 1 , r 2 , and hence the order of h is a rational number. Recent results by Liao-Yang-Zhang (2013) and Liao (2015) follow as special cases of the main results.