6.5.7 4.1

6.5.7.1 [1249] Problem 1

problem number 1249

Added April 3, 2019.

Problem Chapter 5.4.1.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 w + \sinh ^k(\lambda x) \sinh ^n(\beta y) \]

Mathematica ✓

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

\[\left \{\left \{w(x,y)\to e^{\frac {c x}{a}} \left (\int _1^x\frac {e^{-\frac {c K[1]}{a}} \sinh ^k(\lambda K[1]) \sinh ^n\left (\beta \left (y+\frac {b (K[1]-x)}{a}\right )\right )}{a}dK[1]+c_1\left (y-\frac {b x}{a}\right )\right )\right \}\right \}\]

Maple ✓

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

\[w \left ( x,y \right ) = \left ( \int ^{x}\!{\frac { \left ( \sinh \left ( {\it \_a}\,\lambda \right ) \right ) ^{k}}{a} \left ( \sinh \left ( {\frac {\beta \, \left ( ay-b \left ( x-{\it \_a} \right ) \right ) }{a}} \right ) \right ) ^{n}{{\rm e}^{-{\frac {c{\it \_a}}{a}}}}}{d{\it \_a}}+{\it \_F1} \left ( {\frac {ay-xb}{a}} \right ) \right ) {{\rm e}^{{\frac {cx}{a}}}}\]

____________________________________________________________________________________

6.5.7.2 [1250] Problem 2

problem number 1250

Added April 3, 2019.

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

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

\[ a w_x + b w_y = c \sinh ^k(\lambda x) w + s \sinh ^n(\beta x) \]

Mathematica ✓

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

\[\left \{\left \{w(x,y)\to \exp \left (\frac {c \sqrt {\cosh ^2(\lambda x)} \text {sech}(\lambda x) \sinh ^{k+1}(\lambda x) \, _2F_1\left (\frac {1}{2},\frac {k+1}{2};\frac {k+3}{2};-\sinh ^2(\lambda x)\right )}{a k \lambda +a \lambda }\right ) \left (\int _1^x\frac {\exp \left (-\frac {c \sqrt {\cosh ^2(\lambda K[1])} \, _2F_1\left (\frac {1}{2},\frac {k+1}{2};\frac {k+3}{2};-\sinh ^2(\lambda K[1])\right ) \text {sech}(\lambda K[1]) \sinh ^{k+1}(\lambda K[1])}{a \lambda +a k \lambda }\right ) s \sinh ^n(\beta K[1])}{a}dK[1]+c_1\left (y-\frac {b x}{a}\right )\right )\right \}\right \}\]

Maple ✓

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

\[w \left ( x,y \right ) = \left ( \int \!{\frac {s \left ( \sinh \left ( \beta \,x \right ) \right ) ^{n}}{a}{{\rm e}^{-{\frac {c\int \! \left ( \sinh \left ( x\lambda \right ) \right ) ^{k}\,{\rm d}x}{a}}}}}\,{\rm d}x+{\it \_F1} \left ( {\frac {ay-xb}{a}} \right ) \right ) {{\rm e}^{\int \!{\frac { \left ( \sinh \left ( x\lambda \right ) \right ) ^{k}c}{a}}\,{\rm d}x}}\]

____________________________________________________________________________________

6.5.7.3 [1251] Problem 3

problem number 1251

Added April 3, 2019.

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

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

\[ a w_x + b w_y = \left (c_1 \sinh ^{n_1}(\lambda _1 x)+ c_2 \sinh ^{n_2}(\lambda _2 y) \right ) w + s_1 \sinh ^{k_1}(\beta _1 x)+ s_2 \sinh ^{k_2}(\beta _2 y) \]

Mathematica ✓

ClearAll["Global`*"]; 
pde = a*D[w[x, y], x] + b*D[w[x, y], y] == (c1*Sinh[lambda1*x]^n1 + c2*Sinh[lambda2*y]^n2)*w[x,y] + s1*Sinh[beta1*x]^k1+ s2*Sinh[beta2*y]^k2; 
sol =  AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]];
 

\[\left \{\left \{w(x,y)\to \exp \left (\frac {\text {c1} \sqrt {\cosh ^2(\text {lambda1} x)} \text {sech}(\text {lambda1} x) \sinh ^{\text {n1}+1}(\text {lambda1} x) \, _2F_1\left (\frac {1}{2},\frac {\text {n1}+1}{2};\frac {\text {n1}+3}{2};-\sinh ^2(\text {lambda1} x)\right )}{a \text {lambda1} \text {n1}+a \text {lambda1}}+\frac {\text {c2} \sqrt {\cosh ^2(\text {lambda2} y)} \text {sech}(\text {lambda2} y) \sinh ^{\text {n2}+1}(\text {lambda2} y) \, _2F_1\left (\frac {1}{2},\frac {\text {n2}+1}{2};\frac {\text {n2}+3}{2};-\sinh ^2(\text {lambda2} y)\right )}{b \text {lambda2} \text {n2}+b \text {lambda2}}\right ) \left (\int _1^x\frac {\exp \left (-\frac {\text {c1} \sqrt {\cosh ^2(\text {lambda1} K[1])} \, _2F_1\left (\frac {1}{2},\frac {\text {n1}+1}{2};\frac {\text {n1}+3}{2};-\sinh ^2(\text {lambda1} K[1])\right ) \text {sech}(\text {lambda1} K[1]) \sinh ^{\text {n1}+1}(\text {lambda1} K[1])}{a \text {lambda1}+a \text {n1} \text {lambda1}}-\frac {\text {c2} \sqrt {\cosh ^2\left (\text {lambda2} \left (y+\frac {b (K[1]-x)}{a}\right )\right )} \, _2F_1\left (\frac {1}{2},\frac {\text {n2}+1}{2};\frac {\text {n2}+3}{2};-\sinh ^2\left (\text {lambda2} \left (y+\frac {b (K[1]-x)}{a}\right )\right )\right ) \text {sech}\left (\text {lambda2} \left (y+\frac {b (K[1]-x)}{a}\right )\right ) \sinh ^{\text {n2}+1}\left (\text {lambda2} \left (y+\frac {b (K[1]-x)}{a}\right )\right )}{b \text {lambda2}+b \text {n2} \text {lambda2}}\right ) \left (\text {s1} \sinh ^{\text {k1}}(\text {beta1} K[1])+\text {s2} \sinh ^{\text {k2}}\left (\text {beta2} \left (y+\frac {b (K[1]-x)}{a}\right )\right )\right )}{a}dK[1]+c_1\left (y-\frac {b x}{a}\right )\right )\right \}\right \}\]

Maple ✓

restart; 
pde :=  a*diff(w(x,y),x)+ b*diff(w(x,y),y) = (c1*sinh(lambda1*x)^n1 + c2*sinh(lambda2*y)^n2)*w(x,y) + s1*sinh(beta1*x)^k1+ s2*sinh(beta2*y)^k2; 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime'));
 

\[w \left ( x,y \right ) = \left ( \int ^{x}\!{\frac {1}{a} \left ( {\it s2}\, \left ( \sinh \left ( {\frac {\beta 2\, \left ( ay-b \left ( x-{\it \_b} \right ) \right ) }{a}} \right ) \right ) ^{{\it k2}}+{\it s1}\, \left ( \sinh \left ( \beta 1\,{\it \_b} \right ) \right ) ^{{\it k1}} \right ) {{\rm e}^{-{\frac {1}{a}\int \!{\it c1}\, \left ( \sinh \left ( \lambda 1\,{\it \_b} \right ) \right ) ^{{\it n1}}+{\it c2}\, \left ( \sinh \left ( {\frac {\lambda 2\, \left ( ay-b \left ( x-{\it \_b} \right ) \right ) }{a}} \right ) \right ) ^{{\it n2}}\,{\rm d}{\it \_b}}}}}{d{\it \_b}}+{\it \_F1} \left ( {\frac {ay-xb}{a}} \right ) \right ) {{\rm e}^{\int ^{x}\!{\frac {1}{a} \left ( {\it c1}\, \left ( \sinh \left ( \lambda 1\,{\it \_a} \right ) \right ) ^{{\it n1}}+{\it c2}\, \left ( \sinh \left ( {\frac {\lambda 2\, \left ( ay-b \left ( x-{\it \_a} \right ) \right ) }{a}} \right ) \right ) ^{{\it n2}} \right ) }{d{\it \_a}}}}\]

____________________________________________________________________________________

6.5.7.4 [1252] Problem 4

problem number 1252

Added April 3, 2019.

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

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

\[ a \sinh ^n(\lambda x) w_x + b \sinh ^m(\mu x) w_y = c \sinh ^k(\nu x) w + p \sinh ^s(\beta y) \]

Mathematica ✓

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

\[\left \{\left \{w(x,y)\to \exp \left (\frac {c 2^{n-1} e^{-\nu x} \left (e^{\lambda x}-e^{-\lambda x}\right )^{-n} \left (1-e^{2 \lambda x}\right )^n \left ((\lambda n+\nu ) \, _2F_1\left (n,\frac {\lambda n-\nu }{2 \lambda };\frac {1}{2} \left (n-\frac {\nu }{\lambda }+2\right );e^{2 \lambda x}\right )+e^{2 \nu x} (\nu -\lambda n) \, _2F_1\left (n,\frac {\lambda n+\nu }{2 \lambda };\frac {\lambda n+\nu }{2 \lambda }+1;e^{2 \lambda x}\right )\right )}{a (\nu -\lambda n) (\lambda n+\nu )}\right ) \left (\int _1^x\frac {\exp \left (-\frac {2^{n-1} c e^{-\nu K[2]} \left (-e^{-\lambda K[2]}+e^{\lambda K[2]}\right )^{-n} \left (1-e^{2 \lambda K[2]}\right )^n \left ((\lambda n+\nu ) \, _2F_1\left (n,\frac {\lambda n-\nu }{2 \lambda };\frac {1}{2} \left (n-\frac {\nu }{\lambda }+2\right );e^{2 \lambda K[2]}\right )+e^{2 \nu K[2]} (\nu -\lambda n) \, _2F_1\left (n,\frac {\lambda n+\nu }{2 \lambda };\frac {\lambda n+\nu }{2 \lambda }+1;e^{2 \lambda K[2]}\right )\right )}{a (\nu -\lambda n) (\lambda n+\nu )}\right ) p \sinh ^{-n}(\lambda K[2]) \sinh ^s\left (\beta \left (y-\int _1^x\frac {b \sinh ^{-n}(\lambda K[1]) \sinh ^m(\mu K[1])}{a}dK[1]+\int _1^{K[2]}\frac {b \sinh ^{-n}(\lambda K[1]) \sinh ^m(\mu K[1])}{a}dK[1]\right )\right )}{a}dK[2]+c_1\left (y-\int _1^x\frac {b \sinh ^{-n}(\lambda K[1]) \sinh ^m(\mu K[1])}{a}dK[1]\right )\right )\right \}\right \}\]

Maple ✓

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

\[w \left ( x,y \right ) ={\it \_F1} \left ( {\frac {-b\int \! \left ( \sinh \left ( \mu \,x \right ) \right ) ^{m} \left ( \sinh \left ( x\lambda \right ) \right ) ^{-n}\,{\rm d}x+ay}{a}} \right ) +\int ^{x}\!{\frac { \left ( \sinh \left ( {\it \_b}\,\lambda \right ) \right ) ^{-n}}{a} \left ( c\sinh \left ( \nu \,{\it \_b} \right ) w_{{x,y}}+ \left ( \sinh \left ( {\frac {\beta \, \left ( b\int \! \left ( \sinh \left ( \mu \,{\it \_b} \right ) \right ) ^{m} \left ( \sinh \left ( {\it \_b}\,\lambda \right ) \right ) ^{-n}\,{\rm d}{\it \_b}-b\int \! \left ( \sinh \left ( \mu \,x \right ) \right ) ^{m} \left ( \sinh \left ( x\lambda \right ) \right ) ^{-n}\,{\rm d}x+ay \right ) }{a}} \right ) \right ) ^{s}p \right ) }{d{\it \_b}}\]

____________________________________________________________________________________

6.5.7.5 [1253] Problem 5

problem number 1253

Added April 3, 2019.

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

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

\[ a \sinh ^n(\lambda x) w_x + b \sinh ^m(\mu x) w_y = c \sinh ^k(\nu y) w + p \sinh ^s(\beta x) \]

Mathematica ✓

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

\[\left \{\left \{w(x,y)\to \exp \left (\int _1^x\frac {c \sinh ^{-n}(\lambda K[2]) \sinh \left (\nu \left (y-\int _1^x\frac {b \sinh ^{-n}(\lambda K[1]) \sinh ^m(\mu K[1])}{a}dK[1]+\int _1^{K[2]}\frac {b \sinh ^{-n}(\lambda K[1]) \sinh ^m(\mu K[1])}{a}dK[1]\right )\right )}{a}dK[2]\right ) \left (\int _1^x\frac {\exp \left (-\int _1^{K[3]}\frac {c \sinh ^{-n}(\lambda K[2]) \sinh \left (\nu \left (y-\int _1^x\frac {b \sinh ^{-n}(\lambda K[1]) \sinh ^m(\mu K[1])}{a}dK[1]+\int _1^{K[2]}\frac {b \sinh ^{-n}(\lambda K[1]) \sinh ^m(\mu K[1])}{a}dK[1]\right )\right )}{a}dK[2]\right ) p \sinh ^s(\beta K[3]) \sinh ^{-n}(\lambda K[3])}{a}dK[3]+c_1\left (y-\int _1^x\frac {b \sinh ^{-n}(\lambda K[1]) \sinh ^m(\mu K[1])}{a}dK[1]\right )\right )\right \}\right \}\]

Maple ✓

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

\[w \left ( x,y \right ) ={\it \_F1} \left ( {\frac {-b\int \! \left ( \sinh \left ( \mu \,x \right ) \right ) ^{m} \left ( \sinh \left ( x\lambda \right ) \right ) ^{-n}\,{\rm d}x+ay}{a}} \right ) +\int ^{x}\!{\frac { \left ( \sinh \left ( {\it \_b}\,\lambda \right ) \right ) ^{-n}}{a} \left ( \sinh \left ( {\frac {\nu \, \left ( b\int \! \left ( \sinh \left ( \mu \,{\it \_b} \right ) \right ) ^{m} \left ( \sinh \left ( {\it \_b}\,\lambda \right ) \right ) ^{-n}\,{\rm d}{\it \_b}-b\int \! \left ( \sinh \left ( \mu \,x \right ) \right ) ^{m} \left ( \sinh \left ( x\lambda \right ) \right ) ^{-n}\,{\rm d}x+ay \right ) }{a}} \right ) w_{{x,y}}c+p \left ( \sinh \left ( \beta \,{\it \_b} \right ) \right ) ^{s} \right ) }{d{\it \_b}}\]

____________________________________________________________________________________