Numerical modelling of ellipsoidal inclusions