this should be

```python
kinetic_differentials.append(omega1 - alphad[i]) 
```

omega = alpha'