The propagation characteristics of electromagnetic waves in thin films with nonlinear material are investigated numerically using variational methods and finite difference techniques. Both TE and TM polarizations are considered. It is observed that hysteresis and bistability can occur. It it also observed that part of the higher order solutions may be unstable and may follow a route to chaos through period doubling as the guided wave power increases.
