We show how piecewise quadratic Lyapunov functions can be used to estimate regions of attraction for linear systems with saturation. The central issues of how to restrict the analysis region and how to optimize the “size” of the domain of attraction are addressed, and the approach is demonstrated on several examples. We observe that the piecewise quadratic Lyapunov functions yield significant improvements over recently proposed methods based on the Circle and Popov criteria.