Added Feb. 25, 2019.
Problem Chapter 4.5.2.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 x^n+s \ln ^k(\lambda y)) w \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*D[w[x, y], x] + b*D[w[x, y], y] == (c*x^n + s*Log[gamma*y]^k)*w[x, 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 ) \exp \left (\frac {s \log ^k(\gamma y) (-\log (\gamma y))^{-k} \text {Gamma}(k+1,-\log (\gamma y))}{b \gamma }+\frac {c x^{n+1}}{a n+a}\right )\right \}\right \}\]
Maple ✓
restart; pde := a*diff(w(x,y),x)+b*diff(w(x,y),y) = (c*x^n+s*ln(gamma*y)^k)*w(x,y); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime'));
\[w \left ( x,y \right ) ={\it \_F1} \left ( {\frac {ay-xb}{a}} \right ) {{\rm e}^{\int ^{x}\!{\frac {1}{a} \left ( c{{\it \_a}}^{n}+s \left ( \ln \left ( \gamma \right ) +\ln \left ( {\frac {ay-b \left ( x-{\it \_a} \right ) }{a}} \right ) \right ) ^{k} \right ) }{d{\it \_a}}}}\]
____________________________________________________________________________________
Added Feb. 25, 2019.
Problem Chapter 4.5.2.2, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ w_x + a w_y = (b y^2+c x^n y+ s \ln ^k(\lambda x)) w \]
Mathematica ✓
ClearAll["Global`*"]; pde = D[w[x, y], x] + a*D[w[x, y], y] == (b*y^2 + c*x^n*y + s*Log[lambda*x]^k)*w[x, y]; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to c_1(y-a x) \exp \left (\frac {s \log ^k(\lambda x) (-\log (\lambda x))^{-k} \text {Gamma}(k+1,-\log (\lambda x))}{\lambda }+\frac {1}{3} a^2 b x^3+a b x^2 (y-a x)+b x (y-a x)^2+\frac {c x^{n+1} ((n+2) y-a x)}{(n+1) (n+2)}\right )\right \}\right \}\]
Maple ✓
restart; pde := diff(w(x,y),x)+a*diff(w(x,y),y) = (b*y^2+c*x^n*y+ s*ln(lambda*x)^k)*w(x,y); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime'));
\[w \left ( x,y \right ) ={\it \_F1} \left ( -ax+y \right ) {{\rm e}^{\int ^{x}\!s \left ( \ln \left ( {\it \_a}\,\lambda \right ) \right ) ^{k}+c{{\it \_a}}^{1+n}a-c \left ( ax-y \right ) {{\it \_a}}^{n}+b \left ( \left ( x-{\it \_a} \right ) a-y \right ) ^{2}{d{\it \_a}}}}\]
____________________________________________________________________________________
Added March 9, 2019.
Problem Chapter 4.5.2.3, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ w_x + a w_y = b \ln ^k(\lambda x) \ln ^n(\beta y) w \]
Mathematica ✓
ClearAll["Global`*"]; pde = D[w[x, y], x] + a*D[w[x, y], y] == b*Log[lambda*x]^k*Log[beta*y]^n*w[x, y]; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to c_1(y-a x) \exp \left (\int _1^xb \log ^k(\lambda K[1]) \log ^n(\beta (y+a (K[1]-x)))dK[1]\right )\right \}\right \}\]
Maple ✓
restart; pde := diff(w(x,y),x)+a*diff(w(x,y),y) = b*ln(lambda*x)^k*ln(beta*y)^n*w(x,y); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime'));
\[w \left ( x,y \right ) ={\it \_F1} \left ( -ax+y \right ) {{\rm e}^{\int ^{x}\!b \left ( \ln \left ( {\it \_a}\,\lambda \right ) \right ) ^{k} \left ( \ln \left ( - \left ( \left ( x-{\it \_a} \right ) a-y \right ) \beta \right ) \right ) ^{n}{d{\it \_a}}}}\]
____________________________________________________________________________________
Added March 9, 2019.
Problem Chapter 4.5.2.4, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ w_x + (a y+b x^n) w_y = c \ln ^k(\lambda x) w \]
Mathematica ✓
ClearAll["Global`*"]; pde = D[w[x, y], x] + (a*y + b*x^n)*D[w[x, y], y] == c*Log[lambda*x]^k*w[x, y]; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to c_1\left (b a^{-n-1} \text {Gamma}(n+1,a x)+y e^{-a x}\right ) \exp \left (\int _1^xc \log ^k(\lambda K[1])dK[1]\right )\right \}\right \}\]
Maple ✓
restart; pde := diff(w(x,y),x)+(a*y+b*x^n)*diff(w(x,y),y) = c*ln(lambda*x)^k*w(x,y); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime'));
\[w \left ( x,y \right ) ={\it \_F1} \left ( {\frac {{{\rm e}^{-ax}}}{ \left ( 1+n \right ) a} \left ( -{x}^{n} \left ( ax \right ) ^{-{\frac {n}{2}}} \WhittakerM \left ( {\frac {n}{2}},{\frac {n}{2}}+{\frac {1}{2}},ax \right ) {{\rm e}^{{\frac {ax}{2}}}}b+ay \left ( 1+n \right ) \right ) } \right ) {{\rm e}^{\int \!c \left ( \ln \left ( x\lambda \right ) \right ) ^{k}\,{\rm d}x}}\]
____________________________________________________________________________________
Added March 9, 2019.
Problem Chapter 4.5.2.5, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ a x w_x + b y w_y = x^k (n \ln x+ m \ln y) w \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*x*D[w[x, y], x] + b*y*D[w[x, y], y] == x^k*(n*Log[x] + m*Log[y])*w[x, y]; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to c_1\left (y x^{-\frac {b}{a}}\right ) \exp \left (\frac {x^k (a k m \log (y)+a k n \log (x)-a n-b m)}{a^2 k^2}\right )\right \}\right \}\]
Maple ✓
restart; pde := a*x*diff(w(x,y),x)+ b*y*diff(w(x,y),y) = x^k*(n*ln(x)+m*ln(y))*w(x,y); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime'));
\[w \left ( x,y \right ) ={\it \_F1} \left ( y{x}^{-{\frac {b}{a}}} \right ) \left ( {x}^{{\frac {b}{a}}} \right ) ^{{\frac {m{x}^{k}}{ak}}} \left ( y{x}^{-{\frac {b}{a}}} \right ) ^{{\frac {m{x}^{k}}{ak}}}{x}^{{\frac {{x}^{k}n}{ak}}}{{\rm e}^{-{\frac {{x}^{k}}{2\,{a}^{2}{k}^{2}} \left ( -im\pi \,k \left ( {\it csgn} \left ( iy \right ) -{\it csgn} \left ( i{x}^{{\frac {b}{a}}} \right ) \right ) a{\it csgn} \left ( iy \right ) {\it csgn} \left ( iy{x}^{-{\frac {b}{a}}} \right ) +i\pi \,m \left ( {\it csgn} \left ( iy \right ) \right ) ^{3}ak-i\pi \,m{\it csgn} \left ( i{x}^{{\frac {b}{a}}} \right ) \left ( {\it csgn} \left ( iy \right ) \right ) ^{2}ak+2\,na+2\,mb \right ) }}}\]
____________________________________________________________________________________
Added March 9, 2019.
Problem Chapter 4.5.2.6, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ a x^k w_x + b y^n w_y = (c \ln ^m(\lambda x)+s \ln ^t(\beta y)) w \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*x^k*D[w[x, y], x] + b*y^n*D[w[x, y], y] == (c*Log[lambda*x]^m + s*Log[beta*y]^t)*w[x, y]; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to c_1\left (\frac {b x^{1-k}}{a (k-1)}-\frac {y^{1-n}}{n-1}\right ) \exp \left (\int _1^x\frac {K[1]^{-k} \left (c \log ^m(\lambda K[1])+s \log ^t\left (\beta \left (\frac {a (k-1) x^k y^n K[1]^k}{a (k-1) x^k y K[1]^k-b (n-1) y^n \left (x K[1]^k-x^k K[1]\right )}\right )^{\frac {1}{n-1}}\right )\right )}{a}dK[1]\right )\right \}\right \}\]
Maple ✓
restart; pde := a*x^k*diff(w(x,y),x)+ b*y^n*diff(w(x,y),y) = (c*ln(lambda*x)^m+s*ln(beta*y)^t)*w(x,y); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime'));
\[w \left ( x,y \right ) ={\it \_F1} \left ( {\frac {-{x}^{-k+1}b \left ( n-1 \right ) +{y}^{1-n}a \left ( k-1 \right ) }{a \left ( k-1 \right ) }} \right ) {{\rm e}^{\int ^{x}\!{\frac {{{\it \_a}}^{-k}}{a} \left ( c \left ( \ln \left ( {\it \_a}\,\lambda \right ) \right ) ^{m}+s \left ( \ln \left ( \beta \, \left ( {\frac {-{x}^{-k+1}b \left ( n-1 \right ) +{y}^{1-n}a \left ( k-1 \right ) +{{\it \_a}}^{-k+1}b \left ( n-1 \right ) }{a \left ( k-1 \right ) }} \right ) ^{- \left ( n-1 \right ) ^{-1}} \right ) \right ) ^{t} \right ) }{d{\it \_a}}}}\]
____________________________________________________________________________________