I can’t understand how this is used. Eg his expression for r_t in the inverse at the start. How would you use this in a program to get the inverse? Guess just pointing out that you can iteratively find one coefficient at a time.
Also the differential Eqn looks integrable.