とりあえずy座標がボケてた
( x + ( a/2 )( 1 + sinθ ) , y + ( a/2 )( 1 - 2cosθ + sinθ ) )

あと90度回転なら計算しなくても
手作業で確かめれば100になるに決まってるし
0〜90度でx座標は増加一方だろうに
なんで30度でそんな微増、45度でほぼ最大、60度で一旦減少なん?

もしかして:ラジアンと度数法を取り違えている可能性