Aproximative solutions to the neutrino oscillation problem in matter

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.