This note presents a new approach to deriving stability conditions for continuous-time linear systems interconnected with a saturation. The method presented here can be extended to handle a deadzone, or in general, nonlinearities in the form of piecewise linear functions. By representing the saturation as a constrained optimization problem, the necessary (Kuhn-Tucker) conditions for optimality are used to derive linear and quadratic constraints which characterize the saturation. After selecting a candidate Lyapunov function, we pose the question of whether the Lyapunov function is decreasing along trajectories of the system as an implication between the necessary conditions derived from the saturation optimization, and the time derivative of the Lyapunov function. This leads to stability conditions in terms of linear matrix inequalities (LMIs), which are obtained by an application of the S-procedure to the implication. An example is provided where the proposed technique is compared and contrasted with previous analysis methods.