Linearize model

[1-Bart-1]

Add a function linearize(solver, body_aero) that calculates the jacobian of the aerodynamics from inputs (deformation) to outputs (aerodynamic coefficients).