The flow between eccentric rotating cylinders is studied numerically using a curvilinear coordinate system. The outer cylinder is fixed and the inner cylinder rotates at a constant speed. The three-dimensional Navier-Stokes equations and the equation of continuity are solved. A body-fitted grid system is generated and the basic equations in a physical domain are transformed into a calculational domain and solved by means of a finite difference method. A series of computations is carried out at Re=130 to 500 and the eccentric ratio is 0.1 to 0.5. The flow patterns, velocity vectors and vorticities are obtained. The gap width between two eccentric cylinders changes in a circumferential direction and the maximum Taylor vortex is found in the region where the gap width becomes narrow in the flow direction and the vortex becomes weakest in the region where the gap width is expanding. When the Re number is slightly larger than the critical Re number, the Taylor vortex flow and parallel stable flow coexist in a fluid between two eccentric rotating cylinders.