\[ \boxed { h \left ( y \left ( x \right ) \right ) {\frac {{\rm d}^{2}}{{\rm d}{x}^{2}}}y \left ( x \right ) +a\mbox {D} \left ( h \right ) \left ( y \left ( x \right ) \right ) \left ( {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) \right ) ^{2}+j \left ( y \left ( x \right ) \right ) =0} \]
Mathematica: cpu = 13.253683 (sec), leaf count = 116 \[ \left \{\left \{y(x)\to \text {InverseFunction}\left [\int _1^{\text {$\#$1}} -\frac {e^{a K[2]}}{\sqrt {2 \int _1^{K[2]} -\frac {e^{2 a K[1]} j(K[1])}{h(K[1])} \, dK[1]+c_1}} \, dK[2]\& \right ]\left [c_2+x\right ]\right \},\left \{y(x)\to \text {InverseFunction}\left [\int _1^{\text {$\#$1}} \frac {e^{a K[3]}}{\sqrt {2 \int _1^{K[3]} -\frac {e^{2 a K[1]} j(K[1])}{h(K[1])} \, dK[1]+c_1}} \, dK[3]\& \right ]\left [c_2+x\right ]\right \}\right \} \]
Maple: cpu = 0.125 (sec), leaf count = 90 \[ \left \{ \int ^{y \left ( x \right ) }\!{\frac {1}{ \left ( h \left ( { \it \_b} \right ) \right ) ^{-a}}{\frac {1}{\sqrt {-2\,\int \!{\frac { \left ( \left ( h \left ( {\it \_b} \right ) \right ) ^{a} \right ) ^{2}j \left ( {\it \_b} \right ) }{h \left ( {\it \_b} \right ) }}\,{\rm d}{ \it \_b}+{\it \_C1}}}}}{d{\it \_b}}-x-{\it \_C2}=0,\int ^{y \left ( x \right ) }\!-{\frac {1}{ \left ( h \left ( {\it \_b} \right ) \right ) ^{ -a}}{\frac {1}{\sqrt {-2\,\int \!{\frac { \left ( \left ( h \left ( { \it \_b} \right ) \right ) ^{a} \right ) ^{2}j \left ( {\it \_b} \right ) }{h \left ( {\it \_b} \right ) }}\,{\rm d}{\it \_b}+{\it \_C1}}}}}{d{ \it \_b}}-x-{\it \_C2}=0 \right \} \]