I am thinking about computing the (complex-valued) *incomplete* integral for a specific value of *m* and its inverse (as in the code below in this notebook), including with good precision near the singular corners, where many of the existing approximations have some difficulty. Computing the complete integral (as in the papers above) is a slightly different problem (though I think a rational approximation can potentially be helpful there too). My notebook above is just the disk-to-square case; I thought about but didn’t ever get around to implementing a disk-to-arbitrary-rectangle code (and computing the coefficients in Javascript would take porting a bunch of other building blocks).