These angular velocity definitions are not consistent with your Body fixed rotations. These should be:

```python
A.set_ang_vel(I, omega[0] * I.z + delta[0] * (cos(alpha[0]) * I.x + sin(alpha[0[) * I.y))
```

Or you can change your orientations to Space fixed and leave these definitions as is.