I am trying to numerically integrate the following equation where omega and theta are (3X1) vectors and W(w) is 4x4 skew symmetric matrix
def skew_omega(w):