y=(1+x)*arctan(1/(1-x^2))