Effective viscosity of magnetic fluid is calculated numerically based on the diffusion equation for the spherical particles with permanent magnetic moment. The result is compared with that of the classical constitutive equation proposed by Shliomis. It is found that unlike Shliomis's equation, the viscosity remains monovalued as a function of the shear rate. A new constitutive equation which reproduces the numerical result is proposed.