I will consider classical optical systems where dynamics of the electromagnetic waves can be effectively described by equations that are at home in quantum mechanics. In these settings, the balanced gain and loss of light is reflected by the complex-valued, PT-symmetric refractive index. I will present the confluent Crum-Darboux transformation as an effective tool for construction of exactly solvable models. Its application will be illustrated on explicit examples.