\[ a y''(x) \wp (x;\text {g2},\text {g3})+b y'(x) \wp '(x;\text {g2},\text {g3})+y(x) \left (c \left (6 \wp (x;\text {g2},\text {g3})^2-\frac {\text {g2}}{2}\right )+d\right )+y^{(4)}(x)=0 \] ✗ Mathematica : cpu = 0.0245882 (sec), leaf count = 0
DSolve[(d + c*(-1/2*g2 + 6*WeierstrassP[x, {g2, g3}]^2))*y[x] + b*WeierstrassPPrime[x, {g2, g3}]*Derivative[1][y][x] + a*WeierstrassP[x, {g2, g3}]*Derivative[2][y][x] + Derivative[4][y][x] == 0,y[x],x]
, could not solve
DSolve[(d + c*(-1/2*g2 + 6*WeierstrassP[x, {g2, g3}]^2))*y[x] + b*WeierstrassPPrime[x, {g2, g3}]*Derivative[1][y][x] + a*WeierstrassP[x, {g2, g3}]*Derivative[2][y][x] + Derivative[4][y][x] == 0, y[x], x]
✗ Maple : cpu = 0. (sec), leaf count = 0
dsolve(diff(diff(diff(diff(y(x),x),x),x),x)+a*WeierstrassP(x,g2,g3)*diff(diff(y(x),x),x)+b*WeierstrassPPrime(x,g2,g3)*diff(y(x),x)+(c*(6*WeierstrassP(x,g2,g3)^2-1/2*g2)+d)*y(x)=0,y(x))
, result contains DESol or ODESolStruc
\[y \left (x \right ) = \mathit {DESol}\left (\left \{\frac {d^{4}}{d x^{4}}\textit {\_Y} \left (x \right )+a \WeierstrassP \left (x , \mathit {g2} , \mathit {g3}\right ) \left (\frac {d^{2}}{d x^{2}}\textit {\_Y} \left (x \right )\right )+b \mathcal {P}^{\prime }\left (x ;\mathit {g2} ,\mathit {g3} \right ) \left (\frac {d}{d x}\textit {\_Y} \left (x \right )\right )+\left (6 c \WeierstrassP \left (x , \mathit {g2} , \mathit {g3}\right )^{2}-\frac {c \mathit {g2}}{2}+d \right ) \textit {\_Y} \left (x \right )\right \}, \left \{\textit {\_Y} \left (x \right )\right \}\right )\]