[[_3rd_order, _linear, _nonhomogeneous]]
Book solution method
TO DO
Mathematica ✓
cpu = 0.291977 (sec), leaf count = 1185
Maple ✓
cpu = 0.714 (sec), leaf count = 1277
DSolve[4*y[x] + 2*y'[x] + y''[x] + y'''[x] == Sin[2*x],y[x],x]
Mathematica raw output
{{y[x] -> E^(x*Root[4 + 2*#1 + #1^2 + #1^3 & , 1, 0])*C[1] + E^(x*Root[4 + 2*#1
+ #1^2 + #1^3 & , 2, 0])*C[2] + E^(x*Root[4 + 2*#1 + #1^2 + #1^3 & , 3, 0])*C[3]
+ (8*Cos[2*x]*(Root[4 + 2*#1 + #1^2 + #1^3 & , 2, 0]*Root[4 + 2*#1 + #1^2 + #1^
3 & , 3, 0]*(-252 + (-70 + 26*Root[4 + 2*#1 + #1^2 + #1^3 & , 1, 0] + Root[4 + 2
*#1 + #1^2 + #1^3 & , 1, 0]^2)*Root[4 + 2*#1 + #1^2 + #1^3 & , 3, 0]) + 2*(166 -
(64 + 118*Root[4 + 2*#1 + #1^2 + #1^3 & , 1, 0] + 3*Root[4 + 2*#1 + #1^2 + #1^3
& , 1, 0]^2)*Root[4 + 2*#1 + #1^2 + #1^3 & , 3, 0] + (-84 - 15*Root[4 + 2*#1 +
#1^2 + #1^3 & , 1, 0] + 5*Root[4 + 2*#1 + #1^2 + #1^3 & , 1, 0]^2)*Root[4 + 2*#1
+ #1^2 + #1^3 & , 3, 0]^2) + 2*Root[4 + 2*#1 + #1^2 + #1^3 & , 2, 0]^2*(-26 - 3
*Root[4 + 2*#1 + #1^2 + #1^3 & , 3, 0] + 18*Root[4 + 2*#1 + #1^2 + #1^3 & , 3, 0
]^2 + Root[4 + 2*#1 + #1^2 + #1^3 & , 1, 0]*(10 + 25*Root[4 + 2*#1 + #1^2 + #1^3
& , 3, 0] + 8*Root[4 + 2*#1 + #1^2 + #1^3 & , 3, 0]^2))) + (Root[4 + 2*#1 + #1^
2 + #1^3 & , 2, 0]*(128 + 24*Root[4 + 2*#1 + #1^2 + #1^3 & , 3, 0] - 5*(-24 - 50
*Root[4 + 2*#1 + #1^2 + #1^3 & , 1, 0] + 3*Root[4 + 2*#1 + #1^2 + #1^3 & , 1, 0]
^2)*Root[4 + 2*#1 + #1^2 + #1^3 & , 3, 0]^2) + 4*(24 + (116 - 30*Root[4 + 2*#1 +
#1^2 + #1^3 & , 1, 0] + 23*Root[4 + 2*#1 + #1^2 + #1^3 & , 1, 0]^2)*Root[4 + 2*
#1 + #1^2 + #1^3 & , 3, 0] + (36 + 6*Root[4 + 2*#1 + #1^2 + #1^3 & , 1, 0] + 11*
Root[4 + 2*#1 + #1^2 + #1^3 & , 1, 0]^2)*Root[4 + 2*#1 + #1^2 + #1^3 & , 3, 0]^2
) + Root[4 + 2*#1 + #1^2 + #1^3 & , 2, 0]^2*(Root[4 + 2*#1 + #1^2 + #1^3 & , 1,
0]*(-88 + 96*Root[4 + 2*#1 + #1^2 + #1^3 & , 3, 0] + 25*Root[4 + 2*#1 + #1^2 + #
1^3 & , 3, 0]^2) + 4*(52 + 25*Root[4 + 2*#1 + #1^2 + #1^3 & , 3, 0] + 51*Root[4
+ 2*#1 + #1^2 + #1^3 & , 3, 0]^2)))*Sin[2*x])/(Root[4 + 2*#1 + #1^2 + #1^3 & , 1
, 0]^2*(Root[4 + 2*#1 + #1^2 + #1^3 & , 1, 0] - Root[4 + 2*#1 + #1^2 + #1^3 & ,
2, 0])*Root[4 + 2*#1 + #1^2 + #1^3 & , 2, 0]*(4 + Root[4 + 2*#1 + #1^2 + #1^3 &
, 2, 0]^2)*(Root[4 + 2*#1 + #1^2 + #1^3 & , 1, 0] - Root[4 + 2*#1 + #1^2 + #1^3
& , 3, 0])*(Root[4 + 2*#1 + #1^2 + #1^3 & , 2, 0] - Root[4 + 2*#1 + #1^2 + #1^3
& , 3, 0])*((1 - 2*I) + Root[4 + 2*#1 + #1^2 + #1^3 & , 2, 0] + Root[4 + 2*#1 +
#1^2 + #1^3 & , 3, 0])*((1 + 2*I) + Root[4 + 2*#1 + #1^2 + #1^3 & , 2, 0] + Root
[4 + 2*#1 + #1^2 + #1^3 & , 3, 0])*(-2 + Root[4 + 2*#1 + #1^2 + #1^3 & , 2, 0]^2
+ 2*Root[4 + 2*#1 + #1^2 + #1^3 & , 1, 0]*Root[4 + 2*#1 + #1^2 + #1^3 & , 3, 0]
)*(4 + Root[4 + 2*#1 + #1^2 + #1^3 & , 3, 0]^2))}}
Maple raw input
dsolve(diff(diff(diff(y(x),x),x),x)+diff(diff(y(x),x),x)+2*diff(y(x),x)+4*y(x) = sin(2*x), y(x),'implicit')
Maple raw output
y(x) = (-207500*exp(-1/150*(3*249^(1/2)*(46+3*249^(1/2))^(2/3)-46*(46+3*249^(1/2
))^(2/3)-25*(46+3*249^(1/2))^(1/3)+50)*x)*cos(23/75*(46+3*3^(1/2)*83^(1/2))^(1/3
)*x*((3^(1/2)-9/46*83^(1/2))*(46+3*3^(1/2)*83^(1/2))^(1/3)-25/46*3^(1/2)))*((46+
3*3^(1/2)*83^(1/2))^(1/3)+(-3/25*3^(1/2)*83^(1/2)+117/50)*(46+3*3^(1/2)*83^(1/2)
)^(2/3)+27/2+(-138/625*3^(1/2)*83^(1/2)+4357/1250)*(46+3*3^(1/2)*83^(1/2))^(4/3)
)*Int(3/4150*83^(1/2)*((25/9*3^(1/2)*(46+3*3^(1/2)*83^(1/2))^(1/3)-46/9*(3^(1/2)
-9/46*83^(1/2))*(46+3*3^(1/2)*83^(1/2))^(2/3))*cos(23/75*(46+3*3^(1/2)*83^(1/2))
^(1/3)*x*((3^(1/2)-9/46*83^(1/2))*(46+3*3^(1/2)*83^(1/2))^(1/3)-25/46*3^(1/2)))+
((249^(1/2)-46/3)*(46+3*249^(1/2))^(2/3)-25/3*(46+3*249^(1/2))^(1/3))*sin(23/75*
(46+3*3^(1/2)*83^(1/2))^(1/3)*x*((3^(1/2)-9/46*83^(1/2))*(46+3*3^(1/2)*83^(1/2))
^(1/3)-25/46*3^(1/2))))*sin(2*x)*exp(1/150*(3*249^(1/2)*(46+3*249^(1/2))^(2/3)-4
6*(46+3*249^(1/2))^(2/3)-25*(46+3*249^(1/2))^(1/3)+50)*x),x)-207500*exp(-1/150*(
3*249^(1/2)*(46+3*249^(1/2))^(2/3)-46*(46+3*249^(1/2))^(2/3)-25*(46+3*249^(1/2))
^(1/3)+50)*x)*sin(23/75*(46+3*3^(1/2)*83^(1/2))^(1/3)*x*((3^(1/2)-9/46*83^(1/2))
*(46+3*3^(1/2)*83^(1/2))^(1/3)-25/46*3^(1/2)))*((46+3*3^(1/2)*83^(1/2))^(1/3)+(-
3/25*3^(1/2)*83^(1/2)+117/50)*(46+3*3^(1/2)*83^(1/2))^(2/3)+27/2+(-138/625*3^(1/
2)*83^(1/2)+4357/1250)*(46+3*3^(1/2)*83^(1/2))^(4/3))*Int(1/498*83^(1/2)*sin(2*x
)*exp(1/150*(3*249^(1/2)*(46+3*249^(1/2))^(2/3)-46*(46+3*249^(1/2))^(2/3)-25*(46
+3*249^(1/2))^(1/3)+50)*x)*(((-9/25*249^(1/2)+138/25)*(46+3*249^(1/2))^(2/3)+3*(
46+3*249^(1/2))^(1/3))*cos(23/75*(46+3*3^(1/2)*83^(1/2))^(1/3)*x*((3^(1/2)-9/46*
83^(1/2))*(46+3*3^(1/2)*83^(1/2))^(1/3)-25/46*3^(1/2)))+sin(23/75*(46+3*3^(1/2)*
83^(1/2))^(1/3)*x*((3^(1/2)-9/46*83^(1/2))*(46+3*3^(1/2)*83^(1/2))^(1/3)-25/46*3
^(1/2)))*(3^(1/2)*(46+3*3^(1/2)*83^(1/2))^(1/3)-46/25*(3^(1/2)-9/46*83^(1/2))*(4
6+3*3^(1/2)*83^(1/2))^(2/3))),x)+(-3750*83^(1/2)*3^(1/2)*(cos(2*x)-1/6*sin(2*x))
*(46+3*3^(1/2)*83^(1/2))^(1/3)+(6900*3^(1/2)*83^(1/2)*cos(2*x)-525*3^(1/2)*83^(1
/2)*sin(2*x)-112050*cos(2*x)+18675*sin(2*x))*(46+3*3^(1/2)*83^(1/2))^(2/3)-4357*
sin(2*x)*(3^(1/2)*83^(1/2)-68724/4357)*(46+3*3^(1/2)*83^(1/2))^(4/3))*exp(-1/75*
x*(3*3^(1/2)*83^(1/2)*(46+3*3^(1/2)*83^(1/2))^(2/3)-3*249^(1/2)*(46+3*249^(1/2))
^(2/3)-46*(46+3*3^(1/2)*83^(1/2))^(2/3)+46*(46+3*249^(1/2))^(2/3)-25*(46+3*3^(1/
2)*83^(1/2))^(1/3)+25*(46+3*249^(1/2))^(1/3)))+207500*(_C2*exp(-1/150*(3*249^(1/
2)*(46+3*249^(1/2))^(2/3)-46*(46+3*249^(1/2))^(2/3)-25*(46+3*249^(1/2))^(1/3)+50
)*x)*cos(1/150*(46*3^(1/2)-9*83^(1/2))*x*(46+3*3^(1/2)*83^(1/2))^(2/3)-1/6*(46+3
*3^(1/2)*83^(1/2))^(1/3)*3^(1/2)*x)+_C3*exp(-1/150*(3*249^(1/2)*(46+3*249^(1/2))
^(2/3)-46*(46+3*249^(1/2))^(2/3)-25*(46+3*249^(1/2))^(1/3)+50)*x)*sin(1/150*(46*
3^(1/2)-9*83^(1/2))*x*(46+3*3^(1/2)*83^(1/2))^(2/3)-1/6*(46+3*3^(1/2)*83^(1/2))^
(1/3)*3^(1/2)*x)+_C1*exp(1/75*(3*249^(1/2)*(46+3*249^(1/2))^(2/3)-46*(46+3*249^(
1/2))^(2/3)-25*(46+3*249^(1/2))^(1/3)-25)*x))*((46+3*3^(1/2)*83^(1/2))^(1/3)+(-3
/25*3^(1/2)*83^(1/2)+117/50)*(46+3*3^(1/2)*83^(1/2))^(2/3)+27/2+(-138/625*3^(1/2
)*83^(1/2)+4357/1250)*(46+3*3^(1/2)*83^(1/2))^(4/3)))/(207500*(46+3*3^(1/2)*83^(
1/2))^(1/3)+(-24900*3^(1/2)*83^(1/2)+485550)*(46+3*3^(1/2)*83^(1/2))^(2/3)+28012
50+(-45816*3^(1/2)*83^(1/2)+723262)*(46+3*3^(1/2)*83^(1/2))^(4/3))