Recreate the symbolic function, main file and differential equation function to compute/integrate the Euler Angles 321 sequence. The 313 case is posted already, start from that.