diff(1-exp(-1.56704*atan((x/0.940794)^2.19271)^5.13827),x) where x=3.1