This paper develops approximations of the Gram-Charlier type to the cumulative distribution function of the instrumental variables estimator on classical assumptions. In the special case where there are only two endogenous variables in the estimated equation, exact values of the cumulative distribution function are computed by numerical integration and compared with the approximations. Although the error in the approximation depends critically on the parameters of the stochastic model, the approximation is good for the special case even for small sample size over a wide range of values of the parameters.