Numerical integration of the contravariant form of the two phase flow motion equations