plot pi + 2i -4arcsin((1/1.333)*sin(i)) between 0 and pi/2