Computes derivative of unsteady aerodynamic force with respect to changes in circulation.

Note: the function also checks that the first derivative of the circulation at the linearisation point is null. If not, a further contribution to the added mass, depending on the changes in panel area and normal, arises and needs to be implemented.