In this paper we present a method to recover a time-homogeneous piecewise constant volatility from a finite set of perpetual put option prices. The whole calculation process of the volatility is decomposed into easy computations in many fixed disjoint intervals. In each interval, the volatility is obtained by solving a system of nonlinear equations.