A procedure and a FORTRAN program are presented which enable numerical modelling of viscous flow around a rotating rigid ellipsoidal inclusion. Examples of modelling were visualized using the MATLAB Package.