In this paper we classify the positive solutions of the divergent equation with Neumann boundary on the upper half space{-div(tα▽u) =tβf(u),(y,t) ∈ Rn+1+,limt→0+ tα(e)u/(e)t =0 by the method of moving spheres and Kelvin transformations,where n ≥ 1,α > 0,β >-1,n-1/n+1β ≤α < β + 2,and f :(0,∞) → (0,o∞) is non-negative continuous function satisfying some conditions.This equation arises from a weighed Sobolev inequality involving divergent operator div(tα▽u) on the upper half space.