  Neutrino flavor oscillations are discussed in terms of two coupled Dirac
fields. The interacting Lagrangian is diagonalized to obtain the exact
eigenvalues and eigenfunctions. Flavor wave functions are then derived directly
from the quantized neutrino fields. Probability density obtained by squaring
these wave functions upon taking into account the neutrino chirality are in
agreement with the standard neutrino oscillation probabilities.
