plot y = arctan(x/(1+x^2)) + pi for -4 < x < 4