6.8.1 2.1
6.8.1.1 [1750] Problem 1
problem number 1750
Added June 27, 2019.
Problem Chapter 8.2.2.1, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y,z)\)
\[ a w_x + b w_y + c w_z = (\alpha x+\beta y+\gamma z+ \delta ) w \]
Mathematica ✓
ClearAll["Global`*"];
pde = a*D[w[x, y,z], x] + b*D[w[x, y,z], y] + c*D[w[x,y,z],z]== (alpha*x+beta*y+gamma*z+delta)*w[x,y,z];
sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x, y,z}], 60*10]];
\[\left \{\left \{w(x,y,z)\to c_1\left (y-\frac {b x}{a},z-\frac {c x}{a}\right ) \exp \left (\frac {x (a (\alpha x+2 \beta y+2 \delta +2 \gamma z)-x (b \beta +c \gamma ))}{2 a^2}\right )\right \}\right \}\]
Maple ✓
restart;
local gamma;
pde := a*diff(w(x,y,z),x)+ b*diff(w(x,y,z),y)+ c*diff(w(x,y,z),z)= (alpha*x+beta*y+gamma*z+delta)*w(x,y,z);
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
\[w \left (x , y , z\right ) = f_{1} \left (\frac {y a -b x}{a}, \frac {z a -c x}{a}\right ) {\mathrm e}^{\frac {x \left (\frac {\left (a \alpha -b \beta -c \gamma \right ) x}{2}+a \left (\beta y +\gamma z +\delta \right )\right )}{a^{2}}}\]
___________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________
6.8.1.2 [1751] Problem 2
problem number 1751
Added June 27, 2019.
Problem Chapter 8.2.2.2, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y,z)\)
\[ w_x + a z w_y + b y w_z = (c x+s) w \]
Mathematica ✓
ClearAll["Global`*"];
pde = D[w[x, y,z], x] + a*z*D[w[x, y,z], y] + b*y*D[w[x,y,z],z]== (c*x+s)*w[x,y,z];
sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x, y,z}], 60*10]];
\[\left \{\left \{w(x,y,z)\to e^{\frac {c x^2}{2}+s x} c_1\left (\frac {e^{-\sqrt {a} \sqrt {b} x} \left (\sqrt {b} y \left (e^{2 \sqrt {a} \sqrt {b} x}+1\right )-\sqrt {a} z \left (e^{2 \sqrt {a} \sqrt {b} x}-1\right )\right )}{2 \sqrt {b}},\frac {e^{-\sqrt {a} \sqrt {b} x} \left (\sqrt {a} z \left (e^{2 \sqrt {a} \sqrt {b} x}+1\right )-\sqrt {b} y \left (e^{2 \sqrt {a} \sqrt {b} x}-1\right )\right )}{2 \sqrt {a}}\right )\right \}\right \}\]
Maple ✓
restart;
local gamma;
pde := diff(w(x,y,z),x)+ a*z*diff(w(x,y,z),y)+ b*y*diff(w(x,y,z),z)= (c*x+s)*w(x,y,z);
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
\[w \left (x , y , z\right ) = f_{1} \left (\frac {a \,z^{2}-y^{2} b}{a}, \frac {x \sqrt {a b}-\ln \left (\frac {a b y +\sqrt {a^{2} z^{2}}\, \sqrt {a b}}{\sqrt {a b}}\right )}{\sqrt {a b}}\right ) {\mathrm e}^{\frac {\int _{}^{y}\frac {\left (c x +s \right ) \sqrt {a b}+\ln \left (\frac {a b \textit {\_a} +\sqrt {a \left (\left (\textit {\_a}^{2}-y^{2}\right ) b +a \,z^{2}\right )}\, \sqrt {a b}}{\sqrt {a b}}\right ) c -\ln \left (\frac {a b y +\sqrt {a^{2} z^{2}}\, \sqrt {a b}}{\sqrt {a b}}\right ) c}{\sqrt {a \left (\left (\textit {\_a}^{2}-y^{2}\right ) b +a \,z^{2}\right )}}d \textit {\_a}}{\sqrt {a b}}}\]
___________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________
6.8.1.3 [1752] Problem 3
problem number 1752
Added June 27, 2019.
Problem Chapter 8.2.2.3, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y,z)\)
\[ w_x + (a_1 x+a_0) w_y + (b_1 x+b_0) w_z = (\alpha x+\beta y+\gamma z+\delta ) w \]
Mathematica ✓
ClearAll["Global`*"];
pde = D[w[x, y,z], x] + (a1*x+a0)*D[w[x, y,z], y] + (b1*x+b0)*D[w[x,y,z],z]== (alpha*x+beta*y+gamma*z+delta)*w[x,y,z];
sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x, y,z}], 60*10]];
\[\left \{\left \{w(x,y,z)\to c_1\left (-\text {a0} x-\frac {\text {a1} x^2}{2}+y,-\text {b0} x-\frac {\text {b1} x^2}{2}+z\right ) \exp \left (\frac {1}{6} x \left (-3 \text {a0} \beta x-2 \text {a1} \beta x^2+3 \alpha x-3 \text {b0} \gamma x-2 \text {b1} \gamma x^2+6 \beta y+6 \delta +6 \gamma z\right )\right )\right \}\right \}\]
Maple ✓
restart;
local gamma;
pde := diff(w(x,y,z),x)+ (a1*x+a0)*diff(w(x,y,z),y)+ (b1*x+b0)*diff(w(x,y,z),z)= (alpha*x+beta*y+gamma*z+delta)*w(x,y,z);
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
\[w \left (x , y , z\right ) = f_{1} \left (-\frac {1}{2} \operatorname {a1} \,x^{2}-\operatorname {a0} x +y , -\frac {1}{2} \operatorname {b1} \,x^{2}-\operatorname {b0} x +z \right ) {\mathrm e}^{-\frac {x \left (\left (\operatorname {b1} \,x^{2}+\frac {3}{2} \operatorname {b0} x -3 z \right ) \gamma +\left (\operatorname {a1} \,x^{2}+\frac {3}{2} \operatorname {a0} x -3 y \right ) \beta -\frac {3 \alpha x}{2}-3 \delta \right )}{3}}\]
___________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________
6.8.1.4 [1753] Problem 4
problem number 1753
Added June 27, 2019.
Problem Chapter 8.2.2.4, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y,z)\)
\[ w_x + (a_2 y+a_1 x+a_0) w_y + (b_2 y+b_1 x+b_0) w_z = (c_2 y+c_1 z+c_0 x+s) w \]
Mathematica ✓
ClearAll["Global`*"];
pde = D[w[x, y,z], x] + (a2*y+a1*x+a0)*D[w[x, y,z], y] + (b2*y+b1*x+b0)*D[w[x,y,z],z]== (c2*y+c1*z+c0*x+s)*w[x,y,z];
sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x, y,z}], 60*10]];
\[\left \{\left \{w(x,y,z)\to c_1\left (\frac {e^{-\text {a2} x} (\text {a2} (\text {a0}+\text {a2} y)+\text {a1} \text {a2} x+\text {a1})}{\text {a2}^2},\frac {e^{-\text {a2} x} \left (\text {a2} \left (2 \text {a0} \text {b2} \left (\text {a2} x e^{\text {a2} x}+1\right )-\text {a2} \left (\text {a2} e^{\text {a2} x} \left (2 \text {b0} x+\text {b1} x^2-2 z\right )+2 \text {b2} y \left (e^{\text {a2} x}-1\right )\right )\right )+\text {a1} \text {b2} \left (\text {a2}^2 x^2 e^{\text {a2} x}+2 \text {a2} x+2\right )\right )}{2 \text {a2}^3}\right ) \exp \left (\frac {\text {a2} \left (3 \text {a0} \text {b2} \text {c1} \left (\text {a2}^2 x^2-2 \text {a2} x+2\right )-6 \text {a0} \text {a2} \text {c2} (\text {a2} x-1)+\text {a2}^3 x \left (-3 \text {b0} \text {c1} x-2 \text {b1} \text {c1} x^2+3 \text {c0} x+6 \text {c1} z+6 s\right )+6 \text {a2}^2 y (\text {c2}-\text {b2} \text {c1} x)+6 \text {a2} \text {b2} \text {c1} y\right )+\text {a1} \left (\text {b2} \text {c1} \left (2 \text {a2}^3 x^3-3 \text {a2}^2 x^2+6\right )-3 \text {a2} \text {c2} \left (\text {a2}^2 x^2-2\right )\right )}{6 \text {a2}^4}\right )\right \}\right \}\]
Maple ✓
restart;
local gamma;
pde := diff(w(x,y,z),x)+ (a2*y+a1*x+a0)*diff(w(x,y,z),y)+ (b2*y+b1*x+b0)*diff(w(x,y,z),z)= (c2*y+c1*z+c0*x+s)*w(x,y,z);
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
\[w \left (x , y , z\right ) = f_{1} \left (\frac {{\mathrm e}^{-\operatorname {a2} x} \left (y \,\operatorname {a2}^{2}+\operatorname {a2} \left (\operatorname {a1} x +\operatorname {a0} \right )+\operatorname {a1} \right )}{\operatorname {a2}^{2}}, \frac {\left (-\operatorname {b1} \,x^{2}-2 \operatorname {b0} x +2 z \right ) \operatorname {a2}^{3}+\operatorname {b2} \left (\operatorname {a1} \,x^{2}+2 \operatorname {a0} x -2 y \right ) \operatorname {a2}^{2}-2 \operatorname {a0} \operatorname {a2} \operatorname {b2} -2 \operatorname {a1} \operatorname {b2}}{2 \operatorname {a2}^{3}}\right ) {\mathrm e}^{\frac {\left (\left (2 x^{3} \operatorname {a1} \,\operatorname {a2}^{3}+\left (3 \operatorname {a0} \,\operatorname {a2}^{3}-3 \operatorname {a1} \,\operatorname {a2}^{2}\right ) x^{2}+\left (-6 \operatorname {a2}^{3} y -6 \operatorname {a0} \,\operatorname {a2}^{2}\right ) x +6 y \,\operatorname {a2}^{2}+6 \operatorname {a0} \operatorname {a2} +6 \operatorname {a1} \right ) \operatorname {b2} -2 x \,\operatorname {a2}^{4} \left (\operatorname {b1} \,x^{2}+\frac {3}{2} \operatorname {b0} x -3 z \right )\right ) \operatorname {c1} +6 \operatorname {a2} \left (\frac {\operatorname {a2}^{2} \left (-\operatorname {a1} \operatorname {c2} +\operatorname {a2} \operatorname {c0} \right ) x^{2}}{2}+\left (-\operatorname {a0} \,\operatorname {a2}^{2} \operatorname {c2} +\operatorname {a2}^{3} s \right ) x +\operatorname {c2} \left (y \,\operatorname {a2}^{2}+\operatorname {a0} \operatorname {a2} +\operatorname {a1} \right )\right )}{6 \operatorname {a2}^{4}}}\]
___________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________
6.8.1.5 [1754] Problem 5
problem number 1754
Added June 27, 2019.
Problem Chapter 8.2.2.5, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y,z)\)
\[ w_x + (a y+k_1 x+k_0) w_y + (b z+s_1 x+s_0) w_z = (c_1 x+c_0) w \]
Mathematica ✓
ClearAll["Global`*"];
pde = D[w[x, y,z], x] + (a*y+k1*x+k0)*D[w[x, y,z], y] + (b*z+s1*x+s0)*D[w[x,y,z],z]== (c1*x+c0)*w[x,y,z];
sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x, y,z}], 60*10]];
\[\left \{\left \{w(x,y,z)\to e^{\text {c0} x+\frac {\text {c1} x^2}{2}} c_1\left (\frac {e^{-a x} \left (a^2 y+a (\text {k0}+\text {k1} x)+\text {k1}\right )}{a^2},\frac {e^{-b x} \left (b^2 z+b (\text {s0}+\text {s1} x)+\text {s1}\right )}{b^2}\right )\right \}\right \}\]
Maple ✓
restart;
local gamma;
pde := diff(w(x,y,z),x)+(a*y+k1*x+k0)*diff(w(x,y,z),y)+(b*z+s1*x+s0)*diff(w(x,y,z),z)= (c1*x+c0)*w(x,y,z);
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
\[w \left (x , y , z\right ) = f_{1} \left (\frac {\left (y \,a^{2}+a \left (\operatorname {k1} x +\operatorname {k0} \right )+\operatorname {k1} \right ) {\mathrm e}^{-a x}}{a^{2}}, \frac {\left (z \,b^{2}+b \left (\operatorname {s1} x +\operatorname {s0} \right )+\operatorname {s1} \right ) {\mathrm e}^{-b x}}{b^{2}}\right ) {\mathrm e}^{\frac {x \left (\operatorname {c1} x +2 \operatorname {c0} \right )}{2}}\]
___________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________
6.8.1.6 [1755] Problem 6
problem number 1755
Added June 27, 2019.
Problem Chapter 8.2.2.6, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y,z)\)
\[ a x w_x + b y w_y + c z w_z = (\alpha x+\beta y+\gamma z+\delta ) w \]
Mathematica ✓
ClearAll["Global`*"];
pde = a*x*D[w[x, y,z], x] + b*y*D[w[x, y,z], y] + c*z*D[w[x,y,z],z]== (alpha*x+beta*y+gamma*z+beta)*w[x,y,z];
sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x, y,z}], 60*10]];
\[\left \{\left \{w(x,y,z)\to x^{\frac {\beta }{a}} c_1\left (y x^{-\frac {b}{a}},z x^{-\frac {c}{a}}\right ) e^{\frac {\alpha x}{a}+\frac {\beta y}{b}+\frac {\gamma z}{c}}\right \}\right \}\]
Maple ✓
restart;
local gamma;
pde := a*x*diff(w(x,y,z),x)+ b*y*diff(w(x,y,z),y)+c*z*diff(w(x,y,z),z)= (alpha*x+beta*y+gamma*z+beta)*w(x,y,z);
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
\[w \left (x , y , z\right ) = f_{1} \left (y \,x^{-\frac {b}{a}}, z \,x^{-\frac {c}{a}}\right ) x^{\frac {\beta }{a}} {\mathrm e}^{\frac {\left (a \beta y +\alpha b x \right ) c +\gamma z a b}{a b c}}\]
___________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________
6.8.1.7 [1756] Problem 7
problem number 1756
Added June 27, 2019.
Problem Chapter 8.2.2.7, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y,z)\)
\[ x w_x + a z w_y + b y w_z = c w \]
Mathematica ✓
ClearAll["Global`*"];
pde = x*D[w[x, y,z], x] + a*z*D[w[x, y,z], y] + b*y*D[w[x,y,z],z]== c*w[x,y,z];
sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x, y,z}], 60*10]];
\[\left \{\left \{w(x,y,z)\to x^c c_1\left (i y \sinh \left (\sqrt {a} \sqrt {b} \log (x)\right )-\frac {i \sqrt {a} z \cosh \left (\sqrt {a} \sqrt {b} \log (x)\right )}{\sqrt {b}},y \cosh \left (\sqrt {a} \sqrt {b} \log (x)\right )-\frac {\sqrt {a} z \sinh \left (\sqrt {a} \sqrt {b} \log (x)\right )}{\sqrt {b}}\right )\right \}\right \}\]
Maple ✓
restart;
local gamma;
pde := x*diff(w(x,y,z),x)+ a*z*y*diff(w(x,y,z),y)+b*y*diff(w(x,y,z),z)= c*w(x,y,z);
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
\[w \left (x , y , z\right ) = f_{1} \left (\frac {z^{2} a -2 b y}{a}, x \,{\mathrm e}^{\frac {2 \,\operatorname {arctanh}\left (\frac {z}{\sqrt {\frac {z^{2} a -2 b y}{a}}}\right )}{a \sqrt {\frac {z^{2} a -2 b y}{a}}}}\right ) {\mathrm e}^{-\frac {2 c \,\operatorname {arctanh}\left (\frac {\sqrt {a^{2} z^{2}}}{\sqrt {\frac {z^{2} a -2 b y}{a}}\, a}\right )}{\sqrt {\frac {z^{2} a -2 b y}{a}}\, a}}\]
___________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________
6.8.1.8 [1757] Problem 8
problem number 1757
Added June 27, 2019.
Problem Chapter 8.2.2.8, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y,z)\)
\[ a b x w_x + b(a y + b z) w_y + a(a y -b z) w_z = c w \]
Mathematica ✓
ClearAll["Global`*"];
pde = a*b*x*D[w[x, y,z], x] + b*(a*y+b*z)*D[w[x, y,z], y] + a*(a*y-b*z)*D[w[x,y,z],z]== c*w[x,y,z];
sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x, y,z}], 60*10]];
\[\left \{\left \{w(x,y,z)\to x^{\frac {c}{a b}} c_1\left (-\frac {x^{\sqrt {2}} \left (\left (\sqrt {2}-2\right ) a y+\sqrt {2} b z\right )}{4 a},\frac {x^{-\sqrt {2}} \left (\left (2+\sqrt {2}\right ) a y+\sqrt {2} b z\right )}{4 a}\right )\right \}\right \}\]
Maple ✓
restart;
local gamma;
pde := a*b*x*diff(w(x,y,z),x)+ b*(a*y+b*z)*diff(w(x,y,z),y)+a*(a*y-b*z)*diff(w(x,y,z),z)= c*w(x,y,z);
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
\[w \left (x , y , z\right ) = f_{1} \left (-\frac {1}{\sqrt {-a^{2} y^{2}+2 a b y z +b^{2} z^{2}}}, x {\left (-\frac {\left (a^{3} y^{3}-a^{2} b \,y^{2} z -3 a \,b^{2} y \,z^{2}-b^{3} z^{3}\right ) \sqrt {-\frac {a^{2}}{a^{2} y^{2}-2 a b y z -b^{2} z^{2}}}-a^{2} y \sqrt {-2 a^{2} y^{2}+4 a b y z +2 b^{2} z^{2}}}{\left (-a^{2} y^{2}+2 a b y z +b^{2} z^{2}\right )^{{3}/{2}} \sqrt {-\frac {a^{2}}{a^{2} y^{2}-2 a b y z -b^{2} z^{2}}}}\right )}^{-\frac {a}{\sqrt {-2 a^{2} y^{2}+4 a b y z +2 b^{2} z^{2}}\, \sqrt {-\frac {a^{2}}{a^{2} y^{2}-2 a b y z -b^{2} z^{2}}}}}\right ) {\left (-\frac {\left (a^{3} y^{3}-a^{2} b \,y^{2} z -3 a \,b^{2} y \,z^{2}-b^{3} z^{3}\right ) \sqrt {-\frac {a^{2}}{a^{2} y^{2}-2 a b y z -b^{2} z^{2}}}-a^{2} y \sqrt {-2 a^{2} y^{2}+4 a b y z +2 b^{2} z^{2}}}{\left (-a^{2} y^{2}+2 a b y z +b^{2} z^{2}\right )^{{3}/{2}} \sqrt {-\frac {a^{2}}{a^{2} y^{2}-2 a b y z -b^{2} z^{2}}}}\right )}^{\frac {c}{\sqrt {-2 a^{2} y^{2}+4 a b y z +2 b^{2} z^{2}}\, b \sqrt {-\frac {a^{2}}{a^{2} y^{2}-2 a b y z -b^{2} z^{2}}}}}\]
___________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________
6.8.1.9 [1758] Problem 9
problem number 1758
Added June 27, 2019.
Problem Chapter 8.2.2.9, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y,z)\)
\[ (a_1 x +a_0) w_x + (b_1 y + b_0) w_y + (c_1 z +c_0) w_z = (\alpha x+\beta y+\gamma z+\delta ) w \]
Mathematica ✓
ClearAll["Global`*"];
pde = (a1*x+a0)*D[w[x, y,z], x] + (b1*y+b0)*D[w[x, y,z], y] +(c1*z+c0)*D[w[x,y,z],z]== (alpha*x+beta*y+gamma*z+delta)*w[x,y,z];
sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x, y,z}], 60*10]];
\[\left \{\left \{w(x,y,z)\to (\text {a0}+\text {a1} x)^{-\frac {\text {a0} \alpha \text {b1} \text {c1}+\text {a1} \text {b0} \beta \text {c1}+\text {a1} \text {b1} \text {c0} \gamma -\text {a1} \text {b1} \text {c1} \delta }{\text {a1}^2 \text {b1} \text {c1}}} c_1\left (\frac {(\text {b0}+\text {b1} y) (\text {a0}+\text {a1} x)^{-\frac {\text {b1}}{\text {a1}}}}{\text {b1}},\frac {(\text {c0}+\text {c1} z) (\text {a0}+\text {a1} x)^{-\frac {\text {c1}}{\text {a1}}}}{\text {c1}}\right ) \exp \left (\frac {\alpha x}{\text {a1}}+\frac {\beta (\text {b0}+\text {b1} y)}{\text {b1}^2}+\frac {\gamma (\text {c0}+\text {c1} z)}{\text {c1}^2}\right )\right \}\right \}\]
Maple ✓
restart;
local gamma;
pde := (a1*x+a0)*diff(w(x,y,z),x)+(b1*y+b0)*diff(w(x,y,z),y)+(c1*z+c0)*diff(w(x,y,z),z)= (alpha*x+beta*y+gamma*z+delta)*w(x,y,z);
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
\[w \left (x , y , z\right ) = f_{1} \left (\frac {\left (\operatorname {b1} y +\operatorname {b0} \right ) \left (\operatorname {a1} x +\operatorname {a0} \right )^{-\frac {\operatorname {b1}}{\operatorname {a1}}}}{\operatorname {b1}}, \frac {\left (\operatorname {c1} z +\operatorname {c0} \right ) \left (\operatorname {a1} x +\operatorname {a0} \right )^{-\frac {\operatorname {c1}}{\operatorname {a1}}}}{\operatorname {c1}}\right ) \left (\operatorname {a1} x +\operatorname {a0} \right )^{\frac {\left (\left (-\operatorname {c0} \gamma +\operatorname {c1} \delta \right ) \operatorname {b1} -\operatorname {b0} \beta \operatorname {c1} \right ) \operatorname {a1} -\operatorname {a0} \alpha \operatorname {b1} \operatorname {c1}}{\operatorname {a1}^{2} \operatorname {b1} \operatorname {c1}}} {\mathrm e}^{\frac {\alpha x}{\operatorname {a1}}+\frac {\beta \left (\operatorname {b1} y +\operatorname {b0} \right )}{\operatorname {b1}^{2}}+\frac {\left (\operatorname {c1} z +\operatorname {c0} \right ) \gamma }{\operatorname {c1}^{2}}}\]
___________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________