An autonomous fourth order model of vocal fold vibrations is proposed. Each fold is represented by a lower and upper mass, and the aerodynamic forces are derived from a modified Bernoulli equation. The model exhibits many features of normal phonation in a wide parameter region. At the borderlines of this region coexistence of limit cycles, period-doubling and chaos are observed. Implications for an understanding of pathological voices are discussed.