In the paper, we study the plane Couette flow of a rarefied gas between two parallel infinite plates at y = ±L moving relative to each other with opposite velocities (±αL, 0, 0) along the x-direction. Assuming that the stationary state takes the specific form of F (y, vx −αy, vy, vz) with the x-component of the molecular velocity sheared linearly along the y-direction, such steady flow is governed by a boundary value problem on a steady nonlinear Boltzmann equation driven by an external shear force under the homogeneous non-moving diffuse reflection boundary condition. In case of the Maxwell molecule collisions, we establish the existence of spatially inhomogeneous non-equilibrium stationary solutions to the steady problem for any small enough shear rate α > 0 via an elaborate perturbation approach using Caflisch's decomposition together with Guo's L ∞ ∩ L 2 theory. The result indicates the polynomial tail at large velocities for the stationary distribution. Moreover, the large time asymptotic stability of the stationary solution with an exponential convergence is also obtained and as a consequence the nonnegativity of the steady profile is justified. Contents 1. Intoduction 1 2. Basic estimates 9 3. A trace theorem 12 4. Steady problem: the first order correction 15 5. Steady problem: remainder 22 5.1. Caflisch's decomposition 22 5.2. A priori estimates with parameters ǫ and σ 23 5.3. Existence for the linear problem with σ = 1 and ǫ > 0 32 5.4. Estimates on remainder 34 6. Unsteady problem: local existence 41 7. Unsteady problem: asymptotic stability and positivity 45 7.1. L ∞ estimates 46 7.2. L 2 estimates 49 8. Appendix 53 References 54 2020 Mathematics Subject Classification. 35Q20.