add test_rot2.m to check formula for double rotation