Added April 1, 2019.
Problem Chapter 5.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}+s e^{\mu y}) w + k e^{\nu x} \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*D[w[x, y], x] + b*D[w[x, y], y] == (c*Exp[lambda*x]+s*Exp[mu*y])*w[x,y] + k*Exp[nu*x]; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to e^{\frac {c e^{\lambda x}}{a \lambda }+\frac {s e^{\mu y}}{b \mu }} \left (\int _1^x\frac {\exp \left (-\frac {e^{\lambda K[1]} c}{a \lambda }-\frac {e^{\mu \left (y+\frac {b (K[1]-x)}{a}\right )} s}{b \mu }+\nu K[1]\right ) k}{a}dK[1]+c_1\left (y-\frac {b x}{a}\right )\right )\right \}\right \}\]
Maple ✓
restart; pde := a*diff(w(x,y),x)+ b*diff(w(x,y),y) = (c*exp(lambda*x)+s*exp(mu*y))*w(x,y)+ k*exp(nu*x); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime'));
\[w \left ( x,y \right ) = \left ( \int ^{x}\!{\frac {k}{a}{{\rm e}^{{\frac {1}{a\lambda \,b\mu } \left ( -as\lambda \,{{\rm e}^{{\frac { \left ( ya-b \left ( x-{\it \_a} \right ) \right ) \mu }{a}}}}+b\mu \, \left ( {\it \_a}\,\nu \,\lambda \,a-{{\rm e}^{{\it \_a}\,\lambda }}c \right ) \right ) }}}}{d{\it \_a}}+{\it \_F1} \left ( {\frac {ya-bx}{a}} \right ) \right ) {{\rm e}^{{\frac {{{\rm e}^{\lambda \,x}}cb\mu +as\lambda \,{{\rm e}^{\mu \,y}}}{a\lambda \,b\mu }}}}\]
____________________________________________________________________________________
Added April 1, 2019.
Problem Chapter 5.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^{\alpha x+\beta y} w+ k e^{\gamma x} \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*D[w[x, y], x] + b*D[w[x, y], y] == c*Exp[alpha*x+beta*y]*w[x,y] + k*Exp[gamma*x]; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to e^{\frac {c e^{\alpha x+\beta y}}{a \alpha +b \beta }} \left (\int _1^x\frac {\exp \left (\gamma K[1]-\frac {c e^{\beta y+\alpha K[1]+\frac {b \beta (K[1]-x)}{a}}}{a \alpha +b \beta }\right ) k}{a}dK[1]+c_1\left (y-\frac {b x}{a}\right )\right )\right \}\right \}\]
Maple ✓
restart; pde := a*diff(w(x,y),x)+ b*diff(w(x,y),y) = c*exp(alpha*x+beta*y)*w(x,y)+ k*exp(gamma*x); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime'));
\[w \left ( x,y \right ) = \left ( \int ^{x}\!{\frac {k}{a}{{\rm e}^{{\frac {1}{8\,a\alpha +8\,b\beta } \left ( -8\,c{{\rm e}^{{\frac { \left ( {\it \_a}\,\alpha +\beta \,y \right ) a-b\beta \, \left ( x-{\it \_a} \right ) }{a}}}}+{\it \_a}\, \left ( a\alpha +b\beta \right ) \right ) }}}}{d{\it \_a}}+{\it \_F1} \left ( {\frac {ya-bx}{a}} \right ) \right ) {{\rm e}^{{\frac {c{{\rm e}^{\alpha \,x+\beta \,y}}}{a\alpha +b\beta }}}}\]
____________________________________________________________________________________
Added April 1, 2019.
Problem Chapter 5.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 x} w_y = c e^{\gamma y} w+ s e^{\mu x+\delta y} \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*Exp[lambda*x]*D[w[x, y], x] + b*Exp[beta*x]*D[w[x, y], y] == c*Exp[gamma*y]*w[x,y] + s*Exp[mu*x+delta*y]; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to \exp \left (\int _1^x\frac {c \exp \left (y \gamma -\frac {b \left (e^{(\beta -\lambda ) x}-e^{(\beta -\lambda ) K[1]}\right ) \gamma }{a (\beta -\lambda )}-\lambda K[1]\right )}{a}dK[1]\right ) \left (\int _1^x\frac {\exp \left (-\frac {b \delta \left (e^{(\beta -\lambda ) x}-e^{(\beta -\lambda ) K[2]}\right )}{a (\beta -\lambda )}+\delta y+(\mu -\lambda ) K[2]-\int _1^{K[2]}\frac {c \exp \left (y \gamma -\frac {b \left (e^{(\beta -\lambda ) x}-e^{(\beta -\lambda ) K[1]}\right ) \gamma }{a (\beta -\lambda )}-\lambda K[1]\right )}{a}dK[1]\right ) s}{a}dK[2]+c_1\left (\frac {b e^{x (\beta -\lambda )}}{a (\lambda -\beta )}+y\right )\right )\right \}\right \}\]
Maple ✓
restart; pde := a*exp(lambda*x)*diff(w(x,y),x)+ b*exp(beta*x)*diff(w(x,y),y) = c*exp(gamma*y)*w(x,y)+ s*exp(mu*x+delta*y); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime'));
\[w \left ( x,y \right ) = \left ( \int ^{x}\!{\frac {s}{a}{{\rm e}^{{\frac {1}{a \left ( -\beta +\lambda \right ) } \left ( -c \left ( -\beta +\lambda \right ) \int \!{{\rm e}^{1/8\,{\frac {-b{{\rm e}^{x \left ( \beta -\lambda \right ) }}+{{\rm e}^{{\it \_b}\, \left ( \beta -\lambda \right ) }}b+a \left ( \beta -\lambda \right ) \left ( -8\,{\it \_b}\,\lambda +y \right ) }{ \left ( \beta -\lambda \right ) a}}}}\,{\rm d}{\it \_b}-{{\rm e}^{{\it \_b}\, \left ( \beta -\lambda \right ) }}b\delta +{{\rm e}^{x \left ( \beta -\lambda \right ) }}b\delta -a \left ( -\beta +\lambda \right ) \left ( {\it \_b}\,\lambda -\mu \,{\it \_b}-\delta \,y \right ) \right ) }}}}{d{\it \_b}}+{\it \_F1} \left ( {\frac {-b{{\rm e}^{x \left ( \beta -\lambda \right ) }}+ay \left ( \beta -\lambda \right ) }{ \left ( \beta -\lambda \right ) a}} \right ) \right ) {{\rm e}^{\int ^{x}\!{\frac {c}{a}{{\rm e}^{1/8\,{\frac {-b{{\rm e}^{x \left ( \beta -\lambda \right ) }}+{{\rm e}^{{\it \_a}\, \left ( \beta -\lambda \right ) }}b+a \left ( \beta -\lambda \right ) \left ( -8\,{\it \_a}\,\lambda +y \right ) }{ \left ( \beta -\lambda \right ) a}}}}}{d{\it \_a}}}}\]
____________________________________________________________________________________
Added April 1, 2019.
Problem Chapter 5.3.1.4, 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 w+k e^{\mu x+\delta y} \]
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*w[x,y] + k*Exp[mu*x+delta*y]; 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*w(x,y)+ k*exp(mu*x+delta*y); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime'));
\[w \left ( x,y \right ) = \left ( \int ^{x}\!{\frac {k}{a} \left ( {a \left ( 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\lambda -c\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}\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 }}}{{\rm e}^{{\frac {-8\,{{\rm e}^{-\beta \,{\it \_b}+{\it \_b}/8}}b\beta \,\delta +8\, \left ( s{{\rm e}^{-\beta \,{\it \_b}}}+a{\it \_b}\,\beta \, \left ( \mu -\beta \right ) \right ) \left ( -1/8+\beta \right ) }{\beta \, \left ( 8\,\beta -1 \right ) a}}}}}{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 ) \right ) {{\rm e}^{-{\frac {s{{\rm e}^{-\beta \,x}}}{a\beta }}}}\]
____________________________________________________________________________________
Added April 1, 2019.
Problem Chapter 5.3.1.5, 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} w + 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]*w[x,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)*w(x,y)+k; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime'));
\[w \left ( x,y \right ) = \left ( \int ^{x}\!{\frac {k}{a}{{\rm e}^{{\frac {1}{a} \left ( -\beta \,{\it \_f}\,a-s\int \! \left ( {a \left ( 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\lambda -c\int \!{{\rm e}^{{\frac {-8\,b\lambda \,{{\rm e}^{-{\it \_f}\,\beta +{\it \_f}/8}}-8\,\beta \, \left ( -1/8+\beta \right ) {\it \_f}\,a}{ \left ( 8\,\beta -1 \right ) a}}}}\,{\rm d}{\it \_f}\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 }}}{{\rm e}^{{\frac {-8\,\delta \,b{{\rm e}^{-{\it \_f}\,\beta +{\it \_f}/8}}+8\, \left ( -1/8+\beta \right ) {\it \_f}\, \left ( \mu -\beta \right ) a}{ \left ( 8\,\beta -1 \right ) a}}}}\,{\rm d}{\it \_f} \right ) }}}}{d{\it \_f}}+{\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 ) \right ) {{\rm e}^{\int ^{x}\!{\frac {s}{a} \left ( {a \left ( 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\lambda -c\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}\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 }}}{{\rm e}^{{\frac {-8\,\delta \,b{{\rm e}^{-\beta \,{\it \_b}+{\it \_b}/8}}+8\, \left ( -1/8+\beta \right ) {\it \_b}\, \left ( \mu -\beta \right ) a}{ \left ( 8\,\beta -1 \right ) a}}}}}{d{\it \_b}}}}\]
____________________________________________________________________________________
Added April 1, 2019.
Problem Chapter 5.3.1.6, 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^{\sigma y} w + k e^{\mu x+delta y} + d \]
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[sigma*y]*w[x,y]+k*Exp[mu*x+delta*y]+d; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to \exp \left (\frac {c (\gamma -\beta ) \left (e^{\lambda y}\right )^{\frac {\sigma }{\lambda }} e^{-\gamma x-\lambda y} \, _2F_1\left (1,-\frac {\gamma }{\beta -\gamma };\frac {\beta \sigma -\gamma (\lambda +\sigma )}{(\beta -\gamma ) \lambda };1-\frac {a e^{\beta x-\gamma x-\lambda y} (\beta -\gamma )}{b \lambda }\right )}{b (\beta (\lambda -\sigma )+\gamma \sigma )}\right ) \left (\int _1^x\frac {\exp \left (\frac {c e^{-\gamma K[1]} (\beta -\gamma ) \left (\frac {a e^{\lambda y+\beta (x+K[1])} (\beta -\gamma )}{a e^{\beta (x+K[1])} (\beta -\gamma )-b e^{\lambda y} \left (e^{\gamma x+\beta K[1]}-e^{\beta x+\gamma K[1]}\right ) \lambda }\right )^{\frac {\sigma }{\lambda }-1} \, _2F_1\left (1,-\frac {\gamma }{\beta -\gamma };\frac {\beta \sigma -\gamma (\lambda +\sigma )}{(\beta -\gamma ) \lambda };e^{(\beta -\gamma ) K[1]} \left (\frac {a e^{-\lambda y} (\gamma -\beta )}{b \lambda }+e^{(\gamma -\beta ) x}\right )\right )}{b (\beta (\lambda -\sigma )+\gamma \sigma )}-\beta K[1]\right ) \left (e^{\mu K[1]} k \left (\frac {a e^{\lambda y+\beta (x+K[1])} (\gamma -\beta )}{b e^{\lambda y} \left (e^{\gamma x+\beta K[1]}-e^{\beta x+\gamma K[1]}\right ) \lambda -a e^{\beta (x+K[1])} (\beta -\gamma )}\right )^{\delta /\lambda }+d\right )}{a}dK[1]+c_1\left (\frac {b e^{\gamma x-\beta x}}{a \beta -a \gamma }-\frac {e^{-\lambda y}}{\lambda }\right )\right )\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(sigma*y)*w(x,y)+k*exp(mu*x+delta*y)+d; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime'));
\[w \left ( x,y \right ) = \left ( \int ^{x}\!{\frac {1}{a} \left ( k \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+x/8-y\lambda }}+b\lambda \,{{\rm e}^{-\beta \,{\it \_b}+{\it \_b}/8}} \right ) }} \right ) ^{{\frac {\delta }{\lambda }}}{{\rm e}^{{\frac {1}{a} \left ( -c\int \! \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+x/8-y\lambda }}+b\lambda \,{{\rm e}^{-\beta \,{\it \_b}+{\it \_b}/8}} \right ) }} \right ) ^{{\frac {\sigma }{\lambda }}}{{\rm e}^{-\beta \,{\it \_b}}}\,{\rm d}{\it \_b}+a{\it \_b}\, \left ( \mu -\beta \right ) \right ) }}}+d{{\rm e}^{{\frac {1}{a} \left ( -\beta \,{\it \_b}\,a-c\int \! \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+x/8-y\lambda }}+b\lambda \,{{\rm e}^{-\beta \,{\it \_b}+{\it \_b}/8}} \right ) }} \right ) ^{{\frac {\sigma }{\lambda }}}{{\rm e}^{-\beta \,{\it \_b}}}\,{\rm d}{\it \_b} \right ) }}} \right ) }{d{\it \_b}}+{\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+x/8-y\lambda }}}{b\lambda \, \left ( 8\,\beta -1 \right ) }} \right ) \right ) {{\rm e}^{\int ^{x}\!{\frac {c{{\rm e}^{-\beta \,{\it \_a}}}}{a} \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+x/8-y\lambda }}+b\lambda \,{{\rm e}^{-\beta \,{\it \_a}+{\it \_a}/8}} \right ) }} \right ) ^{{\frac {\sigma }{\lambda }}}}{d{\it \_a}}}}\]
____________________________________________________________________________________
Added April 1, 2019.
Problem Chapter 5.3.1.7, 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 w + s e^{\gamma x} \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*Exp[lambda*x]*D[w[x, y], x] + b*Exp[beta*x]*D[w[x, y], y] == c*w[x,y]+s*Exp[gamma*x]; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to e^{-\frac {c e^{-\lambda x}}{a \lambda }} \left (\int _1^x\frac {e^{\frac {e^{-\lambda K[1]} c}{a \lambda }+(\gamma -\lambda ) K[1]} s}{a}dK[1]+c_1\left (\frac {b e^{x (\beta -\lambda )}}{a (\lambda -\beta )}+y\right )\right )\right \}\right \}\]
Maple ✓
restart; pde := a*exp(lambda*x)*diff(w(x,y),x)+ b*exp(beta*x)*diff(w(x,y),y) = c*w(x,y)+s*exp(gamma*x); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime'));
\[w \left ( x,y \right ) = \left ( \int \!{\frac {s}{a}{{\rm e}^{1/8\,{\frac {8\,c{{\rm e}^{-\lambda \,x}}-8\,\lambda \, \left ( \lambda -1/8 \right ) xa}{a\lambda }}}}}\,{\rm d}x+{\it \_F1} \left ( {\frac {-b{{\rm e}^{x \left ( \beta -\lambda \right ) }}+ay \left ( \beta -\lambda \right ) }{ \left ( \beta -\lambda \right ) a}} \right ) \right ) {{\rm e}^{-{\frac {c{{\rm e}^{-\lambda \,x}}}{a\lambda }}}}\]
____________________________________________________________________________________
Added April 1, 2019.
Problem Chapter 5.3.1.8, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ a e^{\lambda y} w_x + b x^{\beta x} w_y = c e^{\gamma x} w + s \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*Exp[lambda*x]*D[w[x, y], x] + b*x^(beta*x)*D[w[x, y], y] == c*Exp[gamma*x]*w[x,y]+s; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to e^{\frac {c e^{x (\gamma -\lambda )}}{a (\gamma -\lambda )}} \left (c_1\left (y-\int _1^x\frac {b e^{-\lambda K[1]} K[1]^{\beta K[1]}}{a}dK[1]\right )+\int _1^x\frac {\exp \left (-\frac {e^{(\gamma -\lambda ) K[2]} c}{a (\gamma -\lambda )}-\lambda K[2]\right ) s}{a}dK[2]\right )\right \}\right \}\]
Maple ✓
restart; pde := a*exp(lambda*x)*diff(w(x,y),x)+ b*x^(beta*x)*diff(w(x,y),y) = c*exp(gamma*x)*w(x,y)+s; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime'));
\[w \left ( x,y \right ) = \left ( \int \!{\frac {s}{a}{{\rm e}^{{\frac {8\,c{{\rm e}^{x/8-\lambda \,x}}-8\,\lambda \, \left ( \lambda -1/8 \right ) xa}{8\,a\lambda -a}}}}}\,{\rm d}x+{\it \_F1} \left ( {\frac {ya-b\int \!{x}^{\beta \,x}{{\rm e}^{-\lambda \,x}}\,{\rm d}x}{a}} \right ) \right ) {{\rm e}^{-8\,{\frac {c{{\rm e}^{x/8-\lambda \,x}}}{ \left ( 8\,\lambda -1 \right ) a}}}}\]
____________________________________________________________________________________