plot 1/2*arcsin((1+tan(t/2)^4)/2) - (arctan(sqrt(cos(t)+sec(t)+1))-pi/4) from t=-pi/2 to pi/2