To solve the problem of light scattering by multilayer scatterers of an arbitrary axisymmetric shape, a separation of variables method that involves special scalar potentials and their expansions in spherical functions is developed. The approach is shown to yield highly exact results even for particles that have 100 layers or more. A graphic library that illustrates the optical properties of layered and homogeneous (with an effective refractive index) spheroids, spheres, and Chebyshev particles of various shapes and sizes (about 650 figures) is created and is put on the Internet. It is noted that the linear polarization of radiation transmitted forward through a polydisperse medium containing partially oriented nonspherical porous particles strongly depends on the structure of scatterers. It is shown that the difference between the degrees of polarization of layered and corresponding homogeneous scatterers can exceed 200-300%.