Before revision:  We devise a numerical method to find  the accessory parameters , solve the equation and to compute the modulus of the quadrilateral.
After revision:  We devise a numerical method to find the accessory parameters . This method uses the Schwarz equation to determine conformal mapping of the unit disk onto a given circular quadrilateral.