Added Feb. 9, 2019.
Problem Chapter 3.3.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 y e^{\lambda x} + k x e^{\mu y} \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*D[w[x, y], x] + b*D[w[x, y], y] == c*y*Exp[lambda*x] + k*x*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 {b c e^{\lambda x}}{a^2 \lambda ^2}-\frac {a k e^{\mu y}}{b^2 \mu ^2}+\frac {c y e^{\lambda x}}{a \lambda }+\frac {k x e^{\mu y}}{b \mu }\right \}\right \}\]
Maple ✓
restart; pde :=a*diff(w(x,y),x) +b*diff(w(x,y),y) =c*y*exp(lambda*x)+k*x*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 {a^{2} b^{2} \lambda ^{2} \mu ^{2} \textit {\_F1} \left (\frac {a y -b x}{a}\right )-\left (-b \mu x +a \right ) a^{2} k \,\lambda ^{2} {\mathrm e}^{\mu y}+\left (a \lambda y -b \right ) b^{2} c \,\mu ^{2} {\mathrm e}^{\lambda x}}{a^{2} b^{2} \lambda ^{2} \mu ^{2}}\]
____________________________________________________________________________________
Added Feb. 9, 2019.
Problem Chapter 3.3.2.2 from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ w_x + a w_y = a x^k e^{\lambda y} \]
Mathematica ✓
ClearAll["Global`*"]; pde = D[w[x, y], x] + a*D[w[x, y], y] == a*x^k*Exp[lambda*y]; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to \frac {x^k (-a \lambda x)^{-k} e^{\lambda (y-a x)} \operatorname {Gamma}(k+1,-a \lambda x)}{\lambda }+c_1(y-a x)\right \}\right \}\]
Maple ✓
restart; pde := diff(w(x,y),x) +a*diff(w(x,y),y) =a*x^k*exp(lambda*y); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
\[w \left (x , y\right ) = \frac {\left (-\Gamma \left (k \right )+\Gamma \left (k , -a \lambda x \right )\right ) k \,x^{k} \left (-a \lambda x \right )^{-k} {\mathrm e}^{\left (-a x +y \right ) \lambda }+\lambda \textit {\_F1} \left (-a x +y \right )+x^{k} {\mathrm e}^{\lambda y}}{\lambda }\]
____________________________________________________________________________________
Added Feb. 9, 2019.
Problem Chapter 3.3.2.3 from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ w_x + (a y+b e^{\lambda x}) w_y = c e^{\beta x} \]
Mathematica ✓
ClearAll["Global`*"]; pde = D[w[x, y], x] + (a*y + b*Exp[lambda*x])*D[w[x, y], y] == c*Exp[beta*x]; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to \frac {c e^{\beta x}}{\beta }+c_1\left (e^{-a x} \left (\frac {b e^{\lambda x}}{a-\lambda }+y\right )\right )\right \}\right \}\]
Maple ✓
restart; pde := diff(w(x,y),x) +(a*y+b*exp(lambda*x))*diff(w(x,y),y) =c*exp(beta*x); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
\[w \left (x , y\right ) = \frac {\beta \textit {\_F1} \left (\frac {\left (b \,{\mathrm e}^{a x}+\left (a -\lambda \right ) y \,{\mathrm e}^{\left (a -\lambda \right ) x}\right ) {\mathrm e}^{-\left (2 a -\lambda \right ) x}}{a -\lambda }\right )+c \,{\mathrm e}^{\beta x}}{\beta }\]
____________________________________________________________________________________
Added Feb. 9, 2019.
Problem Chapter 3.3.2.4 from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ w_x + (a y e^{\lambda x}+b e^{\beta x} y^k) w_y = c e^{\mu x} \]
Mathematica ✓
ClearAll["Global`*"]; pde = D[w[x, y], x] + (a*y*Exp[lambda*x] + b*Exp[beta*x]*y^k)*D[w[x, y], y] == c*Exp[mu*x]; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to \frac {c e^{\mu x}}{\mu }+c_1\left ((k-1) \int _1^xb e^{\frac {a e^{\lambda K[1]} (k-1)}{\lambda }+\beta K[1]}dK[1]+y^{1-k} e^{\frac {a (k-1) e^{\lambda x}}{\lambda }}\right )\right \}\right \}\]
Maple ✓
restart; pde := diff(w(x,y),x) +(a*y*exp(lambda*x)+b*exp(beta*x)*y^k)*diff(w(x,y),y) =c*exp(mu*x); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
\[w \left (x , y\right ) = \frac {c \,{\mathrm e}^{\mu x}+\mu \textit {\_F1} \left (\left (k -1\right ) b \left (\int {\mathrm e}^{\frac {\beta \lambda x +\left (k -1\right ) a \,{\mathrm e}^{\lambda x}}{\lambda }}d x \right )+y^{-k +1} {\mathrm e}^{\frac {\left (k -1\right ) a \,{\mathrm e}^{\lambda x}}{\lambda }}\right )}{\mu }\]
____________________________________________________________________________________
Added Feb. 9, 2019.
Problem Chapter 3.3.2.5 from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ w_x + (a x^k+b x^n e^{\lambda y}) w_y = c e^{\beta x} \]
Mathematica ✓
ClearAll["Global`*"]; pde = D[w[x, y], x] + (a*x^k + b*x^n*Exp[lambda*y])*D[w[x, y], y] == c*Exp[beta*x]; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to \frac {c e^{\beta x}}{\beta }+c_1\left (\frac {b \lambda x^{n+1} \left (-\frac {a \lambda x^{k+1}}{k+1}\right )^{-\frac {n+1}{k+1}} \operatorname {Gamma}\left (\frac {n+1}{k+1},-\frac {a \lambda x^{k+1}}{k+1}\right )-(k+1) e^{-\frac {\lambda \left (-a x^{k+1}+k y+y\right )}{k+1}}}{a b (k+1) \lambda ^2 (k-n)}\right )\right \}\right \}\]
Maple ✓
restart; pde := diff(w(x,y),x) +(a*x^k+b*x^n*exp(lambda*y))*diff(w(x,y),y) =c*exp(beta*x); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
\[w \left (x , y\right ) = \frac {\beta \textit {\_F1} \left (\frac {-\left (a \lambda \,x^{n +1}+\left (-k -n -2\right ) x^{-k +n}\right ) \left (k +1\right )^{2} b \left (-\frac {a \lambda \,x^{k +1}}{k +1}\right )^{\frac {-k -n -2}{2 k +2}} \WhittakerM \left (\frac {-k +n}{2 k +2}, \frac {2 k +n +3}{2 k +2}, -\frac {a \lambda \,x^{k +1}}{k +1}\right ) {\mathrm e}^{\frac {a \lambda \,x^{k +1}}{2 k +2}}-2 \left (k +n +2\right ) \left (-\frac {\left (k +1\right ) \left (k +n +2\right ) b \,x^{-k +n} \left (-\frac {a \lambda \,x^{k +1}}{k +1}\right )^{\frac {-k -n -2}{2 k +2}} \WhittakerM \left (\frac {k +n +2}{2 k +2}, \frac {2 k +n +3}{2 k +2}, -\frac {a \lambda \,x^{k +1}}{k +1}\right ) {\mathrm e}^{\frac {a \lambda \,x^{k +1}}{2 k +2}}}{2}+\left (n +1\right ) \left (k +\frac {n}{2}+\frac {3}{2}\right ) a \,{\mathrm e}^{\frac {\left (a \,x^{k +1}-\left (k +1\right ) y \right ) \lambda }{k +1}}\right )}{\left (n +1\right ) \left (k +n +2\right ) \left (2 k +n +3\right ) a \lambda }\right )+c \,{\mathrm e}^{\beta x}}{\beta }\]
____________________________________________________________________________________
Added Feb. 9, 2019.
Problem Chapter 3.3.2.6 from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ x w_x + y w_y = a x e^{\lambda x+ \mu y} \]
Mathematica ✓
ClearAll["Global`*"]; pde = x*D[w[x, y], x] + y*D[w[x, y], y] == a*x*Exp[lambda*x + mu*y]; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to \frac {a x e^{\lambda x+\mu y}}{\lambda x+\mu y}+c_1\left (\frac {y}{x}\right )\right \}\right \}\]
Maple ✓
restart; pde := x*diff(w(x,y),x) +y*diff(w(x,y),y) =a*x*exp(lambda *x+ mu* y); cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
\[w \left (x , y\right ) = \frac {a \,{\mathrm e}^{\lambda x +\mu y}}{\lambda +\frac {\mu y}{x}}+\textit {\_F1} \left (\frac {y}{x}\right )\]
____________________________________________________________________________________
Added Feb. 9, 2019.
Problem Chapter 3.3.2.7 from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ x w_x + y w_y = a y e^{\lambda x} + b x e^{\mu y} \]
Mathematica ✓
ClearAll["Global`*"]; pde = x*D[w[x, y], x] + y*D[w[x, y], y] == a*y*Exp[lambda*x] + b*x*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 {y}{x}\right )+\frac {a y e^{\lambda x}}{\lambda x}+\frac {b x e^{\mu y}}{\mu y}\right \}\right \}\]
Maple ✓
restart; pde := x*diff(w(x,y),x) +y*diff(w(x,y),y) =a*y*exp(lambda*x) + b*x*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 {a y \,{\mathrm e}^{\lambda x}}{\lambda x}+\frac {b x \,{\mathrm e}^{\mu y}}{\mu y}+\textit {\_F1} \left (\frac {y}{x}\right )\]
____________________________________________________________________________________
Added Feb. 9, 2019.
Problem Chapter 3.3.2.8 from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ a x^k w_x + b e^{\lambda y} w_y = c x^n+s \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*x^k*D[w[x, y], x] + b*Exp[lambda*y]*D[w[x, y], y] == c*x^n + s; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to x^{1-k} \left (\frac {c x^n}{a (-k)+a n+a}+\frac {s}{a-a k}\right )+c_1\left (\frac {b x^{1-k}}{a (k-1)}-\frac {e^{-\lambda y}}{\lambda }\right )\right \}\right \}\]
Maple ✓
restart; pde :=a*x^k*diff(w(x,y),x) +b*exp(lambda*y)*diff(w(x,y),y) =c*x^n+s; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
\[w \left (x , y\right ) = \frac {c \,x^{-k +n +1}}{\left (-k +n +1\right ) a}-\frac {s \,x^{-k +1}}{\left (k -1\right ) a}+\textit {\_F1} \left (\frac {b \lambda \,x^{-k +1}-\left (k -1\right ) a \,{\mathrm e}^{-\lambda y}}{\left (k -1\right ) b \lambda }\right )\]
____________________________________________________________________________________
Added Feb. 9, 2019.
Problem Chapter 3.3.2.9 from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ a y^k w_x + b e^{\lambda x} w_y = c e^{\mu x}+s \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*y^k*D[w[x, y], x] + b*Exp[lambda*x]*D[w[x, y], y] == c*Exp[mu*x] + s; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to \frac {y^{k+1} \left (\left (y^{k+1}\right )^{\frac {1}{k+1}}\right )^{-k} \left ((k+1) \mu s \, _2F_1\left (1,\frac {1}{k+1};\frac {k+2}{k+1};\frac {a \lambda y^{k+1}}{a \lambda y^{k+1}-b e^{\lambda x} (k+1)}\right )-c \lambda e^{\mu x} \, _2F_1\left (1,\frac {\lambda +k \mu +\mu }{k \lambda +\lambda };\frac {\lambda +\mu }{\lambda };\frac {b e^{\lambda x} (k+1)}{b e^{\lambda x} (k+1)-a \lambda y^{k+1}}\right )\right )}{\mu \left (b (k+1) e^{\lambda x}-a \lambda y^{k+1}\right )}+c_1\left (\frac {y^{k+1}}{k+1}-\frac {b e^{\lambda x}}{a \lambda }\right )\right \}\right \}\]
Maple ✓
restart; pde :=a*y^k*diff(w(x,y),x) +b*exp(lambda*x)*diff(w(x,y),y) =c*exp(mu*x)+s; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
\[w \left (x , y\right ) = \int _{}^{x}\frac {\left (c \,{\mathrm e}^{\textit {\_a} \mu }+s \right ) \left (\left (\frac {a \lambda \,y^{k +1}+\left (k +1\right ) b \,{\mathrm e}^{\textit {\_a} \lambda }-\left (k +1\right ) b \,{\mathrm e}^{\lambda x}}{a \lambda }\right )^{\frac {1}{k +1}}\right )^{-k}}{a}d \textit {\_a} +\textit {\_F1} \left (\frac {a \lambda y \,y^{k}-\left (k +1\right ) b \,{\mathrm e}^{\lambda x}}{a \lambda }\right )\]
____________________________________________________________________________________
Added Feb. 9, 2019.
Problem Chapter 3.3.2.10 from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ a e^{\lambda x} w_x + b y^k w_y = c x^n+s \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*Exp[lambda*x]*D[w[x, y], x] + b*y^k*D[w[x, y], y] == c*x^n + s; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to \frac {-c x^n (\lambda x)^{-n} \operatorname {Gamma}(n+1,\lambda x)+a \lambda c_1\left (\frac {b e^{-\lambda x}}{a \lambda }-\frac {y^{1-k}}{k-1}\right )-s e^{-\lambda x}}{a \lambda }\right \}\right \}\]
Maple ✓
restart; pde :=a*exp(lambda*x)*diff(w(x,y),x) +b*y^k*diff(w(x,y),y) = c*x^n+s; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
\[w \left (x , y\right ) = \frac {c \,x^{n} \left (\lambda x \right )^{-\frac {n}{2}} \WhittakerM \left (\frac {n}{2}, \frac {n}{2}+\frac {1}{2}, \lambda x \right ) {\mathrm e}^{-\frac {\lambda x}{2}}+\left (n +1\right ) a \lambda \textit {\_F1} \left (\frac {a \lambda \,y^{-k +1}-\left (k -1\right ) b \,{\mathrm e}^{-\lambda x}}{a \lambda }\right )-\left (n +1\right ) \left ({\mathrm e}^{-\lambda x}-1\right ) s}{\left (n +1\right ) a \lambda }\]
____________________________________________________________________________________
Added Feb. 9, 2019.
Problem Chapter 3.3.2.11 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^k w_y = c Exp[\mu x]+s \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*Exp[lambda*y]*D[w[x, y], x] + b*x^k*D[w[x, y], y] == c*Exp[mu*x] + s; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to \int _1^x\frac {(k+1) \left (e^{\mu K[1]} c+s\right )}{a e^{\lambda y} (k+1)+b \lambda \left (K[1]^{k+1}-x^{k+1}\right )}dK[1]+c_1\left (\frac {e^{\lambda y}}{\lambda }-\frac {b x^{k+1}}{a k+a}\right )\right \}\right \}\]
Maple ✓
restart; pde :=a*exp(lambda*y)*diff(w(x,y),x) +b*x^k*diff(w(x,y),y) = c*exp(mu*x)+s; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
\[w \left (x , y\right ) = \int _{}^{x}\frac {\left (c \,{\mathrm e}^{\textit {\_a} \mu }+s \right ) \left (k +1\right )}{b \lambda \,\textit {\_a}^{k +1}-b \lambda \,x^{k +1}+\left (k +1\right ) a \,{\mathrm e}^{\lambda y}}d \textit {\_a} +\textit {\_F1} \left (\frac {-b \lambda \,x^{k +1}+\left (k +1\right ) a \,{\mathrm e}^{\lambda y}}{\left (k +1\right ) b \lambda }\right )\]
____________________________________________________________________________________