This work is to develop the FENE-P model by using a high-order accuracy finite difference method. We performed an effective numerical simulation scheme with the upwind difference method, to solve the polymer additive solution in the channel flow. In constitutive equation, the 3rd-order accurate upwind difference scheme was applied and steady resolution was achieved