6.8.23 8.1

6.8.23.1 [1891] Problem 1

problem number 1891

Added December 1, 2019.

Problem Chapter 8.8.1.1, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.

Solve for \(w(x,y,z)\) \[ w_x + f(x) w_y + g(x) w_z = \left ( h_2(x) y+h_1(x) z+h_0(x) \right ) w \]

Mathematica ✓

ClearAll["Global`*"]; 
pde =  D[w[x,y,z],x]+f[x]*D[w[x,y,z],y]+g[x]*D[w[x,y,z],z]==(h2[x]*y+h1[x]*z+h0[x])*w[x,y,z]; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x,y,z}], 60*10]];
 

\[\left \{\left \{w(x,y,z)\to c_1\left (y-\int _1^xf(K[1])dK[1],z-\int _1^xg(K[2])dK[2]\right ) \exp \left (\int _1^x\left (\text {h0}(K[3])+\text {h2}(K[3]) \left (y-\int _1^xf(K[1])dK[1]+\int _1^{K[3]}f(K[1])dK[1]\right )+\text {h1}(K[3]) \left (z-\int _1^xg(K[2])dK[2]+\int _1^{K[3]}g(K[2])dK[2]\right )\right )dK[3]\right )\right \}\right \}\]

Maple ✓

restart; 
local gamma; 
pde :=  diff(w(x,y,z),x)+ f(x)*diff(w(x,y,z),y)+ g(x)*diff(w(x,y,z),z)= (h2(x)*y+h1(x)*z+h0(x))*w(x,y,z); 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
 

\[w \left ( x,y,z \right ) ={\it \_F1} \left ( -\int \!f \left ( x \right ) \,{\rm d}x+y,-\int \!g \left ( x \right ) \,{\rm d}x+z \right ) {{\rm e}^{\int ^{x}\!{\it h2} \left ( {\it \_f} \right ) \int \!f \left ( {\it \_f} \right ) \,{\rm d}{\it \_f}+{\it h2} \left ( {\it \_f} \right ) \left ( -\int \!f \left ( x \right ) \,{\rm d}x+y \right ) +{\it h1} \left ( {\it \_f} \right ) \int \!g \left ( {\it \_f} \right ) \,{\rm d}{\it \_f}+{\it h1} \left ( {\it \_f} \right ) \left ( -\int \!g \left ( x \right ) \,{\rm d}x+z \right ) +{\it h0} \left ( {\it \_f} \right ) {d{\it \_f}}}}\]

____________________________________________________________________________________

6.8.23.2 [1892] Problem 2

problem number 1892

Added December 1, 2019.

Problem Chapter 8.8.1.2, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.

Solve for \(w(x,y,z)\) \[ w_x + f(x)(y+a) w_y + g(x)(z+b) w_z = h(x) w \]

Mathematica ✓

ClearAll["Global`*"]; 
pde =  D[w[x,y,z],x]+f[x]*(y+a)*D[w[x,y,z],y]+g[x]*(z+b)*D[w[x,y,z],z]==h[x]*w[x,y,z]; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x,y,z}], 60*10]];
 

\[\left \{\left \{w(x,y,z)\to \exp \left (\int _1^xh(K[5])dK[5]\right ) c_1\left (y \exp \left (-\int _1^xf(K[1])dK[1]\right )-\int _1^xa \exp \left (-\int _1^{K[2]}f(K[1])dK[1]\right ) f(K[2])dK[2],z \exp \left (-\int _1^xg(K[3])dK[3]\right )-\int _1^xb \exp \left (-\int _1^{K[4]}g(K[3])dK[3]\right ) g(K[4])dK[4]\right )\right \}\right \}\]

Maple ✓

restart; 
local gamma; 
pde :=  diff(w(x,y,z),x)+ f(x)*(y+a)*diff(w(x,y,z),y)+ g(x)*(z+b)*diff(w(x,y,z),z)= h(x)*w(x,y,z); 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
 

\[w \left ( x,y,z \right ) ={\it \_F1} \left ( \left ( y+a \right ) {{\rm e}^{-\int \!f \left ( x \right ) \,{\rm d}x}}, \left ( z+b \right ) {{\rm e}^{-\int \!g \left ( x \right ) \,{\rm d}x}} \right ) {{\rm e}^{\int \!h \left ( x \right ) \,{\rm d}x}}\]

____________________________________________________________________________________

6.8.23.3 [1893] Problem 3

problem number 1893

Added December 1, 2019.

Problem Chapter 8.8.1.3, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.

Solve for \(w(x,y,z)\) \[ w_x + \left ( a y+f(x) \right ) w_y + \left ( b z+g(x) \right ) w_z = h(x) w \]

Mathematica ✓

ClearAll["Global`*"]; 
pde =  D[w[x,y,z],x]+(a*y+f[x])*D[w[x,y,z],y]+(b*z+g[x])*D[w[x,y,z],z]==h[x]*w[x,y,z]; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x,y,z}], 60*10]];
 

\[\left \{\left \{w(x,y,z)\to \exp \left (\int _1^xh(K[3])dK[3]\right ) c_1\left (y e^{-a x}-\int _1^xe^{-a K[1]} f(K[1])dK[1],z e^{-b x}-\int _1^xe^{-b K[2]} g(K[2])dK[2]\right )\right \}\right \}\]

Maple ✓

restart; 
local gamma; 
pde :=  diff(w(x,y,z),x)+ (a*y+f(x))*diff(w(x,y,z),y)+ (b*z+g(x))*diff(w(x,y,z),z)= h(x)*w(x,y,z); 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
 

\[w \left ( x,y,z \right ) ={\it \_F1} \left ( -\int \!f \left ( x \right ) {{\rm e}^{-ax}}\,{\rm d}x+y{{\rm e}^{-ax}},-\int \!g \left ( x \right ) {{\rm e}^{-bx}}\,{\rm d}x+z{{\rm e}^{-bx}} \right ) {{\rm e}^{\int \!h \left ( x \right ) \,{\rm d}x}}\]

____________________________________________________________________________________

6.8.23.4 [1894] Problem 4

problem number 1894

Added December 1, 2019.

Problem Chapter 8.8.1.4, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.

Solve for \(w(x,y,z)\) \[ w_x + \left ( f_1(x) y+ f_2(x) \right ) w_y + \left ( g_1(x) y+ g_2(x) \right ) w_z = \left ( h_2(x) y+h_1(x) z+h_0(x) \right ) w \]

Mathematica ✓

ClearAll["Global`*"]; 
pde =  D[w[x,y,z],x]+(f1[x]*y+f2[x])*D[w[x,y,z],y]+(g1[x]*y+g2[x])*D[w[x,y,z],z]==(h2[x]*y+h1[x]*z+h0[x])*w[x,y,z]; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x,y,z}], 60*10]];
 

\[\left \{\left \{w(x,y,z)\to c_1\left (y \exp \left (-\int _1^x\text {f1}(K[1])dK[1]\right )-\int _1^x\exp \left (-\int _1^{K[3]}\text {f1}(K[1])dK[1]\right ) \text {f2}(K[3])dK[3],-\int _1^x\left (\text {g2}(K[4])-\exp \left (-\int _1^{K[4]}\text {f1}(K[1])dK[1]\right ) \text {f2}(K[4]) \int _1^{K[4]}\exp \left (\int _1^{K[2]}\text {f1}(K[1])dK[1]\right ) \text {g1}(K[2])dK[2]\right )dK[4]-y \exp \left (-\int _1^x\text {f1}(K[1])dK[1]\right ) \int _1^x\exp \left (\int _1^{K[2]}\text {f1}(K[1])dK[1]\right ) \text {g1}(K[2])dK[2]+z\right ) \exp \left (\int _1^x\exp \left (-\int _1^x\text {f1}(K[1])dK[1]\right ) \left (\exp \left (\int _1^x\text {f1}(K[1])dK[1]\right ) \text {h0}(K[5])+\exp \left (\int _1^{K[5]}\text {f1}(K[1])dK[1]\right ) \text {h2}(K[5]) \left (y-\exp \left (\int _1^x\text {f1}(K[1])dK[1]\right ) \int _1^x\exp \left (-\int _1^{K[3]}\text {f1}(K[1])dK[1]\right ) \text {f2}(K[3])dK[3]+\exp \left (\int _1^x\text {f1}(K[1])dK[1]\right ) \int _1^{K[5]}\exp \left (-\int _1^{K[3]}\text {f1}(K[1])dK[1]\right ) \text {f2}(K[3])dK[3]\right )+\text {h1}(K[5]) \left (-y \int _1^x\exp \left (\int _1^{K[2]}\text {f1}(K[1])dK[1]\right ) \text {g1}(K[2])dK[2]+\left (y-\exp \left (\int _1^x\text {f1}(K[1])dK[1]\right ) \int _1^x\exp \left (-\int _1^{K[3]}\text {f1}(K[1])dK[1]\right ) \text {f2}(K[3])dK[3]+\exp \left (\int _1^x\text {f1}(K[1])dK[1]\right ) \int _1^{K[5]}\exp \left (-\int _1^{K[3]}\text {f1}(K[1])dK[1]\right ) \text {f2}(K[3])dK[3]\right ) \int _1^{K[5]}\exp \left (\int _1^{K[2]}\text {f1}(K[1])dK[1]\right ) \text {g1}(K[2])dK[2]+\exp \left (\int _1^x\text {f1}(K[1])dK[1]\right ) \left (z-\int _1^x\left (\text {g2}(K[4])-\exp \left (-\int _1^{K[4]}\text {f1}(K[1])dK[1]\right ) \text {f2}(K[4]) \int _1^{K[4]}\exp \left (\int _1^{K[2]}\text {f1}(K[1])dK[1]\right ) \text {g1}(K[2])dK[2]\right )dK[4]+\int _1^{K[5]}\left (\text {g2}(K[4])-\exp \left (-\int _1^{K[4]}\text {f1}(K[1])dK[1]\right ) \text {f2}(K[4]) \int _1^{K[4]}\exp \left (\int _1^{K[2]}\text {f1}(K[1])dK[1]\right ) \text {g1}(K[2])dK[2]\right )dK[4]\right )\right )\right )dK[5]\right )\right \}\right \}\]

Maple ✓

restart; 
local gamma; 
pde :=  diff(w(x,y,z),x)+ (f1(x)*y+f2(x))*diff(w(x,y,z),y)+ (g1(x)*y+g2(x))*diff(w(x,y,z),z)= (h2(x)*y+h1(x)*z+h0(x))*w(x,y,z); 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
 

\[w \left ( x,y,z \right ) ={\it \_F1} \left ( -\int \!{\it f2} \left ( x \right ) {{\rm e}^{-\int \!{\it f1} \left ( x \right ) \,{\rm d}x}}\,{\rm d}x+y{{\rm e}^{-\int \!{\it f1} \left ( x \right ) \,{\rm d}x}},-\int ^{x}\!{{\rm e}^{\int \!{\it f1} \left ( {\it \_f} \right ) \,{\rm d}{\it \_f}}}{{\rm e}^{-\int \!{\it f1} \left ( x \right ) \,{\rm d}x}}y{\it g1} \left ( {\it \_f} \right ) -{\it g1} \left ( {\it \_f} \right ) {{\rm e}^{\int \!{\it f1} \left ( {\it \_f} \right ) \,{\rm d}{\it \_f}}}\int \!{\it f2} \left ( x \right ) {{\rm e}^{-\int \!{\it f1} \left ( x \right ) \,{\rm d}x}}\,{\rm d}x+{\it g1} \left ( {\it \_f} \right ) {{\rm e}^{\int \!{\it f1} \left ( {\it \_f} \right ) \,{\rm d}{\it \_f}}}\int \!{\it f2} \left ( {\it \_f} \right ) {{\rm e}^{-\int \!{\it f1} \left ( {\it \_f} \right ) \,{\rm d}{\it \_f}}}\,{\rm d}{\it \_f}+{\it g2} \left ( {\it \_f} \right ) {d{\it \_f}}+z \right ) {{\rm e}^{\int ^{x}\!{\it h2} \left ( {\it \_g} \right ) y{{\rm e}^{\int \!{\it f1} \left ( {\it \_g} \right ) \,{\rm d}{\it \_g}}}{{\rm e}^{-\int \!{\it f1} \left ( x \right ) \,{\rm d}x}}-{\it h2} \left ( {\it \_g} \right ) {{\rm e}^{\int \!{\it f1} \left ( {\it \_g} \right ) \,{\rm d}{\it \_g}}}\int \!{\it f2} \left ( x \right ) {{\rm e}^{-\int \!{\it f1} \left ( x \right ) \,{\rm d}x}}\,{\rm d}x+{\it h2} \left ( {\it \_g} \right ) {{\rm e}^{\int \!{\it f1} \left ( {\it \_g} \right ) \,{\rm d}{\it \_g}}}\int \!{\it f2} \left ( {\it \_g} \right ) {{\rm e}^{-\int \!{\it f1} \left ( {\it \_g} \right ) \,{\rm d}{\it \_g}}}\,{\rm d}{\it \_g}+{\it h1} \left ( {\it \_g} \right ) z+{\it h1} \left ( {\it \_g} \right ) \int \!{{\rm e}^{\int \!{\it f1} \left ( {\it \_g} \right ) \,{\rm d}{\it \_g}}}{{\rm e}^{-\int \!{\it f1} \left ( x \right ) \,{\rm d}x}}y{\it g1} \left ( {\it \_g} \right ) -{\it g1} \left ( {\it \_g} \right ) {{\rm e}^{\int \!{\it f1} \left ( {\it \_g} \right ) \,{\rm d}{\it \_g}}}\int \!{\it f2} \left ( x \right ) {{\rm e}^{-\int \!{\it f1} \left ( x \right ) \,{\rm d}x}}\,{\rm d}x+{\it g1} \left ( {\it \_g} \right ) {{\rm e}^{\int \!{\it f1} \left ( {\it \_g} \right ) \,{\rm d}{\it \_g}}}\int \!{\it f2} \left ( {\it \_g} \right ) {{\rm e}^{-\int \!{\it f1} \left ( {\it \_g} \right ) \,{\rm d}{\it \_g}}}\,{\rm d}{\it \_g}+{\it g2} \left ( {\it \_g} \right ) \,{\rm d}{\it \_g}-{\it h1} \left ( {\it \_g} \right ) \int ^{x}\!{{\rm e}^{\int \!{\it f1} \left ( {\it \_f} \right ) \,{\rm d}{\it \_f}}}{{\rm e}^{-\int \!{\it f1} \left ( x \right ) \,{\rm d}x}}y{\it g1} \left ( {\it \_f} \right ) -{\it g1} \left ( {\it \_f} \right ) {{\rm e}^{\int \!{\it f1} \left ( {\it \_f} \right ) \,{\rm d}{\it \_f}}}\int \!{\it f2} \left ( x \right ) {{\rm e}^{-\int \!{\it f1} \left ( x \right ) \,{\rm d}x}}\,{\rm d}x+{\it g1} \left ( {\it \_f} \right ) {{\rm e}^{\int \!{\it f1} \left ( {\it \_f} \right ) \,{\rm d}{\it \_f}}}\int \!{\it f2} \left ( {\it \_f} \right ) {{\rm e}^{-\int \!{\it f1} \left ( {\it \_f} \right ) \,{\rm d}{\it \_f}}}\,{\rm d}{\it \_f}+{\it g2} \left ( {\it \_f} \right ) {d{\it \_f}}+{\it h0} \left ( {\it \_g} \right ) {d{\it \_g}}}}\]

____________________________________________________________________________________

6.8.23.5 [1895] Problem 5

problem number 1895

Added December 1, 2019.

Problem Chapter 8.8.1.5, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.

Solve for \(w(x,y,z)\) \[ w_x + \left ( f_1(x) y+ f_2(x) \right ) w_y + \left ( g_1(x) z+ g_2(x) \right ) w_z = \left ( h_2(x) y+h_1(x) z+h_0(x) \right ) w \]

Mathematica ✓

ClearAll["Global`*"]; 
pde =  D[w[x,y,z],x]+(f1[x]*y+f2[x])*D[w[x,y,z],y]+(g1[x]*z+g2[x])*D[w[x,y,z],z]==(h2[x]*y+h1[x]*z+h0[x])*w[x,y,z]; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x,y,z}], 60*10]];
 

\[\left \{\left \{w(x,y,z)\to c_1\left (y \exp \left (-\int _1^x\text {f1}(K[1])dK[1]\right )-\int _1^x\exp \left (-\int _1^{K[2]}\text {f1}(K[1])dK[1]\right ) \text {f2}(K[2])dK[2],z \exp \left (-\int _1^x\text {g1}(K[3])dK[3]\right )-\int _1^x\exp \left (-\int _1^{K[4]}\text {g1}(K[3])dK[3]\right ) \text {g2}(K[4])dK[4]\right ) \exp \left (\int _1^x\left (\text {h0}(K[5])+\exp \left (\int _1^{K[5]}\text {f1}(K[1])dK[1]\right ) \text {h2}(K[5]) \left (\exp \left (-\int _1^x\text {f1}(K[1])dK[1]\right ) y-\int _1^x\exp \left (-\int _1^{K[2]}\text {f1}(K[1])dK[1]\right ) \text {f2}(K[2])dK[2]+\int _1^{K[5]}\exp \left (-\int _1^{K[2]}\text {f1}(K[1])dK[1]\right ) \text {f2}(K[2])dK[2]\right )+\exp \left (\int _1^{K[5]}\text {g1}(K[3])dK[3]\right ) \text {h1}(K[5]) \left (\exp \left (-\int _1^x\text {g1}(K[3])dK[3]\right ) z-\int _1^x\exp \left (-\int _1^{K[4]}\text {g1}(K[3])dK[3]\right ) \text {g2}(K[4])dK[4]+\int _1^{K[5]}\exp \left (-\int _1^{K[4]}\text {g1}(K[3])dK[3]\right ) \text {g2}(K[4])dK[4]\right )\right )dK[5]\right )\right \}\right \}\]

Maple ✓

restart; 
local gamma; 
pde :=  diff(w(x,y,z),x)+ (f1(x)*y+f2(x))*diff(w(x,y,z),y)+ (g1(x)*z+g2(x))*diff(w(x,y,z),z)= (h2(x)*y+h1(x)*z+h0(x))*w(x,y,z); 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
 

\[w \left ( x,y,z \right ) ={\it \_F1} \left ( -\int \!{\it f2} \left ( x \right ) {{\rm e}^{-\int \!{\it f1} \left ( x \right ) \,{\rm d}x}}\,{\rm d}x+y{{\rm e}^{-\int \!{\it f1} \left ( x \right ) \,{\rm d}x}},-\int \!{\it g2} \left ( x \right ) {{\rm e}^{-\int \!{\it g1} \left ( x \right ) \,{\rm d}x}}\,{\rm d}x+z{{\rm e}^{-\int \!{\it g1} \left ( x \right ) \,{\rm d}x}} \right ) {{\rm e}^{\int ^{x}\!{\it h1} \left ( {\it \_h} \right ) {{\rm e}^{\int \!{\it g1} \left ( {\it \_h} \right ) \,{\rm d}{\it \_h}}}{{\rm e}^{-\int \!{\it g1} \left ( x \right ) \,{\rm d}x}}z+{\it h2} \left ( {\it \_h} \right ) y{{\rm e}^{\int \!{\it f1} \left ( {\it \_h} \right ) \,{\rm d}{\it \_h}}}{{\rm e}^{-\int \!{\it f1} \left ( x \right ) \,{\rm d}x}}-{\it h1} \left ( {\it \_h} \right ) {{\rm e}^{\int \!{\it g1} \left ( {\it \_h} \right ) \,{\rm d}{\it \_h}}}\int \!{\it g2} \left ( x \right ) {{\rm e}^{-\int \!{\it g1} \left ( x \right ) \,{\rm d}x}}\,{\rm d}x+{\it h1} \left ( {\it \_h} \right ) {{\rm e}^{\int \!{\it g1} \left ( {\it \_h} \right ) \,{\rm d}{\it \_h}}}\int \!{\it g2} \left ( {\it \_h} \right ) {{\rm e}^{-\int \!{\it g1} \left ( {\it \_h} \right ) \,{\rm d}{\it \_h}}}\,{\rm d}{\it \_h}-{\it h2} \left ( {\it \_h} \right ) {{\rm e}^{\int \!{\it f1} \left ( {\it \_h} \right ) \,{\rm d}{\it \_h}}}\int \!{\it f2} \left ( x \right ) {{\rm e}^{-\int \!{\it f1} \left ( x \right ) \,{\rm d}x}}\,{\rm d}x+{\it h2} \left ( {\it \_h} \right ) {{\rm e}^{\int \!{\it f1} \left ( {\it \_h} \right ) \,{\rm d}{\it \_h}}}\int \!{\it f2} \left ( {\it \_h} \right ) {{\rm e}^{-\int \!{\it f1} \left ( {\it \_h} \right ) \,{\rm d}{\it \_h}}}\,{\rm d}{\it \_h}+{\it h0} \left ( {\it \_h} \right ) {d{\it \_h}}}}\]

____________________________________________________________________________________

6.8.23.6 [1896] Problem 6

problem number 1896

Added December 1, 2019.

Problem Chapter 8.8.1.6, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.

Solve for \(w(x,y,z)\) \[ w_x + \left ( y^2-a^2+a \lambda \sinh (\lambda x)-a^2 \sinh ^2(\lambda x) \right ) w_y + f(x) \sinh (\gamma z) w_z = g(x) w \]

Mathematica ✓

ClearAll["Global`*"]; 
pde =  D[w[x,y,z],x]+(y^2-a^2+a*lambda*Sinh[lambda*x]-a^2*Sinh[lambda*x]^2)*D[w[x,y,z],y]+f[x]*Sinh[gamma*z]*D[w[x,y,z],z]==g[x]*w[x,y,z]; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x,y,z}], 60*10]];
 

\[\left \{\left \{w(x,y,z)\to \exp \left (\int _1^xg(K[3])dK[3]\right ) c_1\left (\frac {\log \left (\tanh \left (\frac {\gamma z}{2}\right )\right )}{\gamma }-\int _1^xf(K[2])dK[2],\frac {2 \lambda e^{\frac {a e^{-\lambda x} \left (e^{2 \lambda x}-1\right )}{\lambda }+\lambda x}}{a e^{2 \lambda x}+a-2 y e^{\lambda x}}-\int _1^{e^{\lambda x}}\frac {e^{\frac {a \left (K[1]^2-1\right )}{\lambda K[1]}}}{K[1]}dK[1]\right )\right \}\right \}\]

Maple ✓

restart; 
local gamma; 
pde :=  diff(w(x,y,z),x)+ (y^2-a^2+a*lambda*sinh(lambda*x)-a^2*sinh(lambda*x)^2)*diff(w(x,y,z),y)+ f(x)*sinh(gamma*z)*diff(w(x,y,z),z)= g(x)*w(x,y,z); 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
 

\[w \left ( x,y,z \right ) ={{\rm e}^{\int \!g \left ( x \right ) \,{\rm d}x}}{\it \_F1} \left ( -2\,{\sqrt {\sinh \left ( \lambda \,x \right ) +i} \left ( \left ( i \left ( \sinh \left ( \lambda \,x \right ) \right ) ^{2}-i+2\,\sinh \left ( \lambda \,x \right ) \right ) \left ( a\cosh \left ( \lambda \,x \right ) +y \right ) \HeunC \left ( {\frac {4\,ia}{\lambda }},-1/2,-1/2,{\frac {2\,ia}{\lambda }},-1/8\,{\frac {8\,ia-3\,\lambda }{\lambda }},1/2-i/2\sinh \left ( \lambda \,x \right ) \right ) - \left ( i\sinh \left ( \lambda \,x \right ) -1/2\, \left ( \sinh \left ( \lambda \,x \right ) \right ) ^{2}+1/2 \right ) \cosh \left ( \lambda \,x \right ) \lambda \,\HeunCPrime \left ( {\frac {4\,ia}{\lambda }},-1/2,-1/2,{\frac {2\,ia}{\lambda }},-1/8\,{\frac {8\,ia-3\,\lambda }{\lambda }},1/2-i/2\sinh \left ( \lambda \,x \right ) \right ) \right ) \left ( \left ( \left ( 2\,i \left ( \sinh \left ( \lambda \,x \right ) \right ) ^{3}a+ \left ( i\lambda +2\,a \right ) \left ( \sinh \left ( \lambda \,x \right ) \right ) ^{2}+ \left ( 2\,ia+2\,\lambda \right ) \sinh \left ( \lambda \,x \right ) -i\lambda +2\,a \right ) \cosh \left ( \lambda \,x \right ) +2\,y \left ( 1+i\sinh \left ( \lambda \,x \right ) \right ) \left ( \left ( \sinh \left ( \lambda \,x \right ) \right ) ^{2}+1 \right ) \right ) \HeunC \left ( {\frac {4\,ia}{\lambda }},1/2,-1/2,{\frac {2\,ia}{\lambda }},-1/8\,{\frac {8\,ia-3\,\lambda }{\lambda }},1/2-i/2\sinh \left ( \lambda \,x \right ) \right ) -\HeunCPrime \left ( {\frac {4\,ia}{\lambda }},1/2,-1/2,{\frac {2\,ia}{\lambda }},-1/8\,{\frac {8\,ia-3\,\lambda }{\lambda }},1/2-i/2\sinh \left ( \lambda \,x \right ) \right ) \left ( -\sinh \left ( \lambda \,x \right ) +i \right ) \cosh \left ( \lambda \,x \right ) \lambda \, \left ( \left ( \sinh \left ( \lambda \,x \right ) \right ) ^{2}+1 \right ) \right ) ^{-1}},{\frac {-\int \!f \left ( x \right ) \,{\rm d}x\gamma -2\,\arctanh \left ( {{\rm e}^{\gamma \,z}} \right ) }{\gamma }} \right ) \]

____________________________________________________________________________________

6.8.23.7 [1897] Problem 7

problem number 1897

Added December 1, 2019.

Problem Chapter 8.8.1.7, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.

Solve for \(w(x,y,z)\) \[ w_x + \left ( f_1(x) y+f_2(x) y^k \right ) w_y + \left ( g_1(x) z+g_2(x) z^m \right ) w_z = h(x) w \]

Mathematica ✓

ClearAll["Global`*"]; 
pde =  D[w[x,y,z],x]+( f1[x]*y+f2[x]*y^k)*D[w[x,y,z],y]+(g1[x]*z+g2[x]*z^m)*D[w[x,y,z],z]==h[x]*w[x,y,z]; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x,y,z}], 60*10]];
 

\[\left \{\left \{w(x,y,z)\to \exp \left (\int _1^xh(K[5])dK[5]\right ) c_1\left ((k-1) \int _1^x\exp \left ((k-1) \int _1^{K[2]}\text {f1}(K[1])dK[1]\right ) \text {f2}(K[2])dK[2]+y^{1-k} \exp \left ((k-1) \int _1^x\text {f1}(K[1])dK[1]\right ),(m-1) \int _1^x\exp \left ((m-1) \int _1^{K[4]}\text {g1}(K[3])dK[3]\right ) \text {g2}(K[4])dK[4]+z^{1-m} \exp \left ((m-1) \int _1^x\text {g1}(K[3])dK[3]\right )\right )\right \}\right \}\]

Maple ✓

restart; 
local gamma; 
pde :=  diff(w(x,y,z),x)+ ( f1(x)*y+f2(x)*y^k)*diff(w(x,y,z),y)+ (g1(x)*z+g2(x)*z^m)*diff(w(x,y,z),z)= h(x)*w(x,y,z); 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
 

\[w \left ( x,y,z \right ) ={\it \_F1} \left ( \left ( k-1 \right ) \int \!{{\rm e}^{ \left ( k-1 \right ) \int \!{\it f1} \left ( x \right ) \,{\rm d}x}}{\it f2} \left ( x \right ) \,{\rm d}x+{y}^{-k+1}{{\rm e}^{ \left ( k-1 \right ) \int \!{\it f1} \left ( x \right ) \,{\rm d}x}}, \left ( -1+m \right ) \int \!{{\rm e}^{ \left ( -1+m \right ) \int \!{\it g1} \left ( x \right ) \,{\rm d}x}}{\it g2} \left ( x \right ) \,{\rm d}x+{z}^{1-m}{{\rm e}^{ \left ( -1+m \right ) \int \!{\it g1} \left ( x \right ) \,{\rm d}x}} \right ) {{\rm e}^{\int \!h \left ( x \right ) \,{\rm d}x}}\]

____________________________________________________________________________________

6.8.23.8 [1898] Problem 8

problem number 1898

Added December 1, 2019.

Problem Chapter 8.8.1.8, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.

Solve for \(w(x,y,z)\) \[ w_x + \left ( f_1(x) y+f_2(x) y^k \right ) w_y + \left ( g_1(x) +g_2(x) e^{\lambda z} \right ) w_z = h(x) w \]

Mathematica ✗

ClearAll["Global`*"]; 
pde =  D[w[x,y,z],x]+( f1[x]*y+f2[x]*y^k)*D[w[x,y,z],y]+(g1[x]+g2[x]*Exp[lambda*z])*D[w[x,y,z],z]==h[x]*w[x,y,z]; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x,y,z}], 60*10]];
 

Failed

Maple ✓

restart; 
local gamma; 
pde :=  diff(w(x,y,z),x)+ ( f1(x)*y+f2(x)*y^k)*diff(w(x,y,z),y)+ (g1(x)+g2(x)*exp(lambda*z))*diff(w(x,y,z),z)= h(x)*w(x,y,z); 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
 

\[w \left ( x,y,z \right ) ={\it \_F1} \left ( \left ( k-1 \right ) \int \!{{\rm e}^{ \left ( k-1 \right ) \int \!{\it f1} \left ( x \right ) \,{\rm d}x}}{\it f2} \left ( x \right ) \,{\rm d}x+{y}^{-k+1}{{\rm e}^{ \left ( k-1 \right ) \int \!{\it f1} \left ( x \right ) \,{\rm d}x}},{\frac {-\lambda \,\int \!{\it g2} \left ( x \right ) {{\rm e}^{\lambda \,\int \!{\it g1} \left ( x \right ) \,{\rm d}x}}\,{\rm d}x-{{\rm e}^{\lambda \, \left ( \int \!{\it g1} \left ( x \right ) \,{\rm d}x-z \right ) }}}{\lambda }} \right ) {{\rm e}^{\int \!h \left ( x \right ) \,{\rm d}x}}\]

____________________________________________________________________________________

6.8.23.9 [1899] Problem 9

problem number 1899

Added December 1, 2019.

Problem Chapter 8.8.1.9, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.

Solve for \(w(x,y,z)\) \[ w_x + \left ( f_1(x)+f_2(x) e^{\lambda y} \right ) w_y + \left ( g_1(x) +g_2(x) e^{\beta z} \right ) w_z = h(x) w \]

Mathematica ✗

ClearAll["Global`*"]; 
pde =  D[w[x,y,z],x]+( f1[x]+f2[x]*Exp[lambda*y])*D[w[x,y,z],y]+(g1[x]+g2[x]*Exp[beta*z])*D[w[x,y,z],z]==h[x]*w[x,y,z]; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x,y,z}], 60*10]];
 

Failed

Maple ✓

restart; 
local gamma; 
pde :=  diff(w(x,y,z),x)+ ( f1(x)+f2(x)*exp(lambda*y))*diff(w(x,y,z),y)+ (g1(x)+g2(x)*exp(beta*z))*diff(w(x,y,z),z)= h(x)*w(x,y,z); 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
 

\[w \left ( x,y,z \right ) ={\it \_F1} \left ( {\frac {-\lambda \,\int \!{\it f2} \left ( x \right ) {{\rm e}^{\lambda \,\int \!{\it f1} \left ( x \right ) \,{\rm d}x}}\,{\rm d}x-{{\rm e}^{-\lambda \, \left ( y-\int \!{\it f1} \left ( x \right ) \,{\rm d}x \right ) }}}{\lambda }},{\frac {-\beta \,\int \!{\it g2} \left ( x \right ) {{\rm e}^{\beta \,\int \!{\it g1} \left ( x \right ) \,{\rm d}x}}\,{\rm d}x-{{\rm e}^{\beta \, \left ( \int \!{\it g1} \left ( x \right ) \,{\rm d}x-z \right ) }}}{\beta }} \right ) {{\rm e}^{\int \!h \left ( x \right ) \,{\rm d}x}}\]

____________________________________________________________________________________