A numerical procedure is used to calculate some of the internal free modes of oscillation in a two‐layer model of Lake Ontario, assuming a uniform equivalent depth. The modes fall into two categories, one set resembling Kelvin‐type waves and the other resembling Poincaré‐type waves. Observational evidence from Lake Ontario agrees qualitatively with the properties of these two types of modes.