Document-Slva072 has a derivation of load regulation. Can you get the conclusion Req=1/G*(R1+R2)/R2 from the derivation? If so, then the assumption of equation (23) in the derivation of line regulation in the next section is somewhat problematic. Bringing the expression of Req into the final result (27) should be divided by 2.
In addition, I have doubts about the line regulation. If the saturation zone I is completely independent of Vds, then even if Vin changes, can Vo and I remain unchanged? However, this conclusion cannot be seen from (27). I think the line regulation is related to the unsatisfactory I-V curve. For example, the current in the saturation region is also related to Vds, which is often the case with actual devices.