The results for the local-to-unity model and the "constructive" results (ie the augmented HPD set) in the Imbens-Manski problem were generated by the fortran code in "LTU.f90" and "IMconst.f90", respectively. The code uses service routines from the other .f90 files. Compilations requires the IMSL Numerical Libraries (we used version 7.01), and we set the default real type to double precision.