plot θ - arcsin(tan(θ*pi/180))/pi*180 from θ=-30 to 30