I want like this plots here two phase are there boundary is -2 to0 and 0 to 2

f(±2)=0,g(±2)=0,f(0)=1,h(±2)=1,H(±2)=1

How to integrate the below term

int((n*p+exp(n*p))*exp(-p)*p^(s-1)*f^s/((n*p+exp(2*n*p))*factorial(s)), p = 0 .. 1);
eval(%);

det_roots.docx

How to find the roots of the equation ?

like if the roots of the equations are m ,n,o,p then the value of these roots in terms of λ₀  λ₁  λ₂  λ₃  µ1  µ2 and  µ3

Please Help me to solve this problem 
 

Download ode_Plots_error.mw

P1=(D1*ga+C1)/(1+ga) and P2= (D1-C1)/(1+ga), How to substitute in the solution of U_exact and Theta_exact

Please help me to solve