6.3.6 3.1

6.3.6.1 [878] Problem 1
6.3.6.2 [879] Problem 2
6.3.6.3 [880] Problem 3
6.3.6.4 [881] Problem 4
6.3.6.5 [882] Problem 5
6.3.6.6 [883] Problem 6
6.3.6.7 [884] Problem 7
6.3.6.8 [885] Problem 8
6.3.6.9 [886] Problem 9
6.3.6.10 [887] Problem 10
6.3.6.11 [888] Problem 11

6.3.6.1 [878] Problem 1

problem number 878

Added Feb. 9, 2019.

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

Solve for \(w(x,y)\) \[ a w_x + b w_y = c e^{\lambda x} + d e^{\mu y} \]

Mathematica

ClearAll["Global`*"]; 
pde =  a*D[w[x, y], x] + b*D[w[x, y], y] == c*Exp[lambda*x] + d*Exp[mu*y]; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
 

\[\left \{\left \{w(x,y)\to c_1\left (y-\frac {b x}{a}\right )+\frac {c e^{\lambda x}}{a \lambda }+\frac {d e^{\mu y}}{b \mu }\right \}\right \}\]

Maple

restart; 
pde :=a*diff(w(x,y),x) +b*diff(w(x,y),y) =c*exp(lambda*x)+d*exp(mu*y); 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
 

\[w \left ( x,y \right ) ={\frac {1}{b\mu \,a\lambda } \left ( {\it \_F1} \left ( {\frac {ya-bx}{a}} \right ) b\mu \,a\lambda +{{\rm e}^{\lambda \,x}}cb\mu +d{{\rm e}^{\mu \,y}}a\lambda \right ) }\]

____________________________________________________________________________________

6.3.6.2 [879] Problem 2

problem number 879

Added Feb. 9, 2019.

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

Solve for \(w(x,y)\)

\[ a w_x + b w_y = c e^{\lambda x + \beta y} \]

Mathematica

ClearAll["Global`*"]; 
pde =  a*D[w[x, y], x] + b*D[w[x, y], y] == c*Exp[lambda*x + beta*y]; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
 

\[\left \{\left \{w(x,y)\to c_1\left (y-\frac {b x}{a}\right )+\frac {c e^{\beta y+\lambda x}}{a \lambda +b \beta }\right \}\right \}\]

Maple

restart; 
pde :=a*diff(w(x,y),x) +b*diff(w(x,y),y) =c*exp(lambda*x+beta*y); 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
 

\[w \left ( x,y \right ) ={\frac {c{{\rm e}^{\beta \,y+\lambda \,x}}}{a\lambda +b\beta }}+{\it \_F1} \left ( {\frac {ya-bx}{a}} \right ) \]

____________________________________________________________________________________

6.3.6.3 [880] Problem 3

problem number 880

Added Feb. 9, 2019.

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

Solve for \(w(x,y)\)

\[ a e^{\lambda x} w_x + b e^{\beta y} w_y = c \]

Mathematica

ClearAll["Global`*"]; 
pde =  a*Exp[lambda*x]*D[w[x, y], x] + b*Exp[beta*y]*D[w[x, y], y] == c; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
 

\[\left \{\left \{w(x,y)\to c_1\left (\frac {b e^{-\lambda x}}{a \lambda }-\frac {e^{-\beta y}}{\beta }\right )-\frac {c e^{-\lambda x}}{a \lambda }\right \}\right \}\]

Maple

restart; 
pde :=a*exp(lambda*x)*diff(w(x,y),x) +b*exp(beta*y)*diff(w(x,y),y) =c; 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
 

\[w \left ( x,y \right ) ={\frac {1}{a\lambda } \left ( {\it \_F1} \left ( -{\frac { \left ( a\lambda \,{{\rm e}^{\lambda \,x}}-{{\rm e}^{\beta \,y}}b\beta \right ) {{\rm e}^{-\beta \,y-\lambda \,x}}}{b\beta \,\lambda }} \right ) a\lambda -c{{\rm e}^{-\lambda \,x}} \right ) }\]

____________________________________________________________________________________

6.3.6.4 [881] Problem 4

problem number 881

Added Feb. 9, 2019.

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

Solve for \(w(x,y)\)

\[ a e^{\lambda y} w_x + b e^{\beta x} w_y = c \]

Mathematica

ClearAll["Global`*"]; 
pde =  a*Exp[lambda*y]*D[w[x, y], x] + b*Exp[beta*x]*D[w[x, y], y] == c; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
 

\[\left \{\left \{w(x,y)\to c_1\left (\frac {e^{\lambda y}}{\lambda }-\frac {b e^{\beta x}}{a \beta }\right )+\frac {c \left (\beta x-\log \left (\frac {a \beta e^{\lambda y}}{\lambda }\right )\right )}{a \beta e^{\lambda y}-b \lambda e^{\beta x}}\right \}\right \}\]

Maple

restart; 
pde :=a*exp(lambda*y)*diff(w(x,y),x) +b*exp(beta*x)*diff(w(x,y),y) =c; 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
 

\[w \left ( x,y \right ) ={\frac {1}{-{{\rm e}^{\beta \,x}}b\lambda +{{\rm e}^{y\lambda }}a\beta } \left ( \left ( -{{\rm e}^{\beta \,x}}b\lambda +{{\rm e}^{y\lambda }}a\beta \right ) {\it \_F1} \left ( {\frac {-{{\rm e}^{\beta \,x}}b\lambda +{{\rm e}^{y\lambda }}a\beta }{b\beta \,\lambda }} \right ) +c \left ( \ln \left ( {{\rm e}^{\beta \,x}} \right ) -\ln \left ( {\frac {{{\rm e}^{y\lambda }}a\beta }{\lambda \,b}} \right ) \right ) \right ) }\]

____________________________________________________________________________________

6.3.6.5 [882] Problem 5

problem number 882

Added Feb. 9, 2019.

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

Solve for \(w(x,y)\)

\[ a e^{\alpha x} w_x + b e^{\beta y} w_y = c e^{\gamma x-\beta y} \]

Mathematica

ClearAll["Global`*"]; 
pde =  a*Exp[alpha*x]*D[w[x, y], x] + b*Exp[beta*y]*D[w[x, y], y] == c*Exp[gamma*x - beta*y]; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
 

\[\left \{\left \{w(x,y)\to \frac {a^2 \left (2 \alpha ^2-3 \alpha \gamma +\gamma ^2\right ) c_1\left (\frac {b e^{-\alpha x}}{a \alpha }-\frac {e^{-\beta y}}{\beta }\right )+a c (\gamma -2 \alpha ) e^{-\alpha x-\beta y+\gamma x}+b \beta c e^{x (\gamma -2 \alpha )}}{a^2 \left (2 \alpha ^2-3 \alpha \gamma +\gamma ^2\right )}\right \}\right \}\]

Maple

restart; 
pde :=a*exp(alpha*x)*diff(w(x,y),x) +b*exp(beta*y)*diff(w(x,y),y) =c*exp(gamma*x-beta*y); 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
 

\[w \left ( x,y \right ) =128\,{\frac { \left ( \left ( -1/16+\alpha \right ) \left ( {{\rm e}^{\beta \,y}}b\beta -a\alpha \,{{\rm e}^{\alpha \,x}} \right ) {{\rm e}^{1/8\, \left ( -16\,\alpha +1 \right ) x-\beta \,y}}-1/2\,b{{\rm e}^{x/8-2\,\alpha \,x}}\beta \, \left ( -1/8+\alpha \right ) \right ) c}{{a}^{2}\alpha \, \left ( -1+8\,\alpha \right ) \left ( -1+16\,\alpha \right ) }}+{\it \_F1} \left ( -{\frac { \left ( a\alpha \,{{\rm e}^{\alpha \,x}}-{{\rm e}^{\beta \,y}}b\beta \right ) {{\rm e}^{-\alpha \,x-\beta \,y}}}{\alpha \,b\beta }} \right ) \]

____________________________________________________________________________________

6.3.6.6 [883] Problem 6

problem number 883

Added Feb. 9, 2019.

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

Solve for \(w(x,y)\)

\[ a e^{\alpha x} w_x + b e^{\beta y} w_y = c e^{\gamma x-2 \beta y} \]

Mathematica

ClearAll["Global`*"]; 
pde =  a*Exp[alpha*x]*D[w[x, y], x] + b*Exp[beta*y]*D[w[x, y], y] == c*Exp[gamma*x - 2*beta*y]; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
 

\[\left \{\left \{w(x,y)\to \frac {e^{-2 (\alpha x+\beta y)} \left (a^3 \left (6 \alpha ^3-11 \alpha ^2 \gamma +6 \alpha \gamma ^2-\gamma ^3\right ) e^{2 \alpha x+2 \beta y} c_1\left (\frac {b e^{-\alpha x}}{a \alpha }-\frac {e^{-\beta y}}{\beta }\right )-c e^{x (\gamma -\alpha )} \left (a^2 \left (6 \alpha ^2-5 \alpha \gamma +\gamma ^2\right ) e^{2 \alpha x}-2 a b \beta (3 \alpha -\gamma ) e^{\alpha x+\beta y}+2 b^2 \beta ^2 e^{2 \beta y}\right )\right )}{a^3 (\alpha -\gamma ) (2 \alpha -\gamma ) (3 \alpha -\gamma )}\right \}\right \}\]

Maple

restart; 
pde :=a*exp(alpha*x)*diff(w(x,y),x) +b*exp(beta*y)*diff(w(x,y),y) =c*exp(gamma*x-2*beta*y); 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
 

\[w \left ( x,y \right ) ={\frac {{b}^{2}{\beta }^{2}c}{{a}^{3}{\alpha }^{2}} \left ( {\frac { \left ( a\alpha \,{{\rm e}^{\alpha \,x}}-{{\rm e}^{\beta \,y}}b\beta \right ) ^{2}{{\rm e}^{1/8\, \left ( -24\,\alpha +1 \right ) x-2\,\beta \,y}}}{{b}^{2}{\beta }^{2} \left ( 1/8-\alpha \right ) }}+2\,{\frac { \left ( a\alpha \,{{\rm e}^{\alpha \,x}}-{{\rm e}^{\beta \,y}}b\beta \right ) {{\rm e}^{1/8\, \left ( -24\,\alpha +1 \right ) x-\beta \,y}}}{\beta \, \left ( 1/8-2\,\alpha \right ) b}}+{\frac {{{\rm e}^{x/8-3\,\alpha \,x}}}{1/8-3\,\alpha }} \right ) }+{\it \_F1} \left ( -{\frac { \left ( a\alpha \,{{\rm e}^{\alpha \,x}}-{{\rm e}^{\beta \,y}}b\beta \right ) {{\rm e}^{-\alpha \,x-\beta \,y}}}{\alpha \,b\beta }} \right ) \]

____________________________________________________________________________________

6.3.6.7 [884] Problem 7

problem number 884

Added Feb. 9, 2019.

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

Solve for \(w(x,y)\)

\[ a e^{\alpha x} w_x + b e^{\beta y} w_y = c e^{\gamma x} + s e^{\mu y} \]

Mathematica

ClearAll["Global`*"]; 
pde =  a*Exp[alpha*x]*D[w[x, y], x] + b*Exp[beta*y]*D[w[x, y], y] == c*Exp[gamma*x] + s*Exp[mu*y]; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
 

\[\left \{\left \{w(x,y)\to c_1\left (\frac {b e^{-\alpha x}}{a \alpha }-\frac {e^{-\beta y}}{\beta }\right )-\frac {c e^{\gamma x-\alpha x}}{a \alpha -a \gamma }-\frac {s \left (e^{-\beta y}\right )^{1-\frac {\mu }{\beta }}}{b \beta -b \mu }\right \}\right \}\]

Maple

restart; 
pde :=a*exp(alpha*x)*diff(w(x,y),x) +b*exp(beta*y)*diff(w(x,y),y) =c*exp(gamma*x) + s*exp(mu*y); 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
 

\[w \left ( x,y \right ) ={\frac {1}{b\alpha \,a \left ( -1+8\,\alpha \right ) \left ( -\beta +\mu \right ) } \left ( 8\,s \left ( \left ( a\alpha \,{{\rm e}^{\alpha \,x}}-{{\rm e}^{\beta \,y}}b\beta \right ) {{\rm e}^{-\alpha \,x-\beta \,y}}+{{\rm e}^{-\alpha \,x}}\beta \,b \right ) \left ( -1/8+\alpha \right ) \left ( {\frac {a\alpha }{ \left ( a\alpha \,{{\rm e}^{\alpha \,x}}-{{\rm e}^{\beta \,y}}b\beta \right ) {{\rm e}^{-\alpha \,x-\beta \,y}}+{{\rm e}^{-\alpha \,x}}\beta \,b}} \right ) ^{{\frac {\mu }{\beta }}}+8\, \left ( -\beta +\mu \right ) b\alpha \, \left ( a \left ( -1/8+\alpha \right ) {\it \_F1} \left ( -{\frac { \left ( a\alpha \,{{\rm e}^{\alpha \,x}}-{{\rm e}^{\beta \,y}}b\beta \right ) {{\rm e}^{-\alpha \,x-\beta \,y}}}{\alpha \,b\beta }} \right ) -c{{\rm e}^{x/8-\alpha \,x}} \right ) \right ) }\]

____________________________________________________________________________________

6.3.6.8 [885] Problem 8

problem number 885

Added Feb. 9, 2019.

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

Solve for \(w(x,y)\)

\[ a e^{\beta x} w_x + (b e^{\gamma x}+ c e^{\lambda y} ) w_y = s e^{\mu x} + k e^{\delta y} + p \]

Mathematica

ClearAll["Global`*"]; 
pde =  a*Exp[beta*x]*D[w[x, y], x] + (b*Exp[gamma*x] + c*Exp[lambda*y])*D[w[x, y], y] == s*Exp[mu*x] + k*Exp[delta*y] + p; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
 

Failed

Maple

restart; 
pde :=a*exp(beta*x)*diff(w(x,y),x) +(b*exp(gamma*x)+c*exp(lambda*y))*diff(w(x,y),y) =s*exp(mu*x) + k*exp(delta*y)+p; 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
 

\[w \left ( x,y \right ) =\int ^{x}\!{\frac {1}{a} \left ( \left ( {a \left ( \lambda \,a\int \!{\frac {c}{a}{{\rm e}^{{\frac {-8\,\lambda \,b{{\rm e}^{-\beta \,x+x/8}}-8\,\beta \, \left ( -1/8+\beta \right ) xa}{ \left ( 8\,\beta -1 \right ) a}}}}}\,{\rm d}x-\int \!{{\rm e}^{{\frac {-8\,b\lambda \,{{\rm e}^{-\beta \,{\it \_b}+{\it \_b}/8}}-8\,\beta \, \left ( -1/8+\beta \right ) {\it \_b}\,a}{ \left ( 8\,\beta -1 \right ) a}}}}\,{\rm d}{\it \_b}c\lambda +a{{\rm e}^{-8\,{\frac { \left ( b{{\rm e}^{-\beta \,x+x/8}}+ay \left ( -1/8+\beta \right ) \right ) \lambda }{ \left ( 8\,\beta -1 \right ) a}}}} \right ) ^{-1}} \right ) ^{{\frac {\delta }{\lambda }}}k{{\rm e}^{{\frac {-8\,\delta \,b{{\rm e}^{-\beta \,{\it \_b}+{\it \_b}/8}}-8\,\beta \, \left ( -1/8+\beta \right ) {\it \_b}\,a}{ \left ( 8\,\beta -1 \right ) a}}}}+s{{\rm e}^{{\it \_b}\, \left ( -\beta +\mu \right ) }}+{{\rm e}^{-\beta \,{\it \_b}}}p \right ) }{d{\it \_b}}+{\it \_F1} \left ( {\frac {1}{\lambda } \left ( -\lambda \,\int \!{\frac {c}{a}{{\rm e}^{{\frac {-8\,\lambda \,b{{\rm e}^{-\beta \,x+x/8}}-8\,\beta \, \left ( -1/8+\beta \right ) xa}{ \left ( 8\,\beta -1 \right ) a}}}}}\,{\rm d}x-{{\rm e}^{-8\,{\frac { \left ( b{{\rm e}^{-\beta \,x+x/8}}+ay \left ( -1/8+\beta \right ) \right ) \lambda }{ \left ( 8\,\beta -1 \right ) a}}}} \right ) } \right ) \]

____________________________________________________________________________________

6.3.6.9 [886] Problem 9

problem number 886

Added Feb. 9, 2019.

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

Solve for \(w(x,y)\)

\[ a e^{\beta x} w_x + (b e^{\gamma x}+ c e^{\lambda y} ) w_y = s e^{\mu x+\delta y} + k \]

Mathematica

ClearAll["Global`*"]; 
pde =  a*Exp[beta*x]*D[w[x, y], x] + (b*Exp[gamma*x] + c*Exp[lambda*y])*D[w[x, y], y] == s*Exp[mu*x + delta*y] + k; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
 

Failed

Maple

restart; 
pde :=a*exp(beta*x)*diff(w(x,y),x) +(b*exp(gamma*x)+c*exp(lambda*y))*diff(w(x,y),y) =s*exp(mu*x+delta*y)+k; 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
 

\[w \left ( x,y \right ) =\int ^{x}\!{\frac {1}{a} \left ( \left ( {a \left ( \lambda \,a\int \!{\frac {c}{a}{{\rm e}^{{\frac {-8\,\lambda \,b{{\rm e}^{-\beta \,x+x/8}}-8\,\beta \, \left ( -1/8+\beta \right ) xa}{ \left ( 8\,\beta -1 \right ) a}}}}}\,{\rm d}x-\int \!{{\rm e}^{{\frac {-8\,b\lambda \,{{\rm e}^{-\beta \,{\it \_b}+{\it \_b}/8}}-8\,\beta \, \left ( -1/8+\beta \right ) {\it \_b}\,a}{ \left ( 8\,\beta -1 \right ) a}}}}\,{\rm d}{\it \_b}c\lambda +a{{\rm e}^{-8\,{\frac { \left ( b{{\rm e}^{-\beta \,x+x/8}}+ay \left ( -1/8+\beta \right ) \right ) \lambda }{ \left ( 8\,\beta -1 \right ) a}}}} \right ) ^{-1}} \right ) ^{{\frac {\delta }{\lambda }}}s{{\rm e}^{{\frac {-8\,\delta \,b{{\rm e}^{-\beta \,{\it \_b}+{\it \_b}/8}}+8\, \left ( -1/8+\beta \right ) \left ( -\beta +\mu \right ) {\it \_b}\,a}{ \left ( 8\,\beta -1 \right ) a}}}}+{{\rm e}^{-\beta \,{\it \_b}}}k \right ) }{d{\it \_b}}+{\it \_F1} \left ( {\frac {1}{\lambda } \left ( -\lambda \,\int \!{\frac {c}{a}{{\rm e}^{{\frac {-8\,\lambda \,b{{\rm e}^{-\beta \,x+x/8}}-8\,\beta \, \left ( -1/8+\beta \right ) xa}{ \left ( 8\,\beta -1 \right ) a}}}}}\,{\rm d}x-{{\rm e}^{-8\,{\frac { \left ( b{{\rm e}^{-\beta \,x+x/8}}+ay \left ( -1/8+\beta \right ) \right ) \lambda }{ \left ( 8\,\beta -1 \right ) a}}}} \right ) } \right ) \]

____________________________________________________________________________________

6.3.6.10 [887] Problem 10

problem number 887

Added Feb. 9, 2019.

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

Solve for \(w(x,y)\)

\[ a e^{\beta x} w_x + b e^{\gamma x+\lambda y} w_y = c e^{\mu x+\delta y} + k \]

Mathematica

ClearAll["Global`*"]; 
pde =  a*Exp[beta*x]*D[w[x, y], x] + b*Exp[gamma*x + lambda*y]*D[w[x, y], y] == c*Exp[mu*x + delta*y] + k; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
 

\[\left \{\left \{w(x,y)\to c_1\left (\frac {b e^{\gamma x-\beta x}}{a \beta -a \gamma }-\frac {e^{-\lambda y}}{\lambda }\right )+\frac {c (\gamma -\beta ) \left (e^{\lambda y}\right )^{\delta /\lambda } e^{-\gamma x-\lambda y+\mu x} \, _2F_1\left (1,\frac {\mu -\gamma }{\beta -\gamma };\frac {\beta \delta -\gamma \delta -\gamma \lambda +\lambda \mu }{\beta \lambda -\gamma \lambda };1-\frac {a e^{\beta x-\gamma x-\lambda y} (\beta -\gamma )}{b \lambda }\right )}{b (\beta (\lambda -\delta )+\delta \gamma -\lambda \mu )}-\frac {k e^{-\beta x}}{a \beta }\right \}\right \}\]

Maple

restart; 
pde :=a*exp(beta*x)*diff(w(x,y),x) +b*exp(gamma*x+lambda*y)*diff(w(x,y),y) =c*exp(mu*x+delta*y)+k; 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
 

\[w \left ( x,y \right ) =\int ^{x}\!{\frac {1}{a} \left ( c \left ( {\frac {a \left ( 8\,\beta -1 \right ) ^{2}}{ \left ( 64\,\beta -8 \right ) \left ( \left ( \left ( -1/8+\beta \right ) a{{\rm e}^{\beta \,x-x/8}}-{{\rm e}^{y\lambda }}\lambda \,b \right ) {{\rm e}^{-\beta \,x-y\lambda +x/8}}+b\lambda \,{{\rm e}^{-\beta \,{\it \_a}+{\it \_a}/8}} \right ) }} \right ) ^{{\frac {\delta }{\lambda }}}{{\rm e}^{{\it \_a}\, \left ( -\beta +\mu \right ) }}+{{\rm e}^{-\beta \,{\it \_a}}}k \right ) }{d{\it \_a}}+{\it \_F1} \left ( -8\,{\frac { \left ( \left ( -1/8+\beta \right ) a{{\rm e}^{\beta \,x-x/8}}-{{\rm e}^{y\lambda }}\lambda \,b \right ) {{\rm e}^{-\beta \,x-y\lambda +x/8}}}{\lambda \,b \left ( 8\,\beta -1 \right ) }} \right ) \]

____________________________________________________________________________________

6.3.6.11 [888] Problem 11

problem number 888

Added Feb. 9, 2019.

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

Solve for \(w(x,y)\)

\[ a e^{\lambda x} w_x + b e^{\beta x} w_y = c e^{\gamma x} + d \]

Mathematica

ClearAll["Global`*"]; 
pde =  a*Exp[lambda*y]*D[w[x, y], x] + b*Exp[beta*x]*D[w[x, y], y] == c*Exp[gamma*y] + d; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
 

\[\left \{\left \{w(x,y)\to c_1\left (\frac {e^{\lambda y}}{\lambda }-\frac {b e^{\beta x}}{a \beta }\right )-\frac {d \gamma \log \left (\frac {a \beta e^{\lambda y}}{\lambda }\right )-\beta d \gamma x+c \lambda \left (e^{\lambda y}\right )^{\gamma /\lambda }}{a \beta \gamma e^{\lambda y}-b \gamma \lambda e^{\beta x}}-\frac {a \beta c \lambda \left (e^{\lambda y}\right )^{\frac {\gamma +\lambda }{\lambda }} \, _2F_1\left (1,\frac {\gamma +\lambda }{\lambda };\frac {\gamma }{\lambda }+2;\frac {a \beta e^{\lambda y}}{a \beta e^{\lambda y}-b e^{\beta x} \lambda }\right )}{(\gamma +\lambda ) \left (a \beta e^{\lambda y}-b \lambda e^{\beta x}\right )^2}\right \}\right \}\]

Maple

restart; 
pde :=a*exp(lambda*y)*diff(w(x,y),x) +b*exp(beta*x)*diff(w(x,y),y) =c*exp(gamma*y)+d; 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
 

\[w \left ( x,y \right ) =\int ^{x}\!{\frac {\beta }{{{\rm e}^{y\lambda }}a\beta -b\lambda \, \left ( {{\rm e}^{\beta \,x}}-{{\rm e}^{\beta \,{\it \_a}}} \right ) } \left ( c \left ( {\frac {{{\rm e}^{y\lambda }}a\beta -b\lambda \, \left ( {{\rm e}^{\beta \,x}}-{{\rm e}^{\beta \,{\it \_a}}} \right ) }{a\beta }} \right ) ^{1/8\,{\lambda }^{-1}}+d \right ) }{d{\it \_a}}+{\it \_F1} \left ( {\frac {-{{\rm e}^{\beta \,x}}b\lambda +{{\rm e}^{y\lambda }}a\beta }{b\beta \,\lambda }} \right ) \]

____________________________________________________________________________________