Added June 26, 2019.
Problem Chapter 7.7.1.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 \arcsin ^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*ArcSin[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 c_1(y-a x,z-b x)+\frac {\left (\sin ^{-1}(\lambda x)^2\right )^{-k} \left (-i c \left (i \sin ^{-1}(\lambda x)\right )^k \sin ^{-1}(\lambda x)^k \operatorname {Gamma}\left (k+1,-i \sin ^{-1}(\lambda x)\right )+i c \left (-i \sin ^{-1}(\lambda x)\right )^k \sin ^{-1}(\lambda x)^k \operatorname {Gamma}\left (k+1,i \sin ^{-1}(\lambda x)\right )+2 \lambda s x \left (\sin ^{-1}(\lambda x)^2\right )^k\right )}{2 \lambda }\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*arcsin(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 \arcsin \left (\lambda x \right )^{k}d x +\textit {\_F1} \left (-a x +y , -b x +z \right )\]
____________________________________________________________________________________
Added June 26, 2019.
Problem Chapter 7.7.1.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 \arcsin (\lambda _1 x)+b_2 \arcsin (\lambda _2 y)+b_3 \arcsin (\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*ArcSin[lambda1*x]+b2*ArcSin[lambda2*y]+b3*ArcSin[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 {\text {b1} \sqrt {1-\text {lambda1}^2 x^2}}{\text {a1} \text {lambda1}}+\frac {\text {b1} x \sin ^{-1}(\text {lambda1} x)}{\text {a1}}+\frac {\text {b2} \sqrt {1-\text {lambda2}^2 y^2}}{\text {a2} \text {lambda2}}+\frac {\text {b2} y \sin ^{-1}(\text {lambda2} y)}{\text {a2}}+\frac {\text {b3} \sqrt {1-\text {lambda3}^2 z^2}}{\text {a3} \text {lambda3}}+\frac {\text {b3} z \sin ^{-1}(\text {lambda3} z)}{\text {a3}}\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*arcsin(lambda1*x)+b2*arcsin(lambda2*y)+b3*arcsin(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 {\mathit {a1} \mathit {a2} \mathit {a3} \lambda 1 \lambda 2 \lambda 3 \textit {\_F1} \left (\frac {\mathit {a1} y -\mathit {a2} x}{\mathit {a1}}, \frac {\mathit {a1} z -\mathit {a3} x}{\mathit {a1}}\right )+\sqrt {-\lambda 1^{2} x^{2}+1}\, \mathit {a2} \mathit {a3} \mathit {b1} \lambda 2 \lambda 3 +\left (\sqrt {-\lambda 2^{2} y^{2}+1}\, \mathit {a1} \mathit {a3} \mathit {b2} \lambda 3 +\left (\mathit {a2} \mathit {a3} \mathit {b1} \lambda 3 x \arcsin \left (\lambda 1 x \right )+\left (\mathit {a3} \mathit {b2} \lambda 3 y \arcsin \left (\lambda 2 y \right )+\left (\lambda 3 z \arcsin \left (\lambda 3 z \right )+\sqrt {-\lambda 3^{2} z^{2}+1}\right ) \mathit {a2} \mathit {b3} \right ) \mathit {a1} \right ) \lambda 2 \right ) \lambda 1}{\mathit {a1} \mathit {a2} \mathit {a3} \lambda 1 \lambda 2 \lambda 3}\]
____________________________________________________________________________________
Added June 26, 2019.
Problem Chapter 7.7.1.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 \arcsin ^n(\lambda x) \arcsin ^k(\beta z) w_z = s \arcsin ^m(\gamma x) \]
Mathematica ✗
ClearAll["Global`*"]; pde = a*D[w[x, y,z], x] + b*D[w[x, y,z], y] + c*ArcSin[lambda*x]^n*ArcSin[beta*z]^k*D[w[x,y,z],z]== s*ArcSin[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*arcsin(lambda*x)^n*arcsin(beta*z)^k*diff(w(x,y,z),z)= s*arcsin(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 \arcsin \left (\gamma x \right )^{m}}{a}d x +\textit {\_F1} \left (\frac {a y -b x}{a}, -\frac {2 \left (\frac {\left (k -1\right ) \left (\arcsin \left (\lambda x \right )^{n}-\frac {\LommelS 1 \left (n +\frac {3}{2}, \frac {1}{2}, \arcsin \left (\lambda x \right )\right )}{\sqrt {\arcsin \left (\lambda x \right )}}\right ) \left (-\lambda ^{2} x^{2}+1\right ) \beta c \lambda x 2^{n} 2^{-n}}{2}+\frac {\left (\lambda x -1\right ) \left (\lambda x +1\right ) \left (k -1\right ) \left (\arcsin \left (\lambda x \right )^{n}-\frac {\LommelS 1 \left (n +\frac {3}{2}, \frac {1}{2}, \arcsin \left (\lambda x \right )\right )}{\sqrt {\arcsin \left (\lambda x \right )}}\right ) \sqrt {-\lambda ^{2} x^{2}+1}\, \beta c 2^{n} 2^{-n} \arcsin \left (\lambda x \right )}{2}+\left (\lambda x -1\right ) \left (-\frac {\left (n +1\right ) \left (-\arcsin \left (\beta z \right )^{-k} \arcsin \left (\beta z \right )^{\frac {3}{2}}+\LommelS 1 \left (-k +\frac {3}{2}, \frac {1}{2}, \arcsin \left (\beta z \right )\right ) \arcsin \left (\beta z \right )\right ) \sqrt {-\beta ^{2} z^{2}+1}\, a 2^{k} 2^{-k}}{2 \sqrt {\arcsin \left (\beta z \right )}}+\left (-\frac {\left (n +1\right ) a k z 2^{k} 2^{-k} \LommelS 1 \left (-k +\frac {1}{2}, \frac {3}{2}, \arcsin \left (\beta z \right )\right ) \sqrt {\arcsin \left (\beta z \right )}}{2}+\frac {\left (k -1\right ) c n x 2^{n} 2^{-n} \LommelS 1 \left (n +\frac {1}{2}, \frac {3}{2}, \arcsin \left (\lambda x \right )\right ) \sqrt {\arcsin \left (\lambda x \right )}}{2}+\left (k -1\right ) c x 2^{n} 2^{-n -1} \arcsin \left (\lambda x \right )^{n}+\frac {\left (n +1\right ) a z 2^{k} 2^{-k} \LommelS 1 \left (-k +\frac {3}{2}, \frac {1}{2}, \arcsin \left (\beta z \right )\right )}{2 \sqrt {\arcsin \left (\beta z \right )}}+\left (n +1\right ) \left (-\frac {2^{k}}{2}+2^{k -1}\right ) a z 2^{-k} \arcsin \left (\beta z \right )^{-k}\right ) \beta \right ) \left (\lambda x +1\right ) \lambda \right )}{\left (n +1\right ) \left (\lambda ^{2} x^{2}-1\right ) \left (k -1\right ) \beta c \lambda }\right )\]
____________________________________________________________________________________
Added June 26, 2019.
Problem Chapter 7.7.1.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 \arcsin ^n(\lambda x) \arcsin ^m(\beta y) \arcsin ^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*ArcSin[lambda*x]^n*ArcSin[beta*y]^m*ArcSin[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*arcsin(lambda*x)^n*arcsin(beta*y)^m*arcsin(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}+\textit {\_F1} \left (\frac {a y -b x}{a}, -\left (\int _{}^{x}\arcsin \left (\textit {\_a} \lambda \right )^{n} \arcsin \left (\frac {\left (a y -\left (-\textit {\_a} +x \right ) b \right ) \beta }{a}\right )^{m}d \textit {\_a} \right )+\frac {\left (\gamma k z 2^{k} \LommelS 1 \left (-k +\frac {1}{2}, \frac {3}{2}, \arcsin \left (\gamma z \right )\right ) \arcsin \left (\gamma z \right )-\gamma z 2^{k} \LommelS 1 \left (-k +\frac {3}{2}, \frac {1}{2}, \arcsin \left (\gamma z \right )\right )+\sqrt {-\gamma ^{2} z^{2}+1}\, 2^{k} \LommelS 1 \left (-k +\frac {3}{2}, \frac {1}{2}, \arcsin \left (\gamma z \right )\right ) \arcsin \left (\gamma z \right )-\sqrt {-\gamma ^{2} z^{2}+1}\, 2^{k} \arcsin \left (\gamma z \right )^{-k +\frac {3}{2}}\right ) a 2^{-k}}{\left (k -1\right ) c \gamma \sqrt {\arcsin \left (\gamma z \right )}}\right )\]
____________________________________________________________________________________
Added June 26, 2019.
Problem Chapter 7.7.1.5, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y,z)\)
\[ a w_x + b \arcsin ^n(\lambda x) w_y + c \arcsin ^k(\beta z) w_z = s \arcsin ^m(\gamma x) \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*D[w[x, y,z], x] + b*ArcSin[lambda*x]^n*D[w[x, y,z], y] + c*ArcSin[beta*z]^k*D[w[x,y,z],z]== s*ArcSin[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 \sin ^{-1}\left (\frac {\gamma \left (i a \left (-i \sin ^{-1}(\beta z)\right )^k \operatorname {Gamma}\left (1-k,-i \sin ^{-1}(\beta z)\right ) \sin ^{-1}(\beta z)^{-k}-i a \left (i \sin ^{-1}(\beta z)\right )^k \operatorname {Gamma}\left (1-k,i \sin ^{-1}(\beta z)\right ) \sin ^{-1}(\beta z)^{-k}+\sin ^{-1}(\beta K[1])^{-k} \left (-i a \operatorname {Gamma}\left (1-k,-i \sin ^{-1}(\beta K[1])\right ) \left (-i \sin ^{-1}(\beta K[1])\right )^k+2 \beta c x \sin ^{-1}(\beta K[1])^k+i a \left (i \sin ^{-1}(\beta K[1])\right )^k \operatorname {Gamma}\left (1-k,i \sin ^{-1}(\beta K[1])\right )\right )\right )}{2 \beta c}\right )^m \sin ^{-1}(\beta K[1])^{-k}}{c}dK[1]+c_1\left (-\frac {c x}{a}-\frac {i \sin ^{-1}(\beta z)^{-k} \left (\left (-i \sin ^{-1}(\beta z)\right )^k \operatorname {Gamma}\left (1-k,-i \sin ^{-1}(\beta z)\right )-\left (i \sin ^{-1}(\beta z)\right )^k \operatorname {Gamma}\left (1-k,i \sin ^{-1}(\beta z)\right )\right )}{2 \beta },\frac {\left (\sin ^{-1}(\lambda x)^2\right )^{-n} \left (i b \left (i \sin ^{-1}(\lambda x)\right )^n \sin ^{-1}(\lambda x)^n \operatorname {Gamma}\left (n+1,-i \sin ^{-1}(\lambda x)\right )-i b \left (-i \sin ^{-1}(\lambda x)\right )^n \sin ^{-1}(\lambda x)^n \operatorname {Gamma}\left (n+1,i \sin ^{-1}(\lambda x)\right )+2 a \lambda y \left (\sin ^{-1}(\lambda x)^2\right )^n\right )}{2 a \lambda }\right )\right \}\right \}\] Generates Solve::incnst: Inconsistent or redundant transcendental equation
Maple ✗
restart; local gamma; pde := a*diff(w(x,y,z),x)+ b*arcsin(lambda*x)^n*diff(w(x,y,z),y)+ c*arcsin(beta*z)^k*diff(w(x,y,z),z)= s*arcsin(gamma*x)^m; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y,z))),output='realtime'));
time expired
____________________________________________________________________________________
Added June 26, 2019.
Problem Chapter 7.7.1.6, from Handbook of first order partial differential equations by Polyanin, Zaitsev, Moussiaux.
Solve for \(w(x,y,z)\)
\[ a w_x + b \arcsin ^n(\lambda x) w_y + c \arcsin ^k(\beta z) w_z = s \]
Mathematica ✓
ClearAll["Global`*"]; pde = a*D[w[x, y,z], x] + b*ArcSin[lambda*x]^n*D[w[x, y,z], y] + c*ArcSin[beta*z]^k*D[w[x,y,z],z]== s; 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 {c x}{a}-\frac {i \sin ^{-1}(\beta z)^{-k} \left (\left (-i \sin ^{-1}(\beta z)\right )^k \operatorname {Gamma}\left (1-k,-i \sin ^{-1}(\beta z)\right )-\left (i \sin ^{-1}(\beta z)\right )^k \operatorname {Gamma}\left (1-k,i \sin ^{-1}(\beta z)\right )\right )}{2 \beta },\frac {\left (\sin ^{-1}(\lambda x)^2\right )^{-n} \left (i b \left (i \sin ^{-1}(\lambda x)\right )^n \sin ^{-1}(\lambda x)^n \operatorname {Gamma}\left (n+1,-i \sin ^{-1}(\lambda x)\right )-i b \left (-i \sin ^{-1}(\lambda x)\right )^n \sin ^{-1}(\lambda x)^n \operatorname {Gamma}\left (n+1,i \sin ^{-1}(\lambda x)\right )+2 a \lambda y \left (\sin ^{-1}(\lambda x)^2\right )^n\right )}{2 a \lambda }\right )-\frac {i s \sin ^{-1}(\beta z)^{-k} \left (\left (-i \sin ^{-1}(\beta z)\right )^k \operatorname {Gamma}\left (1-k,-i \sin ^{-1}(\beta z)\right )-\left (i \sin ^{-1}(\beta z)\right )^k \operatorname {Gamma}\left (1-k,i \sin ^{-1}(\beta z)\right )\right )}{2 \beta c}\right \}\right \}\] Generates Solve::incnst: Inconsistent or redundant transcendental equation
Maple ✗
restart; local gamma; pde := a*diff(w(x,y,z),x)+ b*arcsin(lambda*x)^n*diff(w(x,y,z),y)+ c*arcsin(beta*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'));
time expired
____________________________________________________________________________________