| This research investigates analytic and numeric techniques for solving the Berreman system of ordinary differential equations for liquid crystal display optics in one spatial dimension. Modeling of the optical response of an individual system requires many solutions of the Berreman boundary value problem for various wavelengths and angles of incidence, thus making efficient solutions of the differential equation an important goal. We study techniques which, to some degree, overcome the complex oscillatory stiffness of the differential equation, described by the thickness k0 of the display cell in units of the light wavelength in the external medium. Various methods for solving the problem are compared, indicating both those computational regimes in which they perform favorably, and those situations where each experiences difficulties. |
