I found rotations by modifying the addimpulsivetorque function. instead of modifying omega I modfied L(angular momentum). I used addimpulsivetorque function instead of addtorque