Numerical solution of a linear-elliptic-partial-differential-equation for molecular orientation in shear flow by finite difference methods