t=-pi/2:0.01:pi/2;
x=atan((2*(cos(t)-cos(pi/4)*sin(t)))./(2*(cos(t)-cos(pi/4)).*cos(t)+1));
plot(t,x)