df/dx = x/a*sqrt(4*c^2*x^2+1)