A direct numerical method is presented which calculates the shear rates in concentric cylinder viscometers using an integration formula of interpolatory type. The shear rate is approximated using a piecewise cubic polynomial interpolating at four successive data points. The unequally spaced interpolation points are the shear stresses at either the inner or outer cylinder.