A numerical solution has been developed for the steady state partial differential equations which discribe the flow of non-Newtonian fluids into geometrically and kinematically asymmetric calenders with viscous heating effects. The solution technuque combines the equation of motion at the isothermal condition. the asymmetric calenderign system is simplified by use of bicylindrical coordinates wiich establishes uniform size network for each variable. the initial condition is simply that the calender gap is filled completely with certain fluid whose initial temperature is uniformly constant. Viscous heating is shown to drastically change in the velocity and pressure field near the entrance andexit zone. The temperature distribution exhibits two local maxima and a minimum in the direction of flow, an dthe maximum temperature rise is about 30 ℃ near the roll surface. Numerically stability propblem of the steady state equations associated with high viscosity fluid flow is treated and solved by the use of the successive under relaxation method.