Added Jan 16, 2020.
Problem Chapter 9.2.3.1, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y,z)\) \[ w_x + (a_1 \sqrt x+ a_0) w_y + (b_1 \sqrt x+b_0) w_z = c w + s_1 \sqrt x+s_0 \]
Mathematica ✓
ClearAll["Global`*"]; pde = D[w[x,y,z],x]+(a1*Sqrt[x]+a0)*D[w[x,y,z],y]+(b1*Sqrt[x]+b0)*D[w[x,y,z],z]==c*w[x,y,z]+ s1*Sqrt[x]+s0; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x,y,z}], 60*10]];
\[\left \{\left \{w(x,y,z)\to e^{c x} c_1\left (-\text {a0} x-\frac {2}{3} \text {a1} x^{3/2}+y,-\text {b0} x-\frac {2}{3} \text {b1} x^{3/2}+z\right )+\frac {\sqrt {\pi } \text {s1} e^{c x} \text {erf}\left (\sqrt {c} \sqrt {x}\right )}{2 c^{3/2}}-\frac {\text {s0}+\text {s1} \sqrt {x}}{c}\right \}\right \}\]
Maple ✓
restart; local gamma; pde := diff(w(x,y,z),x)+ (a__1*sqrt(x)+a__0)*diff(w(x,y,z),y)+ (b__1*sqrt(x)+b__0)*diff(w(x,y,z),z)=c*w(x,y,z)+ s__1*sqrt(x)+s__0; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
\[w \left (x , y , z\right ) = \frac {\left (2 \sqrt {\frac {c}{\pi }}\, c \textit {\_F1} \left (-\frac {2 a_{1} x^{\frac {3}{2}}}{3}-a_{0} x +y , -\frac {2 b_{1} x^{\frac {3}{2}}}{3}-b_{0} x +z \right )+s_{1} \erf \left (\sqrt {c}\, \sqrt {x}\right )\right ) {\mathrm e}^{c x}-2 \sqrt {\frac {c}{\pi }}\, \left (s_{1} \sqrt {x}+s_{0} \right )}{2 \sqrt {\frac {c}{\pi }}\, c}\]
____________________________________________________________________________________
Added Jan 16, 2020.
Problem Chapter 9.2.3.2, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y,z)\) \[ w_x + (b_1 x^2+ b_0) w_y + (c_1 y^3+c_0) w_z = a w + s_1 x^3+s_0 \]
Mathematica ✓
ClearAll["Global`*"]; pde = D[w[x,y,z],x]+(b1*x^2+b0)*D[w[x,y,z],y]+(c1*y^3+c0)*D[w[x,y,z],z]==a*w[x,y,z]+ s1*x^3+s0; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x,y,z}], 60*10]];
\[\left \{\left \{w(x,y,z)\to -\frac {a^4 \left (-e^{a x}\right ) c_1\left (-\text {b0} x-\frac {\text {b1} x^3}{3}+y,\text {b0}^2 \left (\frac {19}{60} \text {b1} \text {c1} x^6-\text {c1} x^3 y\right )+\frac {1}{4} \text {b0}^3 \text {c1} x^4+\frac {3}{280} \text {b0} \text {c1} x^2 \left (13 \text {b1}^2 x^6-84 \text {b1} x^3 y+140 y^2\right )-\frac {3}{14} \text {b1}^2 \text {c1} x^7 y+\frac {3}{140} \text {b1}^3 \text {c1} x^{10}+\frac {3}{4} \text {b1} \text {c1} x^4 y^2-\text {c0} x-\text {c1} x y^3+z\right )+a^3 \left (\text {s0}+\text {s1} x^3\right )+3 a^2 \text {s1} x^2+6 a \text {s1} x+6 \text {s1}}{a^4}\right \}\right \}\]
Maple ✓
restart; local gamma; pde := diff(w(x,y,z),x)+ (b__1*x^2+b__0)*diff(w(x,y,z),y)+ (c__1*y^3+c__0)*diff(w(x,y,z),z)=a*w(x,y,z)+ s__1*x^3+s__0; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
\[w \left (x , y , z\right ) = \frac {a^{4} \textit {\_F1} \left (-\frac {1}{3} b_{1} x^{3}-b_{0} x +y , -c_{0} x +\frac {3 \left (b_{1}^{3} x^{9}+\frac {13 b_{0} b_{1}^{2} x^{7}}{2}-10 b_{1}^{2} x^{6} y +\frac {133 b_{0}^{2} b_{1} x^{5}}{9}-42 b_{0} b_{1} x^{4} y -\frac {140 b_{0}^{2} x^{2} y}{3}+70 b_{0} x \,y^{2}+\left (\frac {35 b_{0}^{3}}{3}+35 b_{1} y^{2}\right ) x^{3}-\frac {140 y^{3}}{3}\right ) c_{1} x}{140}+z \right ) {\mathrm e}^{a x}-3 a^{2} s_{1} x^{2}+\left (-s_{1} x^{3}-s_{0} \right ) a^{3}-6 a s_{1} x -6 s_{1}}{a^{4}}\]
____________________________________________________________________________________
Added Jan 16, 2020.
Problem Chapter 9.2.3.3, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y,z)\) \[ w_x + (a y+ k x^3) w_y + (b z+n x^3) w_z = c w + s x^2 \]
Mathematica ✓
ClearAll["Global`*"]; pde = D[w[x,y,z],x]+(a*y+k*x^3)*D[w[x,y,z],y]+(b*z+n*x^3)*D[w[x,y,z],z]==c*w[x,y,z]+ s*x^2; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x,y,z}], 60*10]];
\[\left \{\left \{w(x,y,z)\to -\frac {s \left (c^2 x^2+2 c x+2\right )}{c^3}+e^{c x} c_1\left (\frac {e^{-a x} \left (k \left (a^3 x^3+3 a^2 x^2+6 a x+6\right )+a^4 y\right )}{a^4},\frac {e^{-b x} \left (n \left (b^3 x^3+3 b^2 x^2+6 b x+6\right )+b^4 z\right )}{b^4}\right )\right \}\right \}\]
Maple ✓
restart; local gamma; pde := diff(w(x,y,z),x)+ (a*y+k*x^3)*diff(w(x,y,z),y)+ (b*z+n*x^3)*diff(w(x,y,z),z)=c*w(x,y,z)+ s*x^2; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
\[w \left (x , y , z\right ) = \frac {c^{3} \textit {\_F1} \left (\frac {\left (k \,x^{3} a^{3}+y \,a^{4}+3 k \,x^{2} a^{2}+6 k a x +6 k \right ) {\mathrm e}^{-a x}}{a^{4}}, \frac {\left (n \,x^{3} b^{3}+z \,b^{4}+3 n \,x^{2} b^{2}+6 b n x +6 n \right ) {\mathrm e}^{-b x}}{b^{4}}\right ) {\mathrm e}^{c x}-\left (x^{2} c^{2}+2 c x +2\right ) s}{c^{3}}\]
____________________________________________________________________________________
Added Jan 16, 2020.
Problem Chapter 9.2.3.4, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y,z)\) \[ w_x + (a_1 x y+ a_2 x^3) w_y + (b_1 y z+b_2 y^3) w_z = (c_1 z+c_2 y) w + s_1 x^2 y+ s_2 x z^2 \]
Mathematica ✗
ClearAll["Global`*"]; pde = D[w[x,y,z],x]+(a1*x*y+a2*x^3)*D[w[x,y,z],y]+(b1*y*z+b2*y^3)*D[w[x,y,z],z]==(c1*z+c2*y)*w[x,y,z]+ s1*x^2*y+s2*x*z^2; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x,y,z}], 60*10]];
$Aborted
Maple ✓
restart; local gamma; pde := diff(w(x,y,z),x)+ (a__1*x*y+a__2*x^3)*diff(w(x,y,z),y)+ (b__1*y*z+b__2*y^3)*diff(w(x,y,z),z)=(c__1*z+c__2*y)*w(x,y,z)+ s__1*x^2*y+s__2*x*z^2; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
\[\text {Expression too large to display}\]
____________________________________________________________________________________
Added Jan 16, 2020.
Problem Chapter 9.2.3.5, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y,z)\) \[ a x^3 w_x + b y^3 w_y + c z^3 w_z = x w + k x+ s \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*D[w[x,y,z],x]+b*y^3*D[w[x,y,z],y]+c*z^3*D[w[x,y,z],z]==x*w[x,y,z]+ k*x+s; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x,y,z}], 60*10]];
\[\left \{\left \{w(x,y,z)\to e^{\frac {x^2}{2 a}} c_1\left (-\frac {b x}{a}-\frac {1}{2 y^2},-\frac {c x}{a}-\frac {1}{2 z^2}\right )+\frac {\sqrt {\frac {\pi }{2}} s e^{\frac {x^2}{2 a}} \text {erf}\left (\frac {x}{\sqrt {2} \sqrt {a}}\right )}{\sqrt {a}}-k\right \}\right \}\]
Maple ✓
restart; local gamma; pde := a*diff(w(x,y,z),x)+ b*y^3*diff(w(x,y,z),y)+ c*z^3*diff(w(x,y,z),z)=x*w(x,y,z)+ k*x+s; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
\[w \left (x , y , z\right ) = \frac {\sqrt {2}\, s \erf \left (\frac {\sqrt {2}\, x}{2 \sqrt {a}}\right ) {\mathrm e}^{\frac {x^{2}}{2 a}}+2 \sqrt {\frac {a}{\pi }}\, \textit {\_F1} \left (\frac {2 b x \,y^{2}+a}{a \,y^{2}}, \frac {2 c x \,z^{2}+a}{a \,z^{2}}\right ) {\mathrm e}^{\frac {x^{2}}{2 a}}-2 \sqrt {\frac {a}{\pi }}\, k}{2 \sqrt {\frac {a}{\pi }}}\]
____________________________________________________________________________________