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}])
Out[]=
Out[]=
Out[]=
Out[]=
Out[]=
