Overall feature of CP dependence for neutrino oscillation probability in arbitrary matter profile

We study the CP dependence of neutrino oscillation probability for all channels in arbitrary matter profile within three generations. We show that an oscillation probability for ν_e→ν_μ can be written in the form P(ν_e→ν_μ)=A_eμcosδ+B_eμsinδ+C_eμ without any approximation using the CP phase δ. This result holds not only in constant matter but also in arbitrary matter. Another probability for ν_μ→ν_τ can be written in the form P(ν_μ→ν_τ)=A_μτcosδ+B_μτsinδ+C_μτ+D_μτcos2δ+E_μτsin2δ. The term which is proportional to sin2δ disappear, namely E_μτ=0, in symmetric matter. It means that the probability reduces to the same form as in constant matter. As for other channels, probabilities in arbitrary matter are at most the quadratic polynomials of sinδ and cosδ as in the above two channels. In symmetric matter, the oscillation probability for each channel reduces to the same form with respect to δ as that in constant matter.