140 HFOPDE, chapter 5.6.3

____________________________________________________________________________________

140.1 Problem 1

problem number 1111

Added April 11, 2019.

Problem Chapter 5.6.3.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 + k \tan (\lambda x+\mu y) \]

Mathematica ✓

ClearAll[w, x, y, n, a, b, m, c, k, alpha, beta, gamma, A, C0, s]; 
 ClearAll[lambda, B, mu, d, g, B, v, f, h, q, p, delta, t]; 
 ClearAll[g1, g0, h2, h1, h0, f1, f2,sigma,lambda1,lambda2,n1,n2]; 
 ClearAll[a1, a0, b2, b1, b0, c2, c1, c0, k0, k1, k2, s1, s0, k22, k11, k12, s11, s22, s12, nu]; 
 pde = a*D[w[x, y], x] + b*D[w[x, y], y] == c*w[x,y]+k*Tan[lambda*x+mu*y]; 
 sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]]; 
 sol = Simplify[sol];
 

\[ \left \{\left \{w(x,y)\to e^{\frac {c x}{a}} c_1\left (y-\frac {b x}{a}\right )-\frac {i k \left (-2 \text {Hypergeometric2F1}\left (1,-\frac {i c}{2 a \lambda +2 b \mu },\frac {2 a \lambda +2 b \mu -i c}{2 (a \lambda +b \mu )},-e^{-2 i (\lambda x+\mu y)}\right )+2 e^{\frac {2 i \mu (a y-b x)}{a}} \text {Hypergeometric2F1}\left (1,\frac {i c}{2 a \lambda +2 b \mu },\frac {2 a \lambda +2 b \mu +i c}{2 a \lambda +2 b \mu },-e^{2 i (\lambda x+\mu y)}\right )-e^{\frac {2 i \mu (a y-b x)}{a}}+1\right )}{c \left (1+e^{\frac {2 i \mu (a y-b x)}{a}}\right )}\right \}\right \} \]

Maple ✓

 
unassign('w,x,y,a,b,n,m,c,k,alpha,beta,g,A,f,C,lambda,B,mu,d,s,t'); 
unassign('v,q,p,l,g1,g2,g0,h0,h1,h2,f2,f3,c0,c1,c2,a1,a0,b0,b1,b2'); 
unassign('k0,k1,k2,s0,s1,k22,k12,k11,s22,s12,s11,sigma,lambda1,lambda2,n1,n2,nu'); 
pde :=  a*diff(w(x,y),x)+ b*diff(w(x,y),y) =c*w(x,y)+k*tan(lambda*x+mu*y); 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime')); 
if(not evalb(sol=())) then sol:=simplify(sol,size); fi;
 

\[ w \left ( x,y \right ) ={{\rm e}^{{\frac {cx}{a}}}} \left ( \int ^{x}\!{\frac {k}{a}\tan \left ( {\frac { \left ( \lambda \,{\it \_a}+\mu \,y \right ) a-b\mu \, \left ( x-{\it \_a} \right ) }{a}} \right ) {{\rm e}^{-{\frac {{\it \_a}\,c}{a}}}}}{d{\it \_a}}+{\it \_F1} \left ( {\frac {ya-bx}{a}} \right ) \right ) \]

____________________________________________________________________________________

140.2 Problem 2

problem number 1112

Added April 11, 2019.

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

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

\[ a w_x + b w_y = w + c_1 \tan ^k(\lambda x) + c_2 \tan ^n(\beta y) \]

Mathematica ✗

ClearAll[w, x, y, n, a, b, m, c, k, alpha, beta, gamma, A, C0, s]; 
 ClearAll[lambda, B, mu, d, g, B, v, f, h, q, p, delta, t]; 
 ClearAll[g1, g0, h2, h1, h0, f1, f2,sigma,lambda1,lambda2,n1,n2]; 
 ClearAll[a1, a0, b2, b1, b0, c2, c1, c0, k0, k1, k2, s1, s0, k22, k11, k12, s11, s22, s12, nu]; 
 pde = a*D[w[x, y], x] + b*D[w[x, y], y] == w[x,y]+ c1*Tan[lambda*x]^k + c2*Tan[beta*y]^n; 
 sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]]; 
 sol = Simplify[sol];
 

\[ \text {\$Aborted} \]

Maple ✓

 
unassign('w,x,y,a,b,n,m,c,k,alpha,beta,g,A,f,C,lambda,B,mu,d,s,t'); 
unassign('v,q,p,l,g1,g2,g0,h0,h1,h2,f2,f3,c0,c1,c2,a1,a0,b0,b1,b2'); 
unassign('k0,k1,k2,s0,s1,k22,k12,k11,s22,s12,s11,sigma,lambda1,lambda2,n1,n2,nu'); 
pde :=  a*diff(w(x,y),x)+ b*diff(w(x,y),y) = w(x,y)+ c1*tan(lambda*x)^k + c2*tan(beta*y)^n; 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime')); 
if(not evalb(sol=())) then sol:=simplify(sol,size); fi;
 

\[ w \left ( x,y \right ) ={{\rm e}^{{\frac {x}{a}}}} \left ( {\it \_F1} \left ( {\frac {ya-bx}{a}} \right ) +\int ^{x}\!{\frac {1}{a}{{\rm e}^{-{\frac {{\it \_a}}{a}}}} \left ( {\it c2}\, \left ( \tan \left ( {\frac {\beta \, \left ( ya-b \left ( x-{\it \_a} \right ) \right ) }{a}} \right ) \right ) ^{n}+{\it c1}\, \left ( \tan \left ( \lambda \,{\it \_a} \right ) \right ) ^{k} \right ) }{d{\it \_a}} \right ) \]

____________________________________________________________________________________

140.3 Problem 3

problem number 1113

Added April 11, 2019.

Problem Chapter 5.6.3.3, 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 + \tan ^k(\lambda x) \tan ^n(\beta y) \]

Mathematica ✓

ClearAll[w, x, y, n, a, b, m, c, k, alpha, beta, gamma, A, C0, s]; 
 ClearAll[lambda, B, mu, d, g, B, v, f, h, q, p, delta, t]; 
 ClearAll[g1, g0, h2, h1, h0, f1, f2,sigma,lambda1,lambda2,n1,n2]; 
 ClearAll[a1, a0, b2, b1, b0, c2, c1, c0, k0, k1, k2, s1, s0, k22, k11, k12, s11, s22, s12, nu]; 
 pde = a*D[w[x, y], x] + b*D[w[x, y], y] == c*w[x,y]+ Tan[lambda*x]^k * Tan[beta*y]^n; 
 sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]]; 
 sol = Simplify[sol];
 

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

Maple ✓

 
unassign('w,x,y,a,b,n,m,c,k,alpha,beta,g,A,f,C,lambda,B,mu,d,s,t'); 
unassign('v,q,p,l,g1,g2,g0,h0,h1,h2,f2,f3,c0,c1,c2,a1,a0,b0,b1,b2'); 
unassign('k0,k1,k2,s0,s1,k22,k12,k11,s22,s12,s11,sigma,lambda1,lambda2,n1,n2,nu'); 
pde :=  a*diff(w(x,y),x)+ b*diff(w(x,y),y) = c*w(x,y)+ tan(lambda*x)^k *tan(beta*y)^n; 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime')); 
if(not evalb(sol=())) then sol:=simplify(sol,size); fi;
 

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

____________________________________________________________________________________

140.4 Problem 4

problem number 1114

Added April 11, 2019.

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

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

\[ a w_x + b \tan (\mu y) w_y = c \tan (\lambda x) w + k \tan (\nu x) \]

Mathematica ✓

ClearAll[w, x, y, n, a, b, m, c, k, alpha, beta, gamma, A, C0, s]; 
 ClearAll[lambda, B, mu, d, g, B, v, f, h, q, p, delta, t]; 
 ClearAll[g1, g0, h2, h1, h0, f1, f2,sigma,lambda1,lambda2,n1,n2]; 
 ClearAll[a1, a0, b2, b1, b0, c2, c1, c0, k0, k1, k2, s1, s0, k22, k11, k12, s11, s22, s12, nu]; 
 pde = a*D[w[x, y], x] + b*Tan[mu*y]*D[w[x, y], y] == c*Tan[lambda*x]*w[x,y]+k*Tan[nu*x]; 
 sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]]; 
 sol = Simplify[sol];
 

\[ \left \{\left \{w(x,y)\to \cos ^{-\frac {c}{a \lambda }}(\lambda x) \left (\int _1^x \frac {k \tan (\nu K[1]) \cos ^{\frac {c}{a \lambda }}(\lambda K[1])}{a} \, dK[1]+c_1\left (\frac {\log (\sin (\mu y))}{\mu }-\frac {b x}{a}\right )\right )\right \}\right \} \]

Maple ✓

 
unassign('w,x,y,a,b,n,m,c,k,alpha,beta,g,A,f,C,lambda,B,mu,d,s,t'); 
unassign('v,q,p,l,g1,g2,g0,h0,h1,h2,f2,f3,c0,c1,c2,a1,a0,b0,b1,b2'); 
unassign('k0,k1,k2,s0,s1,k22,k12,k11,s22,s12,s11,sigma,lambda1,lambda2,n1,n2,nu'); 
pde :=  a*diff(w(x,y),x)+ b*tan(mu*y)*diff(w(x,y),y) = c*tan(lambda*x)*w(x,y)+ k*tan(nu*x); 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime')); 
if(not evalb(sol=())) then sol:=simplify(sol,size); fi;
 

\[ w \left ( x,y \right ) = \left ( 1+ \left ( \tan \left ( \lambda \,x \right ) \right ) ^{2} \right ) ^{1/2\,{\frac {c}{a\lambda }}} \left ( \int \!{\frac {k\sin \left ( \nu \,x \right ) }{a\cos \left ( \nu \,x \right ) } \left ( 2\, \left ( \cos \left ( 2\,\lambda \,x \right ) +1 \right ) ^{-1} \right ) ^{-1/2\,{\frac {c}{a\lambda }}}}\,{\rm d}x+{\it \_F1} \left ( {\frac {1}{b\mu } \left ( -b\mu \,x+\ln \left ( {\frac {\tan \left ( \mu \,y \right ) }{\sqrt {1+ \left ( \tan \left ( \mu \,y \right ) \right ) ^{2}}}} \right ) a \right ) } \right ) \right ) \]

____________________________________________________________________________________

140.5 Problem 5

problem number 1115

Added April 11, 2019.

Problem Chapter 5.6.3.5, 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 w + k \tan (\lambda x+\nu y) \]

Mathematica ✓

ClearAll[w, x, y, n, a, b, m, c, k, alpha, beta, gamma, A, C0, s]; 
 ClearAll[lambda, B, mu, d, g, B, v, f, h, q, p, delta, t]; 
 ClearAll[g1, g0, h2, h1, h0, f1, f2,sigma,lambda1,lambda2,n1,n2]; 
 ClearAll[a1, a0, b2, b1, b0, c2, c1, c0, k0, k1, k2, s1, s0, k22, k11, k12, s11, s22, s12, nu]; 
 pde = a*x*D[w[x, y], x] + b*y*D[w[x, y], y] == c*w[x,y]+k*Tan[lambda*x+nu*y]; 
 sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]]; 
 sol = Simplify[sol];
 

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

Maple ✓

 
unassign('w,x,y,a,b,n,m,c,k,alpha,beta,g,A,f,C,lambda,B,mu,d,s,t'); 
unassign('v,q,p,l,g1,g2,g0,h0,h1,h2,f2,f3,c0,c1,c2,a1,a0,b0,b1,b2'); 
unassign('k0,k1,k2,s0,s1,k22,k12,k11,s22,s12,s11,sigma,lambda1,lambda2,n1,n2,nu'); 
pde :=  a*x*diff(w(x,y),x)+ b*y*diff(w(x,y),y) =c*w(x,y)+k*tan(lambda*x+nu*y); 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime')); 
if(not evalb(sol=())) then sol:=simplify(sol,size); fi;
 

\[ w \left ( x,y \right ) ={x}^{{\frac {c}{a}}} \left ( {\it \_F1} \left ( y{x}^{-{\frac {b}{a}}} \right ) +\int ^{x}\!{\frac {k}{a}\tan \left ( \lambda \,{\it \_a}+\nu \,y{x}^{-{\frac {b}{a}}}{{\it \_a}}^{{\frac {b}{a}}} \right ) {{\it \_a}}^{{\frac {-a-c}{a}}}}{d{\it \_a}} \right ) \]

____________________________________________________________________________________

140.6 Problem 6

problem number 1116

Added April 11, 2019.

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

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

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

Mathematica ✗

ClearAll[w, x, y, n, a, b, m, c, k, alpha, beta, gamma, A, C0, s]; 
 ClearAll[lambda, B, mu, d, g, B, v, f, h, q, p, delta, t]; 
 ClearAll[g1, g0, h2, h1, h0, f1, f2,sigma,lambda1,lambda2,n1,n2]; 
 ClearAll[a1, a0, b2, b1, b0, c2, c1, c0, k0, k1, k2, s1, s0, k22, k11, k12, s11, s22, s12, nu]; 
 pde = a*Tan[lambda*x]^n*D[w[x, y], x] + b*Tan[mu*x]^m*D[w[x, y], y] == c*Tan[nu*x]^k*w[x,y]+p*Tan[beta*y]^s; 
 sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]]; 
 sol = Simplify[sol];
 

\[ \text {\$Aborted} \]

Maple ✓

 
unassign('w,x,y,a,b,n,m,c,k,alpha,beta,g,A,f,C,lambda,B,mu,d,s,t'); 
unassign('v,q,p,l,g1,g2,g0,h0,h1,h2,f2,f3,c0,c1,c2,a1,a0,b0,b1,b2'); 
unassign('k0,k1,k2,s0,s1,k22,k12,k11,s22,s12,s11,sigma,lambda1,lambda2,n1,n2,nu'); 
pde :=  a*tan(lambda*x)^n*diff(w(x,y),x)+ b*tan(mu*x)^m*diff(w(x,y),y) =c*tan(nu*x)^k*w(x,y)+p*tan(beta*y)^s; 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime')); 
if(not evalb(sol=())) then sol:=simplify(sol,size); fi;
 

\[ w \left ( x,y \right ) ={{\rm e}^{\int \!{\frac { \left ( \tan \left ( \nu \,x \right ) \right ) ^{k}c \left ( \tan \left ( \lambda \,x \right ) \right ) ^{-n}}{a}}\,{\rm d}x}} \left ( \int ^{x}\!{\frac {p}{a} \left ( {1\sin \left ( {\frac {\beta }{a} \left ( b\int \! \left ( {\frac {\sin \left ( {\it \_f}\,\mu \right ) }{\cos \left ( {\it \_f}\,\mu \right ) }} \right ) ^{m} \left ( {\frac {\sin \left ( \lambda \,{\it \_f} \right ) }{\cos \left ( \lambda \,{\it \_f} \right ) }} \right ) ^{-n}\,{\rm d}{\it \_f}+ya-b\int \! \left ( {\frac {\sin \left ( \mu \,x \right ) }{\cos \left ( \mu \,x \right ) }} \right ) ^{m} \left ( {\frac {\sin \left ( \lambda \,x \right ) }{\cos \left ( \lambda \,x \right ) }} \right ) ^{-n}\,{\rm d}x \right ) } \right ) \left ( \cos \left ( {\frac {\beta }{a} \left ( b\int \! \left ( {\frac {\sin \left ( {\it \_f}\,\mu \right ) }{\cos \left ( {\it \_f}\,\mu \right ) }} \right ) ^{m} \left ( {\frac {\sin \left ( \lambda \,{\it \_f} \right ) }{\cos \left ( \lambda \,{\it \_f} \right ) }} \right ) ^{-n}\,{\rm d}{\it \_f}+ya-b\int \! \left ( {\frac {\sin \left ( \mu \,x \right ) }{\cos \left ( \mu \,x \right ) }} \right ) ^{m} \left ( {\frac {\sin \left ( \lambda \,x \right ) }{\cos \left ( \lambda \,x \right ) }} \right ) ^{-n}\,{\rm d}x \right ) } \right ) \right ) ^{-1}} \right ) ^{s} \left ( {\frac {\sin \left ( \lambda \,{\it \_f} \right ) }{\cos \left ( \lambda \,{\it \_f} \right ) }} \right ) ^{-n}{{\rm e}^{-{\frac {c}{a}\int \! \left ( {\frac {\sin \left ( \nu \,{\it \_f} \right ) }{\cos \left ( \nu \,{\it \_f} \right ) }} \right ) ^{k} \left ( {\frac {\sin \left ( \lambda \,{\it \_f} \right ) }{\cos \left ( \lambda \,{\it \_f} \right ) }} \right ) ^{-n}\,{\rm d}{\it \_f}}}}}{d{\it \_f}}+{\it \_F1} \left ( {\frac {1}{a} \left ( ya-b\int \! \left ( {\frac {\sin \left ( \mu \,x \right ) }{\cos \left ( \mu \,x \right ) }} \right ) ^{m} \left ( {\frac {\sin \left ( \lambda \,x \right ) }{\cos \left ( \lambda \,x \right ) }} \right ) ^{-n}\,{\rm d}x \right ) } \right ) \right ) \]

____________________________________________________________________________________

140.7 Problem 7

problem number 1117

Added April 11, 2019.

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

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

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

Mathematica ✗

ClearAll[w, x, y, n, a, b, m, c, k, alpha, beta, gamma, A, C0, s]; 
 ClearAll[lambda, B, mu, d, g, B, v, f, h, q, p, delta, t]; 
 ClearAll[g1, g0, h2, h1, h0, f1, f2,sigma,lambda1,lambda2,n1,n2]; 
 ClearAll[a1, a0, b2, b1, b0, c2, c1, c0, k0, k1, k2, s1, s0, k22, k11, k12, s11, s22, s12, nu]; 
 pde = a*Tan[lambda*x]^n*D[w[x, y], x] + b*Tan[mu*x]^m*D[w[x, y], y] == c*Tan[nu*y]^k*w[x,y]+p*Tan[beta*x]^s; 
 sol = AbsoluteTiming[TimeConstrained[DSolve[pde, w[x, y], {x, y}], 60*10]]; 
 sol = Simplify[sol];
 

\[ \text {\$Aborted} \]

Maple ✓

 
unassign('w,x,y,a,b,n,m,c,k,alpha,beta,g,A,f,C,lambda,B,mu,d,s,t'); 
unassign('v,q,p,l,g1,g2,g0,h0,h1,h2,f2,f3,c0,c1,c2,a1,a0,b0,b1,b2'); 
unassign('k0,k1,k2,s0,s1,k22,k12,k11,s22,s12,s11,sigma,lambda1,lambda2,n1,n2,nu'); 
pde :=  a*tan(lambda*x)^n*diff(w(x,y),x)+ b*tan(mu*x)^m*diff(w(x,y),y) =c*tan(nu*y)^k*w(x,y)+p*tan(beta*x)^s; 
cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime')); 
if(not evalb(sol=())) then sol:=simplify(sol,size); fi;
 

\[ w \left ( x,y \right ) ={{\rm e}^{\int ^{x}\!{\frac {c \left ( \tan \left ( {\it \_b}\,\lambda \right ) \right ) ^{-n}}{a} \left ( \tan \left ( {\frac {\nu }{a} \left ( -b\int \! \left ( {\frac {\sin \left ( \mu \,x \right ) }{\cos \left ( \mu \,x \right ) }} \right ) ^{m} \left ( {\frac {\sin \left ( \lambda \,x \right ) }{\cos \left ( \lambda \,x \right ) }} \right ) ^{-n}\,{\rm d}x+a \left ( \int \!{\frac {b \left ( \tan \left ( {\it \_b}\,\mu \right ) \right ) ^{m} \left ( \tan \left ( {\it \_b}\,\lambda \right ) \right ) ^{-n}}{a}}\,{\rm d}{\it \_b}+y \right ) \right ) } \right ) \right ) ^{k}}{d{\it \_b}}}} \left ( \int ^{x}\!{\frac {p}{a} \left ( {\frac {\sin \left ( \beta \,{\it \_f} \right ) }{\cos \left ( \beta \,{\it \_f} \right ) }} \right ) ^{s} \left ( {\frac {\sin \left ( \lambda \,{\it \_f} \right ) }{\cos \left ( \lambda \,{\it \_f} \right ) }} \right ) ^{-n}{{\rm e}^{-{\frac {c}{a}\int \! \left ( {1\sin \left ( {\frac {\nu }{a} \left ( b\int \! \left ( {\frac {\sin \left ( {\it \_f}\,\mu \right ) }{\cos \left ( {\it \_f}\,\mu \right ) }} \right ) ^{m} \left ( {\frac {\sin \left ( \lambda \,{\it \_f} \right ) }{\cos \left ( \lambda \,{\it \_f} \right ) }} \right ) ^{-n}\,{\rm d}{\it \_f}+ya-b\int \! \left ( {\frac {\sin \left ( \mu \,x \right ) }{\cos \left ( \mu \,x \right ) }} \right ) ^{m} \left ( {\frac {\sin \left ( \lambda \,x \right ) }{\cos \left ( \lambda \,x \right ) }} \right ) ^{-n}\,{\rm d}x \right ) } \right ) \left ( \cos \left ( {\frac {\nu }{a} \left ( b\int \! \left ( {\frac {\sin \left ( {\it \_f}\,\mu \right ) }{\cos \left ( {\it \_f}\,\mu \right ) }} \right ) ^{m} \left ( {\frac {\sin \left ( \lambda \,{\it \_f} \right ) }{\cos \left ( \lambda \,{\it \_f} \right ) }} \right ) ^{-n}\,{\rm d}{\it \_f}+ya-b\int \! \left ( {\frac {\sin \left ( \mu \,x \right ) }{\cos \left ( \mu \,x \right ) }} \right ) ^{m} \left ( {\frac {\sin \left ( \lambda \,x \right ) }{\cos \left ( \lambda \,x \right ) }} \right ) ^{-n}\,{\rm d}x \right ) } \right ) \right ) ^{-1}} \right ) ^{k} \left ( {\frac {\sin \left ( \lambda \,{\it \_f} \right ) }{\cos \left ( \lambda \,{\it \_f} \right ) }} \right ) ^{-n}\,{\rm d}{\it \_f}}}}}{d{\it \_f}}+{\it \_F1} \left ( {\frac {1}{a} \left ( ya-b\int \! \left ( {\frac {\sin \left ( \mu \,x \right ) }{\cos \left ( \mu \,x \right ) }} \right ) ^{m} \left ( {\frac {\sin \left ( \lambda \,x \right ) }{\cos \left ( \lambda \,x \right ) }} \right ) ^{-n}\,{\rm d}x \right ) } \right ) \right ) \]