plot arctan((a+tan(x/2))/(sqrt(1-a^2))) for a={0,1/3,1/2,0.9} from x=0 to 2pi