In[]:=
Clear["Global`*"]p[x_]:=A1*Exp[-x/l]+A2*Exp[x/l];u[x_]:=A3+A4*x-(h/c)*(A1*Exp[-x/l]+A2*Exp[x/l]);eqn3=μ''[x]+(q/ϵ)*p'[x]-μ[x]/λ^2+(1/λ^2)*(μhat+α*u'[x])==0;μ=FullSimplify[DSolveValue[eqn3,μ,x]];eqn9=c*u'[(w/2)]+h*p'[(w/2)]+(α/β)*(μ[(w/2)]-μhat-α*u'[(w/2)])==σ;eqn4=h*u'[(w/2)]+b*p'[(w/2)]==0;eqn5=h*u'[-(w/2)]+b*p'[-(w/2)]==0;eqn6=u[0]==0;eqn7=μ[(w/2)]==μ1;eqn8=μ[-(w/2)]==μ0;{a1,a2,a3,a4,a5,a6}=SolveValues[{eqn9,eqn4,eqn5,eqn6,eqn7,eqn8},{A1,A2,A3,A4,C[2],C[1]}];mu=μ/.{A1->a1,A2->a2,A3->a3,A4->a4,C[1]->a5,C[2]->a6};c[x_]:=c0+(1/β)*(mu[x]-μhat-α*u'[x]);c[x]K=1/((1/w)*Integrate[1/(γ*c[x]),{x,-(w/2),-(w/2)}])ComplexExpand[K]
w
2l
(1+)b-bcβ+β
Power[2]w
2
α
2
h
w
2l
(1+)b-bcβ+β
1
2
α
2
h
1
27
--w(Times[2]+Power[2])-
w
2l
1
2
w(l+λ)
2lλ
Times[3]+3
2
l
2
α
2c(1+)(-1+)-b+bcβ-βϵ(l-λ)(l+λ)
Power[2]w
2wPower[2]
2
α
2
h
Out[]=
c0+a6+a5-A4-α-μhat+a2+hlαϵ-a1hlαϵ-a1hlαϵ+a2+x-++hlαϵ-a4c+αϵ-a4c+x-++αϵ+a2+hαϵλ+a1hαϵλ-a1hαϵλ-a2+x-++hαϵλ-a2c+lq+a1clq+a1clq-a2c+x-++lq+a4c+αϵ+a4c+x-++αϵ-a2c+q-a1cq+a1cq+a2c+x-++q-c+ϵμhat-c+x-++ϵμhat+c+ϵμhat+c+x-++ϵμhat
1
β
-
x
λ
x/λ
h-+
A1
-
x
l
l
A2
x
l
l
c
1
2cϵ(-l+λ)(l+λ)
--
x
λ
x(l+λ)
lλ
2x
l
2x
λ
2x
λ
x-++
1
l
1
λ
x(l+λ)
lλ
2x
l
1
l
1
λ
x(l+λ)
lλ
x
l
2x
λ
2
l
x
l
1
l
1
λ
x(l+λ)
lλ
2
l
2x
l
2x
λ
2x
λ
x-++
1
l
1
λ
x(l+λ)
lλ
2x
l
1
l
1
λ
x(l+λ)
lλ
2x
l
2x
λ
2
λ
2x
λ
2
λ
x-++
1
l
1
λ
x(l+λ)
lλ
2
λ
2x
l
1
l
1
λ
x(l+λ)
lλ
2
λ
x
l
2x
λ
2
λ
x
l
1
l
1
λ
x(l+λ)
lλ
2
λ
2x
l
2x
λ
3
λ
2x
λ
3
λ
x-++
1
l
1
λ
x(l+λ)
lλ
3
λ
2x
l
1
l
1
λ
x(l+λ)
lλ
3
λ
x
l
2x
λ
2
l
x
l
1
l
1
λ
x(l+λ)
lλ
2
l
x
l
2x
λ
2
λ
x
l
1
l
1
λ
x(l+λ)
lλ
2
λ
1
0
Out[]=
ComplexInfinity
Out[]=
ComplexInfinity
