The fact that there is an analytical solution to the problem is a really nice achievement, but I have to wonder: solving this problem numerically has surely been within reach of computers for a long time, and should therefore have made actually manufacturing such lenses somewhat straightforward.
I assume that if numerically solving for a correcting shape was burdensome you would model an element using a first surface and a 'dummy' second surface that your software just treats as magically correcting. ... then when you are happy with your design, you go ahead and compute the actual shape for that second surface... or similar.
I would expect a closed form solution to more useful for meta analysis of the problem-- looking at its behaviors (and especially derivatives) in various external cases may suggest interesting and novel optical system designs. Like, "oh, foo changes cubically with with an infinite focal length that means if we could make lenses with cherry flavored unobtanium we could bounce a gravitation particle beam off the main deflector dish!".
What am i missing?