Added Feb. 9, 2019.
Problem Chapter 3.2.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^2+d y^2+ k x y+n \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*D[w[x, y], x] + b*D[w[x, y], y] == c*x^2 + d*y^2 + k*x*y + n; 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 {x \left (a^2 \left (2 c x^2+6 d y^2+3 k x y+6 n\right )-a b x (6 d y+k x)+2 b^2 d x^2\right )}{6 a^3}\right \}\right \}\]
Maple ✓
restart; pde := a* diff(w(x,y),x)+b*diff(w(x,y),y) = c*x^2+d*y^2+k*x*y+n; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
\[w \left ( x,y \right ) =1/6\,{\frac {1}{{a}^{3}} \left ( 6\,{\it \_F1} \left ( {\frac {ya-bx}{a}} \right ) {a}^{3}+2\, \left ( \left ( c{x}^{2}+3\,d{y}^{2}+3/2\,kxy+3\,n \right ) {a}^{2}-3\,bx \left ( dy+1/6\,kx \right ) a+{b}^{2}d{x}^{2} \right ) x \right ) }\]
____________________________________________________________________________________
Added Feb. 9, 2019.
Problem Chapter 3.2.2.2 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 = c x^2+d y^2+ k x y+n \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*x*D[w[x, y], x] + b*y*D[w[x, y], y] == c*x^2 + d*y^2 + k*x*y + n; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to \frac {2 a b (a+b) c_1\left (y x^{-\frac {b}{a}}\right )+a^2 d y^2+a b c x^2+a b d y^2+2 a b k x y+2 b n (a+b) \log (x)+b^2 c x^2}{2 a b (a+b)}\right \}\right \}\]
Maple ✓
restart; pde := a*x*diff(w(x,y),x)+b*y*diff(w(x,y),y) = c*x^2+d*y^2+k*x*y+n; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
\[w \left ( x,y \right ) ={\frac {kxy}{a+b}}+1/2\,{\frac {c{x}^{2}}{a}}+{\frac {n\ln \left ( x \right ) }{a}}+1/2\,{\frac {d{y}^{2}}{b}}+{\it \_F1} \left ( y{x}^{-{\frac {b}{a}}} \right ) \]
____________________________________________________________________________________
Added Feb. 9, 2019.
Problem Chapter 3.2.2.3 from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ a y w_x + b x w_y = c x y+d \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*y*D[w[x, y], x] + b*x*D[w[x, y], y] == c*x*y + d; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\begin {align*} & \left \{w(x,y)\to c_1\left (\frac {a y^2-b x^2}{2 a}\right )-\frac {d \tanh ^{-1}\left (\frac {\sqrt {b} x}{\sqrt {a y^2}}\right )}{\sqrt {a} \sqrt {b}}+\frac {c x^2}{2 a}\right \}\\& \left \{w(x,y)\to c_1\left (\frac {a y^2-b x^2}{2 a}\right )+\frac {d \tanh ^{-1}\left (\frac {\sqrt {b} x}{\sqrt {a y^2}}\right )}{\sqrt {a} \sqrt {b}}+\frac {c x^2}{2 a}\right \}\\ \end {align*}
Maple ✓
restart; pde := a*y*diff(w(x,y),x)+b*x*diff(w(x,y),y) = c*x*y+d; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
\[w \left ( x,y \right ) ={\frac {1}{a\sqrt {ab}} \left ( 1/2\,c{x}^{2}\sqrt {ab}+d\ln \left ( {\frac {axb}{\sqrt {ab}}}+\sqrt {{y}^{2}{a}^{2}} \right ) a+{\it \_F1} \left ( {\frac {{y}^{2}a-b{x}^{2}}{a}} \right ) a\sqrt {ab} \right ) }\]
____________________________________________________________________________________
Added Feb. 9, 2019.
Problem Chapter 3.2.2.4 from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ a x^2 w_x + b y^2 w_y = c x^2+d y^2+ k x y+ n x+ m y+s \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*x^2*D[w[x, y], x] + b*y^2*D[w[x, y], y] == c*x^2 + d*y^2 + k*x*y + n*x + m*y + s; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to \frac {a b x (a x-b y) c_1\left (\frac {b}{a x}-\frac {1}{y}\right )-a^2 m x^2 \log \left (\frac {a x}{y}\right )+a b c x^3-a b d x y^2+a b k x^2 y \log \left (\frac {a x}{y}\right )+x \log (x) (a m+b n) (a x-b y)+a b m x y \log \left (\frac {a x}{y}\right )-a b s x-b^2 c x^2 y+b^2 s y}{a b x (a x-b y)}\right \}\right \}\]
Maple ✓
restart; pde := a*x^2*diff(w(x,y),x)+b*y^2*diff(w(x,y),y) =c*x^2+d*y^2+ k*x*y+ n*x+ m*y+s; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
\[w \left ( x,y \right ) ={\frac {cx}{a}}-{\frac {d{y}^{2}}{ax-by}}-{\frac {m}{b}\ln \left ( {\frac {ax}{y}} \right ) }+{\frac {kxy}{ax-by}\ln \left ( {\frac {ax}{y}} \right ) }-{\frac {s}{ax}}+{\frac {m\ln \left ( x \right ) }{b}}+{\frac {n\ln \left ( x \right ) }{a}}+{\it \_F1} \left ( {\frac {ax-by}{axy}} \right ) \]
____________________________________________________________________________________
Added Feb. 9, 2019.
Problem Chapter 3.2.2.5 from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ x^2 w_x + a x y w_y = b y^2 \]
Mathematica ✓
ClearAll["Global`*"]; pde = x^2*D[w[x, y], x] + a*x*y*D[w[x, y], y] == b*y^2; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to c_1\left (y x^{-a}\right )-\frac {b y^2}{x-2 a x}\right \}\right \}\]
Maple ✓
restart; pde := x^2*diff(w(x,y),x)+a*x*y*diff(w(x,y),y) =b*y^2; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
\[w \left ( x,y \right ) ={\frac {b{y}^{2}}{ \left ( 2\,a-1 \right ) x}}+{\it \_F1} \left ( y{x}^{-a} \right ) \]
____________________________________________________________________________________
Added Feb. 9, 2019.
Problem Chapter 3.2.2.6 from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ a y^2 w_x + b x^2 w_y = c x^2+d \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*y^2*D[w[x, y], x] + b*x^2*D[w[x, y], y] == c*x^2 + d; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\begin {align*} & \left \{w(x,y)\to c_1\left (\frac {a y^3-b x^3}{3 a}\right )+\frac {b d x \left (\frac {a y^3}{a y^3-b x^3}\right )^{2/3} \, _2F_1\left (\frac {1}{3},\frac {2}{3};\frac {4}{3};\frac {b x^3}{b x^3-a y^3}\right )+a c y^3}{\sqrt [3]{a} b \left (a y^3\right )^{2/3}}\right \}\\& \left \{w(x,y)\to c_1\left (\frac {a y^3-b x^3}{3 a}\right )-\frac {\sqrt [3]{-1} \left (b d x \left (\frac {a y^3}{a y^3-b x^3}\right )^{2/3} \, _2F_1\left (\frac {1}{3},\frac {2}{3};\frac {4}{3};\frac {b x^3}{b x^3-a y^3}\right )+a c y^3\right )}{\sqrt [3]{a} b \left (a y^3\right )^{2/3}}\right \}\\& \left \{w(x,y)\to c_1\left (\frac {a y^3-b x^3}{3 a}\right )+\frac {(-1)^{2/3} \left (b d x \left (\frac {a y^3}{a y^3-b x^3}\right )^{2/3} \, _2F_1\left (\frac {1}{3},\frac {2}{3};\frac {4}{3};\frac {b x^3}{b x^3-a y^3}\right )+a c y^3\right )}{\sqrt [3]{a} b \left (a y^3\right )^{2/3}}\right \}\\ \end {align*}
Maple ✓
restart; pde := a*y^2*diff(w(x,y),x)+b*x^2*diff(w(x,y),y) =c*x^2+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 { \left ( {{\it \_a}}^{2}c+d \right ) a}{ \left ( \left ( {{\it \_a}}^{3}b+\RootOf \left ( ya-\sqrt [3]{{a}^{2}b{x}^{3}+{a}^{3}{\it \_Z}} \right ) a \right ) {a}^{2} \right ) ^{2/3}}}{d{\it \_a}}+{\it \_F1} \left ( \RootOf \left ( ya-\sqrt [3]{{a}^{2}b{x}^{3}+{a}^{3}{\it \_Z}} \right ) \right ) \] Contains unresolved integral with RootOf
____________________________________________________________________________________
Added Feb. 9, 2019.
Problem Chapter 3.2.2.7 from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y)\)
\[ a y^2 w_x + b x y w_y = c x^2+d y^2 \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*y^2*D[w[x, y], x] + b*x*y*D[w[x, y], y] == c*x^2 + d*y^2; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
\[\left \{\left \{w(x,y)\to c_1\left (\frac {a y^2-b x^2}{2 a}\right )-\frac {\sqrt {a} c \sqrt {y^2-\frac {b x^2}{a}} \tan ^{-1}\left (\frac {\sqrt {b} x}{\sqrt {a} \sqrt {y^2-\frac {b x^2}{a}}}\right )}{b^{3/2}}+\frac {d x}{a}+\frac {c x}{b}\right \}\right \}\]
Maple ✓
restart; pde := a*y^2*diff(w(x,y),x)+b*x*y*diff(w(x,y),y) =c*x^2+d*y^2; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y)) ),output='realtime'));
\[w \left ( x,y \right ) ={\frac {1}{ba\sqrt { \left ( {y}^{2}a-b{x}^{2} \right ) b}} \left ( -ac \left ( {y}^{2}a-b{x}^{2} \right ) \arctan \left ( {\frac {bx}{\sqrt { \left ( {y}^{2}a-b{x}^{2} \right ) b}}} \right ) +\sqrt { \left ( {y}^{2}a-b{x}^{2} \right ) b} \left ( {\it \_F1} \left ( {\frac {{y}^{2}a-b{x}^{2}}{a}} \right ) ab+ \left ( ca+bd \right ) x \right ) \right ) }\]
____________________________________________________________________________________