The local dependence function, introduced by Holland and Wang (1987) and studied by Wang (1993) as a continuous version of the local cross-ratio, describes the local relation between two random variables. Three explicit numerical algorithms are proposed to approximate bivariate densities given the marginal densities and the local dependence function. This approach is suited for simulation purposes, to provide illustrative examples of densities with given marginals, and for estimation of model parameters. The technique involves non-classical integral equation theory. The accuracy of the approximations is investigated.