ODE
\[ y'(x)=x^2+2 y(x)+3 \cosh (x) \] ODE Classification
[[_linear, `class A`]]
Book solution method
Linear ODE
Mathematica ✓
cpu = 0.194809 (sec), leaf count = 46
\[\left \{\left \{y(x)\to -\frac {1}{4} e^{-x} \left (e^x \left (2 x^2+2 x+1\right )+6 e^{2 x}+2\right )+c_1 e^{2 x}\right \}\right \}\]
Maple ✓
cpu = 0.078 (sec), leaf count = 80
\[\left [y \left (x \right ) = {\mathrm e}^{2 x} \textit {\_C1} -\frac {{\mathrm e}^{2 x} \left (-2 x^{2} \sinh \left (2 x \right )+2 x^{2} \cosh \left (2 x \right )-2 x \sinh \left (2 x \right )+2 x \cosh \left (2 x \right )+6 \cosh \left (x \right )-6 \sinh \left (x \right )-2 \sinh \left (3 x \right )+2 \cosh \left (3 x \right )-\sinh \left (2 x \right )+\cosh \left (2 x \right )\right )}{4}\right ]\] Mathematica raw input
DSolve[y'[x] == x^2 + 3*Cosh[x] + 2*y[x],y[x],x]
Mathematica raw output
{{y[x] -> -1/4*(2 + 6*E^(2*x) + E^x*(1 + 2*x + 2*x^2))/E^x + E^(2*x)*C[1]}}
Maple raw input
dsolve(diff(y(x),x) = x^2+3*cosh(x)+2*y(x), y(x))
Maple raw output
[y(x) = exp(2*x)*_C1-1/4*exp(2*x)*(-2*x^2*sinh(2*x)+2*x^2*cosh(2*x)-2*x*sinh(2*x
)+2*x*cosh(2*x)+6*cosh(x)-6*sinh(x)-2*sinh(3*x)+2*cosh(3*x)-sinh(2*x)+cosh(2*x))
]