Numerical solution of the partial differential equations that model the steady three‐dimensional flow and heat transfer of Carreau fluid between two stretchable rotatory disks