I come across a confusing point.
Considering the ODE for the function $y(x)$ $$\frac{1}{a^2} y'' + axy = 0$$ I consider the change of variable $ u = ax$ and, using the chain rule for a linear change of variables $$ \frac{\mathrm{d}^2 y}{\mathrm{d}x^2} = \frac{\mathrm{d}^2 y}{\mathrm{d}u^2} \left( \frac{\mathrm{d} u}{\mathrm{d}x} \right)^2 $$
I get the ODE
$$ y'' + uy = 0$$
If I on the other hand define $u = -ax$, I get the ODE $$ y'' - uy = 0$$
The solutions of the equations should be then identical upon changing the sign of the dependent variable. If I check (see http://mathworld.wolfram.com/AiryDifferentialEquation.html). I find the solution of the first equation is given by
$$ y(x) = \frac 1 3 \sqrt{x} \left[A I_{-1/3} \left(\frac 2 3 k x^{3/2}\right) - B I_{1/3} \left(\frac 2 3 k x^{3/2}\right) \right]$$ where $I$ stands for the modified Bessel function of the first kind, while the solution for the second is given by
$$y(x) = \frac 1 3 \sqrt{x} \left[A J_{-1/3} \left(\frac 2 3 k x^{3/2}\right) - B J_{1/3} \left(\frac 2 3 k x^{3/2}\right) \right]$$ where $J$ stands for the modified Bessel function of the first kind. I know that $$I_\beta(x) = i^{-\beta} J_\beta(ix),$$ but that does not imply
$$ I_\beta(x) = J_\beta(-x) $$ So how come by two different changes of variables I end up with different solutions? Thanks a lot