6.7.23 7.4

6.7.23.1 [1717] Problem 1
6.7.23.2 [1718] Problem 2
6.7.23.3 [1719] Problem 3
6.7.23.4 [1720] Problem 4
6.7.23.5 [1721] Problem 5

6.7.23.1 [1717] Problem 1

problem number 1717

Added June 26, 2019.

Problem Chapter 7.7.4.1, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.

Solve for \(w(x,y,z)\)

\[ w_x + a w_y + b w_z = c \arccot ^k(\lambda x)+s \]

Mathematica

ClearAll["Global`*"]; 
pde =  D[w[x, y,z], x] + a*D[w[x, y,z], y] + b*D[w[x,y,z],z]==c*ArcCot[lambda*x]^k+s; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x, y,z}], 60*10]];
 

\[\left \{\left \{w(x,y,z)\to \int _1^x\left (c \cot ^{-1}(\lambda K[1])^k+s\right )dK[1]+c_1(y-a x,z-b x)\right \}\right \}\]

Maple

restart; 
local gamma; 
pde :=  diff(w(x,y,z),x)+ a*diff(w(x,y,z),y)+ b*diff(w(x,y,z),z)= c*arccot(lambda*x)^k+s; 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
 

\[w \left (x , y , z\right ) = s x +\int c \left (-\arctan \left (\lambda x \right )+\frac {\pi }{2}\right )^{k}d x +\mathit {\_F1} \left (-a x +y , -b x +z \right )\]

____________________________________________________________________________________

6.7.23.2 [1718] Problem 2

problem number 1718

Added June 26, 2019.

Problem Chapter 7.7.4.2, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.

Solve for \(w(x,y,z)\)

\[ a_1 w_x + a_2 w_y + a_3 w_z = b_1 \arccot (\lambda _1 x)+b_2 \arccot (\lambda _2 y)+b_3 \arccot (\lambda _3 z) \]

Mathematica

ClearAll["Global`*"]; 
pde =  a1*D[w[x, y,z], x] + a2*D[w[x, y,z], y] + a3*D[w[x,y,z],z]== b1*ArcCot[lambda1*x]+b2*ArcCot[lambda2*y]+b3*ArcCot[lambda3*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 {\text {a2} x}{\text {a1}},z-\frac {\text {a3} x}{\text {a1}}\right )+\frac {1}{2} \left (\frac {\text {b2} \log \left (\text {a1}^2 \left (\text {lambda2}^2 y^2+1\right )\right )}{\text {a2} \text {lambda2}}+\frac {\text {b3} \log \left (\text {a1}^2 \left (\text {lambda3}^2 z^2+1\right )\right )}{\text {a3} \text {lambda3}}+\frac {\text {b1} \log \left (\text {lambda1}^2 x^2+1\right )}{\text {a1} \text {lambda1}}+\frac {2 \text {b2} x \tan ^{-1}(\text {lambda2} y)}{\text {a1}}+\frac {2 \text {b2} x \cot ^{-1}(\text {lambda2} y)}{\text {a1}}+\frac {2 \text {b3} x \tan ^{-1}(\text {lambda3} z)}{\text {a1}}+\frac {2 \text {b3} x \cot ^{-1}(\text {lambda3} z)}{\text {a1}}-\frac {2 \text {b2} y \tan ^{-1}(\text {lambda2} y)}{\text {a2}}-\frac {2 \text {b3} z \tan ^{-1}(\text {lambda3} z)}{\text {a3}}\right )+\frac {\text {b1} x \cot ^{-1}(\text {lambda1} x)}{\text {a1}}\right \}\right \}\]

Maple

restart; 
local gamma; 
pde :=  a1*diff(w(x,y,z),x)+ a2*diff(w(x,y,z),y)+ a3*diff(w(x,y,z),z)= b1*arccot(lambda1*x)+b2*arccot(lambda2*y)+b3*arccot(lambda3*z); 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
 

\[w \left (x , y , z\right ) = \frac {2 \mathit {a1} \mathit {a2} \mathit {a3} \lambda 1 \lambda 2 \lambda 3 \mathit {\_F1} \left (\frac {\mathit {a1} y -\mathit {a2} x}{\mathit {a1}}, \frac {z \mathit {a1} -\mathit {a3} x}{\mathit {a1}}\right )+\mathit {a2} \mathit {a3} \mathit {b1} \lambda 2 \lambda 3 \ln \left (\lambda 1^{2} x^{2}+1\right )+\left (\mathit {a1} \mathit {a3} \mathit {b2} \lambda 3 \ln \left (\lambda 2^{2} y^{2}+1\right )+\left (\mathit {a1} \mathit {a2} \mathit {b3} \ln \left (\lambda 3^{2} z^{2}+1\right )+\left (-2 \mathit {a1} \mathit {a3} \mathit {b2} y \arctan \left (\lambda 2 y \right )-2 \mathit {a2} \mathit {a3} \mathit {b1} x \arctan \left (\lambda 1 x \right )+\left (-2 \mathit {a1} \mathit {b3} z \arctan \left (\lambda 3 z \right )+\pi \left (\mathit {b1} +\mathit {b2} +\mathit {b3} \right ) \mathit {a3} x \right ) \mathit {a2} \right ) \lambda 3 \right ) \lambda 2 \right ) \lambda 1}{2 \mathit {a1} \mathit {a2} \mathit {a3} \lambda 1 \lambda 2 \lambda 3}\]

____________________________________________________________________________________

6.7.23.3 [1719] Problem 3

problem number 1719

Added June 26, 2019.

Problem Chapter 7.7.4.3, 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 \arccot ^n(\lambda x) \arccot ^k(\beta z) w_z = s \arccot ^m(\gamma x) \]

Mathematica

ClearAll["Global`*"]; 
pde =  a*D[w[x, y,z], x] + b*D[w[x, y,z], y] + c*ArcCot[lambda*x]^n*ArcCot[beta*z]^k*D[w[x,y,z],z]== s*ArcCot[gamma*x]^m; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x, y,z}], 60*10]];
 

Failed

Maple

restart; 
local gamma; 
pde :=  a*diff(w(x,y,z),x)+ b*diff(w(x,y,z),y)+ c*arccot(lambda*x)^n*arccot(beta*z)^k*diff(w(x,y,z),z)= s*arccot(gamma*x)^m; 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
 

\[w \left (x , y , z\right ) = \int \frac {s \left (-\arctan \left (\gamma x \right )+\frac {\pi }{2}\right )^{m}}{a}d x +\mathit {\_F1} \left (\frac {a y -b x}{a}, -\left (\int \left (-\arctan \left (\lambda x \right )+\frac {\pi }{2}\right )^{n}d x \right )+\int \frac {a \left (-\arctan \left (\beta z \right )+\frac {\pi }{2}\right )^{-k}}{c}d z \right )\]

____________________________________________________________________________________

6.7.23.4 [1720] Problem 4

problem number 1720

Added June 26, 2019.

Problem Chapter 7.7.4.4, 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 \arccot ^n(\lambda x) \arccot ^m(\beta y) \arccot ^k(\gamma z) w_z = s \]

Mathematica

ClearAll["Global`*"]; 
pde =  a*D[w[x, y,z], x] + b*D[w[x, y,z], y] + c*ArcCot[lambda*x]^n*ArcCot[beta*y]^m*ArcCot[gamma*z]^k*D[w[x,y,z],z]== s; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x, y,z}], 60*10]];
 

Failed

Maple

restart; 
local gamma; 
pde :=  a*diff(w(x,y,z),x)+ b*diff(w(x,y,z),y)+ c*arccot(lambda*x)^n*arccot(beta*y)^m*arccot(gamma*z)^k*diff(w(x,y,z),z)= s; 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
 

\[w \left (x , y , z\right ) = \frac {s x}{a}+\mathit {\_F1} \left (\frac {a y -b x}{a}, \int \frac {a \left (-\arctan \left (\gamma z \right )+\frac {\pi }{2}\right )^{-k}}{c}d z -\left (\int _{}^{x}\left (-\arctan \left (\frac {\left (a y -\left (-\mathit {\_a} +x \right ) b \right ) \beta }{a}\right )+\frac {\pi }{2}\right )^{m} \left (-\arctan \left (\mathit {\_a} \lambda \right )+\frac {\pi }{2}\right )^{n}d\mathit {\_a} \right )\right )\]

____________________________________________________________________________________

6.7.23.5 [1721] Problem 5

problem number 1721

Added June 26, 2019.

Problem Chapter 7.7.4.5, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.

Solve for \(w(x,y,z)\)

\[ a w_x + b \arccot ^n(\lambda x) w_y + c \arccot ^k(\beta z) w_z = s \arccot ^m(\gamma x) \]

Mathematica

ClearAll["Global`*"]; 
pde =  a*D[w[x, y,z], x] + b*ArcTan[lambda*x]^n*D[w[x, y,z], y] + c*ArcTan[beta*z]^k*D[w[x,y,z],z]== s*ArcCot[gamma*x]^m; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y,z], {x, y,z}], 60*10]];
 

\[\left \{\left \{w(x,y,z)\to \int _1^z\frac {s \cot ^{-1}\left (\frac {\gamma \left (c x-a \int _1^z\tan ^{-1}(\beta K[2])^{-k}dK[2]+a \int _1^{K[3]}\tan ^{-1}(\beta K[2])^{-k}dK[2]\right )}{c}\right ){}^m \tan ^{-1}(\beta K[3])^{-k}}{c}dK[3]+c_1\left (y-\int _1^x\frac {b \tan ^{-1}(\lambda K[1])^n}{a}dK[1],\int _1^z\tan ^{-1}(\beta K[2])^{-k}dK[2]-\frac {c x}{a}\right )\right \}\right \}\]

Maple

restart; 
local gamma; 
pde :=  a*diff(w(x,y,z),x)+ b*arctan(lambda*x)^n*diff(w(x,y,z),y)+ c*arctan(beta*z)^k*diff(w(x,y,z),z)= s*arccot(gamma*x)^m; 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
 

\[w \left (x , y , z\right ) = \int _{}^{y}\frac {s \left (-\arctan \left (\gamma \RootOf \left (\mathit {\_b} -y +\int \frac {b \arctan \left (\lambda x \right )^{n}}{a}d x -\left (\int _{}^{\mathit {\_Z}}\frac {b \arctan \left (\mathit {\_a} \lambda \right )^{n}}{a}d\mathit {\_a} \right )\right )\right )+\frac {\pi }{2}\right )^{m} \arctan \left (\lambda \RootOf \left (\mathit {\_b} -y +\int \frac {b \arctan \left (\lambda x \right )^{n}}{a}d x -\left (\int _{}^{\mathit {\_Z}}\frac {b \arctan \left (\mathit {\_a} \lambda \right )^{n}}{a}d\mathit {\_a} \right )\right )\right )^{-n}}{b}d\mathit {\_b} +\mathit {\_F1} \left (-y +\int \frac {b \arctan \left (\lambda x \right )^{n}}{a}d x , \int \frac {b \arctan \left (\beta z \right )^{-k}}{c}d z -\left (\int _{}^{y}\arctan \left (\lambda \RootOf \left (\mathit {\_b} -y +\int \frac {b \arctan \left (\lambda x \right )^{n}}{a}d x -\left (\int _{}^{\mathit {\_Z}}\frac {b \arctan \left (\mathit {\_b} \lambda \right )^{n}}{a}d\mathit {\_b} \right )\right )\right )^{-n}d\mathit {\_b} \right )\right )\]

____________________________________________________________________________________