expand ((r+(sqrt((cosy+siny)²-1) - (cosy+siny) )) (r - (sqrt((cosy+siny)²-1) +(cosy+siny)