  We present approximative solutions to the neutrino evolution equation
calculated by different methods. In a two neutrino framework, using the
physical parameters which gives the main effects to neutrino oscillations from
nu{e} to another flavors for L=3000Km and E=1GeV, the results for the
transition probability calculated by using series solutions, by to take the
neutrino evolution operator as a product of ordered partial operators and by
numerical methods, for a linearly and sinusoidally varying matter density are
compared. The extension to an arbitrary density profile is discussed and the
evolution operator as a product of partial operators in the three neutrino case
is obtained.
