Fixing bugs in jacobian in f_discrete calculation. Why had Karan rotated body angular rates to body frame for the dq/dt calculation?
Hurtado's book suggests to use body frame angular rates for that calculation (A Matrix, chapter 5, euler parameter kinematics)