  We present neutrino mass matrices which predict the atmospheric neutrino
mixing to be almost maximal, $\sin^2 2\theta_{atm}>0.999$, as well as the large
solar neutrino mixing, $8/9>\sin^2 2\theta_{sol}>0.87$, and the large CP
violation (the CP violation phase in the standard form is maximal
$\delta=\pi/2$), based on the $Z_3$ symmetry.
