I recently posted a similar call (in fact I think identical except theta instead of t for the animate range variable) and had an interesting error on the range. That was solved by forcing an evalf on the realcons Pi/2, (using the Greek Pi tablet click).
I am now interested in the vectorCaculus form. My main interest is to visualize a vector's transformation as it applies to robotic frames and transformations. Watching in 3D a rotation is really key than seeing a start and finished view on how the transforms work and then the math for 'exp' and Log of velocity transforms, screw, transpose, twist, wrench, etc.
At the moment I am trying to see a two rotation sequence in animate where two rotation matrices are applied to the starting vector, [P]. Using simple matrix multiplication of the rotations matrices inside the arrow function for "what I think is the operation's output" of a two column vector of position for the animation - But, did not achieve the final expected position.
I am wondering now if the order of operations is the problem and the frame-by-frame vector is not completing properly.
P.S., love discovering these tidbits, but often wonder the challenges of finding these when debugging problematic outputs. The matrix math for frame transformations is a challenge in and of itself and I am pre-planning a proc() for the 4x4 and other transforms (skew-symmetric). Hoping Maple has this already in a package that I've yet to discover.
All thoughts welcome.