I've got a strange mismatch (see image below) between linear position and velocity when comparing results of the SPART FD and an identical model in Simulink. I'm trying to verify that the mass and Coriolis matrices that are computed are correct for the attached simple example. The intention is to use these same functions to compute the dynamic model for legged robot I'm developing, where I'll use the mass and corilos matrices in the control scheme.
The position and velocity results differ significantly; however, only when I give the model an initial body angular velocity. I've included my code below. I'd immensely appreciate any effort to help resolve these differences.

float3D_2x2R_sim.zip