DSolve[Derivative[2][y][x] == (x*y[x])/(1 + x^3) - ((-1 + x^3)*Derivative[1][y][x])/(x*(1 + x^3)),y[x],x]