Julio, you might want to check your dynamics books again. You have a obvious sign error in your formula.I can see that you, like Mr Kenb, have a great deal of confusion with Gyroscopic forces stuff, I am going to recommend to you to read my post about how to reproduce, http://www.physicsengine.com/forum/viewtopic.php?t=3746
You wrote:
But the last equality is wrong.D(R)/dt * I| * W| = W~ * R * I| * W| = - R * W|~ * I| * W|
Now:
W~ * R = R * W|~ * transpose(R) * R = R * W|~
So instead of:
You should have:T = -R * W|~ * I| * W| + I * A
T = R * W|~ * I| * W| + I * A
You might want to fix your post so that people don't try your solution and incorrectly judge your engine to be unstable.