{VERSION 4 0 "IBM INTEL NT" "4.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 1 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 1 }{CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 1 } {PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple Output" -1 11 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }3 3 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple Output" -1 12 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 3 0 0 0 0 1 0 1 0 2 2 0 1 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 2097 "#Get the series so lutions in tight coupling\n#Start with the vector perturbation equatio ns, in A=0 gauge\n#See astro-ph/0403583\n#AL March 04\n\n# Use conform al time. \n# Includes cdm, baryons, lambda\n\nrestart;\nno_quint:=\{ps i(t)=0,V(t)=0,V1(t)=0,V2(t)=0,rhopsi(t)=0,ppsi(t)=0\};\nno_numassive:= \{pinu(t)=0,qnu(t)=0,rhonu(t)=0,pnu(t)=0\};\nassign(no_quint);\nassign (no_numassive);\n\nH_t:=diff(S(t),t)/S(t);\n\ndS:=S(t)*H(t);\ndH:=-1/6 *S(t)^2*kappa*(rho(t)+3*p(t));\n\nK:=0;\nFriedmann:=H(t)^2=1/3*S(t)^2* kappa*rho(t)-K;\n\nphi(t):=-kappa*S(t)^2/2/k^3*(k*rhopi(t) + 4 * H(t) \+ * rhoq(t));\n\ndrag_t:=opac(t)*(4/3*v(t)-qg(t));\nphotbar_t:=rhog(t)/r hob(t);\n\ndqr:=-1/2*k*pir(t);\ndqg:=-1/2*k*pig(t)+drag(t);\ndv:=-(1-3 *c2(t))*H(t)*v(t)-photbar(t)*drag(t) - k/2*B0*rhog(t)/rhob(t);\n\n\n#d sigma:=simplify(-H(t)*sigma(t)+k*phi(t)-1/2*kappa*S(t)^2/k*(rhopi(t))) ;\ndsigma:=simplify(-2*H(t)*sigma(t)-kappa*S(t)^2/k*(rhopi(t)));\n\n\n rho_t:=rhob(t)+rhoc(t)+rhor(t)+rhog(t)+rhonu(t)+rhopsi(t)+rhov(t);\np_ t:=1/3*(rhor(t)+rhog(t))+pb(t)+pnu(t)+ppsi(t)-rhov(t);\ndpb:=c2(t)*drh ob;\n\n\ndpig:=-opac(t)*(pig(t)-polter(t)) - 8/15*k*J_3(t) + 2/5*k*qg( t) + 8/15*k*sigma(t);\ndJ_3:=k*(3/7*pig(t)-15/28*J_4(t))-opac(t)*J_3(t );\n\ndpir:=- 8/15*k*G_3(t) + 2/5*k*qr(t) + 8/15*k*sigma(t);\ndG_3:=k* (3/7*pir(t)-15/28*G_4(t));\ndG_4:=k*(4/9*G_3(t)-24/45*Kf[4]*G_5(t));\n \nrhopi_t:=rhog(t)*pig(t)+rhor(t)*pir(t) + rhonu(t)*pinu(t) + rhog(t)* B0;\nrhoq_t:=rhog(t)*qg(t)+rhor(t)*qr(t)+(rhob(t)+pb(t))*v(t) + rhonu( t)*qnu(t) +k*diff(psi(t),t)*clv(t)/S(t)^2;\n\nsigma(t):=2*kappa*S(t)^2 *rhoq_t/k^2;\n#rhoq(t):=k^2*sigma(t)/2/kappa/S(t)^2;\n\n\nsubtots:=\{r hopi(t)=rhopi_t,rho(t)=rho_t,p(t)=p_t\};\n\ndrhoq:=-4*H(t)*rhoq(t)-1/2 *k*rhopi(t);\n\npolter_t:=2/15*(3*pig(t)/4 + 9*E2(t)/2);\n\ndE2:=-opac (t)*(E2(t) - polter(t)) - 8/27*k*E3(t) + 1/3*k*B2(t);\ndB2:=-opac(t)*B 2(t) - 8/27*k*B3(t) - 1/3*k*E2(t);\n\nsublist:=\{diff(S(t),t)=dS,diff (qr(t),t)=dqr,diff(qg(t),t)=dqg,diff(v(t),t)=dv,diff(H(t),t)=dH,diff(e xptau(t),t)=g(t),diff(pig(t),t)=dpig\};\n\nsubtotderivs:=\{diff(S(t),t )=dS,diff(H(t),t)=dH\};\n\n#End of main definitions and equations\n### ######################################\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n \n\n\n\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%)no_quintG<(/-%$psiG6#% \"tG\"\"!/-%\"VGF)F+/-%#V1GF)F+/-%#V2GF)F+/-%'rhopsiGF)F+/-%%ppsiGF)F+ " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%-no_numassiveG<&/-%%pinuG6#%\"tG \"\"!/-%$qnuGF)F+/-%&rhonuGF)F+/-%$pnuGF)F+" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$H_tG*&-%%diffG6$-%\"SG6#%\"tGF,\"\"\"F)!\"\"" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%#dSG*&-%\"SG6#%\"tG\"\"\"-%\"HGF(F* " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#dHG,$*()-%\"SG6#%\"tG\"\"#\"\" \"%&kappaGF-,&-%$rhoGF*F-*&\"\"$F--%\"pGF*F-F-F-#!\"\"\"\"'" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"KG\"\"!" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%*FriedmannG/*$)-%\"HG6#%\"tG\"\"#\"\"\",$*()-%\"SGF*F,F-%&kappaGF --%$rhoGF*F-#F-\"\"$" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>-%$phiG6#%\"t G,$*&*(%&kappaG\"\"\")-%\"SGF&\"\"#F,,&*&%\"kGF,-%&rhopiGF&F,F,*(\"\"% F,-%\"HGF&F,-%%rhoqGF&F,F,F,F,*$)F3\"\"$F,!\"\"#F?F0" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'drag_tG*&-%%opacG6#%\"tG\"\"\",&-%\"vGF(#\"\"%\" \"$-%#qgGF(!\"\"F*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%*photbar_tG*&- %%rhogG6#%\"tG\"\"\"-%%rhobGF(!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6 #>%$dqrG,$*&%\"kG\"\"\"-%$pirG6#%\"tGF(#!\"\"\"\"#" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>%$dqgG,&*&%\"kG\"\"\"-%$pigG6#%\"tGF(#!\"\"\"\"#-%%d ragGF+F(" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#dvG,(*(,&\"\"\"F(*&\"\" $F(-%#c2G6#%\"tGF(!\"\"F(-%\"HGF-F(-%\"vGF-F(F/*&-%(photbarGF-F(-%%dra gGF-F(F/*&#F(\"\"#F(*&*(%\"kGF(%#B0GF(-%%rhogGF-F(F(-%%rhobGF-F/F(F/" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'dsigmaG,$*&,&*(-%\"HG6#%\"tG\"\" \"-%&sigmaGF+F-%\"kGF-\"\"#*(%&kappaGF-)-%\"SGF+F1F--%&rhopiGF+F-F-F-F 0!\"\"F9" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&rho_tG,,-%%rhobG6#%\"tG \"\"\"-%%rhocGF(F*-%%rhorGF(F*-%%rhogGF(F*-%%rhovGF(F*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$p_tG,*-%%rhorG6#%\"tG#\"\"\"\"\"$*&F*F+-%%rhogG F(F+F+-%#pbGF(F+-%%rhovGF(!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$ dpbG*&-%#c2G6#%\"tG\"\"\"%&drhobGF*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6 #>%%dpigG,**&-%%opacG6#%\"tG\"\"\",&-%$pigGF)F+-%'polterGF)!\"\"F+F1*& #\"\")\"#:F+*&%\"kGF+-%$J_3GF)F+F+F1*(#\"\"#\"\"&F+F7F+-%#qgGF)F+F+*(# F4F5F+F7F+-%&sigmaGF)F+F+" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%dJ_3G, &*&%\"kG\"\"\",&-%$pigG6#%\"tG#\"\"$\"\"(*&#\"#:\"#GF(-%$J_4GF,F(!\"\" F(F(*&-%%opacGF,F(-%$J_3GF,F(F7" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%% dpirG,(*&%\"kG\"\"\"-%$G_3G6#%\"tGF(#!\")\"#:*(#\"\"#\"\"&F(F'F(-%#qrG F+F(F(*(#\"\")F/F(F'F(-%&sigmaGF+F(F(" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%dG_3G*&%\"kG\"\"\",&-%$pirG6#%\"tG#\"\"$\"\"(*&#\"#:\"#GF'-%$G_4 GF+F'!\"\"F'" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%dG_4G*&%\"kG\"\"\", &-%$G_3G6#%\"tG#\"\"%\"\"**&#\"\")\"#:F'*&&%#KfG6#F.F'-%$G_5GF+F'F'!\" \"F'" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(rhopi_tG,(*&-%%rhogG6#%\"tG \"\"\"-%$pigGF)F+F+*&-%%rhorGF)F+-%$pirGF)F+F+*&F'F+%#B0GF+F+" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%'rhoq_tG,(*&-%%rhogG6#%\"tG\"\"\"-%# qgGF)F+F+*&-%%rhorGF)F+-%#qrGF)F+F+*&,&-%%rhobGF)F+-%#pbGF)F+F+-%\"vGF )F+F+" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>-%&sigmaG6#%\"tG,$*&*(%&kapp aG\"\"\")-%\"SGF&\"\"#F,,(*&-%%rhogGF&F,-%#qgGF&F,F,*&-%%rhorGF&F,-%#q rGF&F,F,*&,&-%%rhobGF&F,-%#pbGF&F,F,-%\"vGF&F,F,F,F,*$)%\"kGF0F,!\"\"F 0" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%(subtotsG<%/-%&rhopiG6#%\"tG,(* &-%%rhogGF)\"\"\"-%$pigGF)F/F/*&-%%rhorGF)F/-%$pirGF)F/F/*&F-F/%#B0GF/ F//-%$rhoGF),,-%%rhobGF)F/-%%rhocGF)F/F3F/F-F/-%%rhovGF)F//-%\"pGF),*F 3#F/\"\"$*&FGF/F-F/F/-%#pbGF)F/FA!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&drhoqG,&*&-%\"HG6#%\"tG\"\"\"-%%rhoqGF)F+!\"%*&#F+\"\"#F+*&% \"kGF+-%&rhopiGF)F+F+!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%)polte r_tG,&-%$pigG6#%\"tG#\"\"\"\"#5*&#\"\"$\"\"&F+-%#E2GF(F+F+" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$dE2G,(*&-%%opacG6#%\"tG\"\"\",&-%#E2GF)F+ -%'polterGF)!\"\"F+F1*&#\"\")\"#FF+*&%\"kGF+-%#E3GF)F+F+F1*(#F+\"\"$F+ F7F+-%#B2GF)F+F+" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$dB2G,(*&-%%opac G6#%\"tG\"\"\"-%#B2GF)F+!\"\"*&#\"\")\"#FF+*&%\"kGF+-%#B3GF)F+F+F.*&#F +\"\"$F+*&F4F+-%#E2GF)F+F+F." }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%(sub listG<)/-%%diffG6$-%'exptauG6#%\"tGF--%\"gGF,/-F(6$-%\"HGF,F-,$*()-%\" SGF,\"\"#\"\"\"%&kappaGF;,&-%$rhoGF,F;*&\"\"$F;-%\"pGF,F;F;F;#!\"\"\" \"'/-F(6$-%$pigGF,F-,**&-%%opacGF,F;,&FJF;-%'polterGF,FEF;FE*&#\"\")\" #:F;*&%\"kGF;-%$J_3GF,F;F;FE*(#F:\"\"&F;FXF;-%#qgGF,F;F;*&**#\"#;FVF;F %-subtotderivsG<$/-%%diffG6$-%\"HG6#%\"tGF-,$*()-%\"SGF,\"\"#\"\"\"% &kappaGF4,&-%$rhoGF,F4*&\"\"$F4-%\"pGF,F4F4F4#!\"\"\"\"'/-F(6$F1F-*&F1 F4F*F4" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 1952 "#Get series sol ution using tight coupling result:\nOrd:=4;\nDoQuint = \"F\";\n\nOrder :=Ord+5;\n\nmakeseries:=proc(v,x,ord)\n local Res,i;\n Res:=0;\n for i from 0 to ord do\n Res:=Res + cat('s',v)[i]*x^i;\n end do;\n v(x)=s eries(Res,x,ord); \nend;\n \n\nexpandfully:=proc(x)\n eval(x);\n %; \n %;\nend;\n\n\nRb:=1-Rc;\nv(t):=3/4*qg(t);\n\nomb:=Rb*omm;\nomc:=Rc* omm;\nomr:=Rv*omm^2*H0^2/omega^2;\nomg:=omm^2*H0^2/omega^2-omr;\n\npb( t):=0;c2(t):=0;\npig(t):=0;\n\n\nrhog(t):=3*omg*H0^2/kappa/S(t)^4;\nrh or(t):=3*omr*H0^2/kappa/S(t)^4;\nrhob(t):=3*omb*H0^2/kappa/S(t)^3;\nrh oc(t):=3*omc*H0^2/kappa/S(t)^3;\nrhov(t):=3*omv*H0^2/kappa;\nH(t):=H_t ;\n\nassign(subtots);\n\nK:=0;\nKfac:=beta[2];\n\n\n\n#Start series so lving\ndsolve(\{expandfully(Friedmann),S(0)=0,D(S)(0)=H0^2*omm/omega\} ,S(t),type=series);\nassign(%);\n\n# S(t):=omm*H0^2/omega^2*(omega*t+( omega*t)^2/4-K/6*omega*t^3-K/48*omega^2*t^4);\n\n\nRb:=1-Rc;\n\nv(t):= 3/4*qg(t);\n\nomb:=Rb*omm;\nomc:=Rc*omm;\nomr:=Rv*omm^2*H0^2/omega^2; \nomg:=omm^2*H0^2/omega^2-omr;\n\npb(t):=0;c2(t):=0;\npig(t):=0;\n\n#U se tight coupling result\nR:=4/3*rhog(t)/rhob(t);\ndrag(t):=1/(1+R)*(- 4/3*H(t)*v(t) - k/2*B0*R);\n\n#Some frame invariant variables\n\nsolv evars:= \{qr,qg,pir,G_3,G_4\};\n\neqs:=\{simplify(dqr-diff(qr(t),t)),s implify(dqg-diff(qg(t),t)),simplify(diff(pir(t),t)-dpir),simplify(diff (G_3(t),t)-dG_3),simplify(diff(G_4(t),t)-dG_4)\};\n\n\nknownvars:=\{H0 ,omm,omv,k,beta[2],Kf[2],Kf[3],Kf[4],omega,Rv,Rc,B0,sqg[0]\};\n\nserie s_subs:=map(makeseries,solvevars,t,Ord+3);\nassign(%);\n\n#spir[0]:=0; \n#sG_3[1]:=0;\n#sG_4[2]:=0;\n\nG_5(t):=sG_5[3]*t^3+sG_5[4]*t^4+sG_ 5*t^5;\nsG_3[0]:=0;\nsG_4[0]:=0;\nsG_4[1]:=0;\n\n\nseries_eqs:=map(exp andfully,eqs);\n\nmap(series,series_eqs,t):\nmap(simplify,%):\nseries_ eqs:=%;\n\ndosolve:=proc(ineqs)\n local tmp,vars;\n map(series,ineqs,t );\n map(simplify,%);\n map(convert,%,polynom);\n map(tcoeff,%,t):\n t mp:=map(simplify,%);\n vars:=indets(tmp) minus knownvars;\n print(vars );\n solve(tmp,vars);\n assign(%);\nend;\n\n\n\nfor i from 1 to Ord+1 \+ do\n dosolve(series_eqs);\nend do;\n" }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#>%$OrdG\"\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/%(DoQuintGQ\"F6\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&OrderG\"\"*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%+makeseriesGR6%%\"vG%\"xG%$ordG6$%$ResG%\"iG6\"F-C%>8 $\"\"!?(8%F1\"\"\"9&%%trueG>F0,&F0F4*&&-%$catG6$.%\"sG9$6#F3F4)9%F3F4F 4/-F@6#FC-%'seriesG6%F0FCF5F-F-F-" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#> %,expandfullyGR6#%\"xG6\"F(F(C%-%%evalG6#9$%\"%GF.F(F(F(" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>%#RbG,&\"\"\"F&%#RcG!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>-%\"vG6#%\"tG,$-%#qgGF&#\"\"$\"\"%" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>%$ombG*&,&\"\"\"F'%#RcG!\"\"F'%$ommGF'" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>%$omcG*&%#RcG\"\"\"%$ommGF'" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>%$omrG*&*(%#RvG\"\"\")%$ommG\"\"#F()%#H0GF+F(F(*$)%& omegaGF+F(!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$omgG,&*&*&)%$omm G\"\"#\"\"\")%#H0GF*F+F+*$)%&omegaGF*F+!\"\"F+*&*(%#RvGF+F(F+F,F+F+*$F /F+F1F1" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>-%#pbG6#%\"tG\"\"!" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>-%#c2G6#%\"tG\"\"!" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>-%$pigG6#%\"tG\"\"!" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>-%%rhogG6#%\"tG,$*&*&,&*&*&)%$ommG\"\"#\"\"\")%#H0GF0F1F1*$)%&omeg aGF0F1!\"\"F1*&*(%#RvGF1F.F1F2F1F1*$F5F1F7F7F1F2F1F1*&%&kappaGF1)-%\"S GF&\"\"%F1F7\"\"$" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>-%%rhorG6#%\"tG, $*&*(%#RvG\"\"\")%$ommG\"\"#F,)%#H0G\"\"%F,F,*()%&omegaGF/F,%&kappaGF, )-%\"SGF&F2F,!\"\"\"\"$" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>-%%rhobG6# %\"tG,$*&*(,&\"\"\"F,%#RcG!\"\"F,%$ommGF,)%#H0G\"\"#F,F,*&%&kappaGF,)- %\"SGF&\"\"$F,F.F8" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>-%%rhocG6#%\"tG ,$*&*(%#RcG\"\"\"%$ommGF,)%#H0G\"\"#F,F,*&%&kappaGF,)-%\"SGF&\"\"$F,! \"\"F6" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>-%%rhovG6#%\"tG,$*&*&%$omvG \"\"\")%#H0G\"\"#F,F,%&kappaG!\"\"\"\"$" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>-%\"HG6#%\"tG*&-%%diffG6$-%\"SGF&F'\"\"\"F,!\"\"" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>%\"KG\"\"!" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#> %%KfacG&%%betaG6#\"\"#" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#/-%\"SG6#%\" tG+1F'*&*&)%#H0G\"\"#\"\"\"%$ommGF.F.%&omegaG!\"\"F.,$*&F+F.F/F.#F.\" \"%F-,$*&*()F,\"\")F.%$omvGF.)F/\"\"$F.F.*$)F0F=F.F1#F.\"#5\"\"&,$*&*( F9F.F;F.F%#RbG,&\"\"\"F&%#RcG!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>-%\"vG6#%\"tG,$-%#qgGF&#\"\"$\"\"%" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>%$ombG*&,&\"\"\"F'%#RcG!\"\"F'%$ommGF'" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>%$omcG*&%#RcG\"\"\"%$ommGF'" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>%$omrG*&*(%#RvG\"\"\")%$ommG\"\"#F()%#H0GF+F(F(*$)%& omegaGF+F(!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$omgG,&*&*&)%$omm G\"\"#\"\"\")%#H0GF*F+F+*$)%&omegaGF*F+!\"\"F+*&*(%#RvGF+F(F+F,F+F+*$F /F+F1F1" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>-%#pbG6#%\"tG\"\"!" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>-%#c2G6#%\"tG\"\"!" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>-%$pigG6#%\"tG\"\"!" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"RG,$*&,&*&*&)%$ommG\"\"#\"\"\")%#H0GF,F-F-*$)%&omegaGF,F-!\"\"F -*&*(%#RvGF-F*F-F.F-F-*$F1F-F3F3F-*(-%\"SG6#%\"tGF-,&F-F-%#RcGF3F-F+F- F3#\"\"%\"\"$" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>-%%dragG6#%\"tG*&,&* &*&-%%diffG6$-%\"SGF&F'\"\"\"-%#qgGF&F1F1F/!\"\"F4*&#\"\"#\"\"$F1*&*(% \"kGF1%#B0GF1,&*&*&)%$ommGF7F1)%#H0GF7F1F1*$)%&omegaGF7F1F4F1*&*(%#RvG F1F@F1FBF1F1*$FEF1F4F4F1F1*(+1F'*&*&FBF1FAF1F1FFF4F1,$*&FBF1FAF1#F1\" \"%F7,$*&*()FC\"\")F1%$omvGF1)FAF8F1F1*$)FFF8F1F4#F1\"#5\"\"&,$*&*(FVF 1FXF1FYF1F1*$FEF1F4#F1\"#?\"\"',$*&*(FVF1FXF1FYF1F1FFF4#F1\"$7\"\"\"(, $*(FVF1FXF1FYF1#F1\"%#z\"FW-%\"OG6#F1\"\"*F1,&F1F1%#RcGF4F1FAF1F4F1F4F 1,&F1F1*&*&#FRF8F1F=F1F1*(FLF1F^pF1FAF1F4F1F4" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%*solvevarsG<'%$G_4G%#qrG%$G_3G%#qgG%$pirG" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%$eqsG<',$*&,2*(+1%\"tG*&*&)%#H0G\"\"#\"\" \"%$ommGF1F1%&omegaG!\"\"\"\"!,$*&F.F1F2F1#F1F0F1,$*&*()F/\"\")F1%$omv GF1)F2\"\"$F1F1*$)F3F@F1F4F8\"\"%,$*&*(FF1F?F1F1*$)F3F0F1F4#F@\"# 5\"\"&,$*&*(FF1F?F1F1F3F4#F1\"#;\"\"',$*(FF1F?F1#F1\"$C#\"\" (-%\"OG6#F1F=F1-%#qgG6#F+F1FHF1!\"$*,F@F1F*F1FZF1FHF1%#RcGF1F1*,F0F1% \"kGF1%#B0GF1F2F1F.F1F4*.F0F1F[oF1F\\oF1F2F1F.F1%#RvGF1F1**F@F1-%%diff G6$FZF+F1FHF1+1F+F,F1,$F7#F1FCF0,$F:#F1FJFK,$FE#F1\"#?FQ,$FM#F1\"$7\"F V,$FS#F1\"%#z\"F=FW\"\"*F1F4*,F@F1F`oF1FHF1FcoF1FinF1F1**FCF1F`oF1F.F1 F2F1F4*,FCF1F`oF1F.F1F2F1F^oF1F1F1,**&FHF1FcoF1Fgn**F@F1FHF1FcoF1FinF1 F1*(FCF1F.F1F2F1F4**FCF1F.F1F2F1F^oF1F1F4F4,&*&F[oF1-%$pirGFfnF1#F4F0- Fao6$-%#qrGFfnF+F4,(-Fao6$-%$G_3GFfnF+F1*&#F@FVF1F[qF1F4*(#\"#:\"#GF1F [oF1-%$G_4GFfnF1F1,(-Fao6$F^rF+F1*&#FCFapF1*&F[oF1FfqF1F1F4**#F=F\\rF1 F[oF1&%#KfG6#FCF1-%$G_5GFfnF1F1,$*&,2**-Fao6$F\\qF+F1F[oF1FHF1)FcoF0F1 F\\r*,F=F1)F[oF0F1FfqF1FHF1FcsF1F1*,FQF1FesF1FaqF1FHF1FcsF1F4**\"#[F1) F/FCF1)F2F0F1FZF1F4*,FhsF1FisF1FjsF1FZF1F^oF1F1*,FhsF1F^oF1FjsF1FisF1F aqF1F4*.\"#OF1F.F1F2F1FZF1FHF1FcoF1F4*0F^tF1F.F1F2F1FZF1FHF1FcoF1FinF1 F1F1*(FcsF1F[oF1FHF1F4#F1F\\r" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%*kn ownvarsG%,series_subsG<'/-%#qgG6#%\"tG+3F*&%$sqgG6#\"\"!F/&F-6 #\"\"\"F2&F-6#\"\"#F5&F-6#\"\"$F8&F-6#\"\"%F;&F-6#\"\"&F>&F-6#\"\"'FA- %\"OGF1\"\"(/-%$pirGF)+3F*&%%spirGF.F/&FJF1F2&FJF4F5&FJF7F8&FJF:F;&FJF =F>&FJF@FAFBFD/-%#qrGF)+3F*&%$sqrGF.F/&FVF1F2&FVF4F5&FVF7F8&FVF:F;&FVF =F>&FVF@FAFBFD/-%$G_3GF)+3F*&%%sG_3GF.F/&F\\oF1F2&F\\oF4F5&F\\oF7F8&F \\oF:F;&F\\oF=F>&F\\oF@FAFBFD/-%$G_4GF)+3F*&%%sG_4GF.F/&FhoF1F2&FhoF4F 5&FhoF7F8&FhoF:F;&FhoF=F>&FhoF@FAFBFD" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>-%$G_5G6#%\"tG,(*&&%%sG_5G6#\"\"$\"\"\")F'F-F.F.*&&F+6#\"\"%F.)F'F 3F.F.*&F+F.)F'\"\"&F.F." }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%%sG_3G6# \"\"!F'" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%%sG_4G6#\"\"!F'" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>&%%sG_4G6#\"\"\"\"\"!" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%+series_eqsG<',(+1%\"tG&%%sG_3G6#\"\"\"\"\"!,$&F *6#\"\"#F1F,,$&F*6#\"\"$F5F1,$&F*6#\"\"%F9F5,$&F*6#\"\"&F=F9,$&F*6#\" \"'FAF=-%\"OGF+FAF,*&#F5\"\"(F,*&%\"kGF,+3F(&%%spirG6#F-F-&FKF+F,&FKF0 F1&FKF4F5&FKF8F9&FKF%+series_eqsG<'+/%\"tG,&&%%sG_4G6#\"\"#F,*&#\"\"%\"\"* \"\"\"*&%\"kGF1&%%sG_3G6#F1F1F1!\"\"F1,&&F*6#\"\"$F;*&#F/F0F1*&F3F1&F5 F+F1F1F7F,,(*&F3F1&F5F:F1#!\"%F0**#\"\")\"#:F1F3F1&%#KfG6#F/F1&%%sG_5G F:F1F1*&F/F1&F*FKF1F1F;,(&F*6#\"\"&FS*&#F/F0F1*&F3F1&F5FKF1F1F7**FFF1F 3F1FIF1&FMFKF1F1F/,(&F*6#\"\"'Fgn**FFF1F3F1FIF1FMF1F1*&#F/F0F1*&F3F1&F 5FRF1F1F7FS-%\"OGF6Fgn+1F',$*&,.*&%&omegaGF1&%$sqgG6#\"\"!F1!\"$*(F/F1 &FfoF6F1%#RvGF1F1*&F/F1F[pF1F7*(F,F1F3F1%#B0GF1F7**F,F1F3F1F_pF1F\\pF1 F1**F;F1FdoF1FeoF1%#RcGF1F1F1,&F7F1F\\pF1F7#F7F/Fho,$*&,B*&&FfoF+F1F\\ pF1\"#k*&\"#KF1FipF1F7*.FgnF1FdoF1F3F1F_pF1F\\pF1FbpF1F1**\"#7F1F[pF1F doF1FbpF1F1*,F_qF1F[pF1FdoF1FbpF1F\\pF1F7*,FgnF1)FdoF,F1FbpF1FeoF1F\\p F1F7*,FgnF1FdoF1F3F1F_pF1FbpF1F7*,FgnF1FdoF1F3F1F_pF1F\\pF1F7*(F_qF1Fd oF1F[pF1F7*(F;F1FbqF1FeoF1F1*(F\\qF1FipF1)F\\pF,F1F7**FgnF1FbqF1FeoF1F \\pF1F1**F_qF1FdoF1F[pF1F\\pF1F1**F_qF1FbqF1FbpF1FeoF1F7**F0F1FbqF1Feo F1)FbpF,F1F1**FgnF1FdoF1F3F1F_pF1F1F1*$)FcpF,F1F7#F1\"#;F1,$*&,\\o*()F doF;F1FeoF1F\\pF1!\"**&FjpF1&FfoF:F1F7*(F/F1F[pF1FbqF1F1*.FgnF1FbqF1F3 F1F_pF1F\\pF1F]rF1F1**F0F1FgrF1FeoF1FbpF1F1*.FGF1FbqF1F3F1F_pF1F\\pF1F bpF1F7*.F,F1FbqF1F3F1F_pF1FhqF1FbpF1F7**F/F1FbqF1F[pF1F\\pF1F1**FGF1Fb qF1F[pF1FhqF1F7**F\\qF1FdoF1FipF1F\\pF1F1**FbrF1FdoF1FipF1FhqF1F7**Fbr F1FbqF1FbpF1F[pF1F7**FbrF1FdoF1FbpF1FipF1F1**F/F1FbqF1F3F1F_pF1F7**\"# =F1FgrF1FeoF1F]rF1F7**F_qF1F[pF1FbqF1F]rF1F1**F0F1FgrF1FeoF1)FbpF;F1F1 *,F\\qF1FdoF1FbpF1FipF1F\\pF1F7*,FbrF1FdoF1FbpF1FipF1FhqF1F1*,FGF1FbqF 1FbpF1F[pF1F\\pF1F1*,FGF1FbqF1FbpF1F[pF1FhqF1F1*,FhsF1FgrF1FeoF1FbpF1F \\pF1F1*,F,F1FbqF1F3F1F_pF1F\\pF1F1*,\"#5F1FbqF1F3F1F_pF1FbpF1F1*,F,F1 FbqF1F3F1F_pF1FhqF1F1*,F0F1FgrF1FeoF1F]rF1F\\pF1F7*,F_qF1F[pF1FbqF1F]r F1F\\pF1F7*,FgnF1FbqF1F3F1F_pF1F]rF1F7*(FbrF1FdoF1FipF1F7*(\"$#>F1FjrF 1FhqF1F7*(FjpF1FjrF1)F\\pF;F1F1*(FjtF1FjrF1F\\pF1F1F1*$)FcpF;F1F7#FioF jpF,,$*&,jq&FfoFK!%C5*(\"#[F1FipF1FbqF1F1*(F0F1)FdoF/F1FeoF1F7*(FjtF1F doF1FjrF1F7*,\"$3\"F1F[pF1FgrF1F[tF1F\\pF1F7*,\"#aF1FgrF1F3F1F_pF1F[tF 1F7*,F\\vF1FgrF1F[pF1FbpF1F\\pF1F1*,\"$_#F1FiuF1FbpF1FeoF1F\\pF1F7*,\" #!*F1FgrF1F3F1F_pF1FbpF1F7*,FhsF1FgrF1F3F1F_pF1F\\pF1F1*,\"$)GF1FbqF1F ipF1F\\pF1FbpF1F1*,\"$;#F1FgrF1F[pF1FbpF1FhqF1F7*,\"#OF1FiuF1FbpF1FeoF 1FhqF1F7*,FjvF1FgrF1F3F1F_pF1FhqF1F7*,F\\vF1FgrF1F[pF1F]rF1FhqF1F1*,Fh sF1FiuF1F]rF1FeoF1FhqF1F1*,F\\vF1FiuF1FeoF1F[tF1F\\pF1F7*,F\\vF1FgrF1F [pF1F]rF1F\\pF1F1*,FfvF1FiuF1F]rF1FeoF1F\\pF1F1*,\"$E\"F1FgrF1F3F1F_pF 1F]rF1F1*,\"#'*F1FbqF1FipF1F\\uF1FbpF1F7*,\"$w&F1FdoF1FbpF1FjrF1F\\pF1 F7*,FfwF1FdoF1FbpF1FjrF1FhqF1F1*,FjtF1FdoF1FbpF1FjrF1F\\uF1F7*.F^vF1Fg rF1F3F1F_pF1F\\pF1F[tF1F1*,FfvF1FbqF1F]rF1FipF1F\\pF1F7*,\"$W\"F1FbqF1 F]rF1FipF1FhqF1F1*(\"%WhF1FduF1FhqF1F7*(\"%'4%F1FduF1F\\uF1F1*(\"%C5F1 FduF1)F\\pF/F1F7*.FhsF1FgrF1F3F1F_pF1F\\pF1FbpF1F1*.FcvF1FgrF1F3F1F_pF 1F\\pF1F]rF1F7*.\"#sF1FgrF1F3F1F_pF1FhqF1FbpF1F1*.FjvF1FgrF1F3F1F_pF1F hqF1F]rF1F7**F\\vF1FgrF1F[pF1FhqF1F1**F\\vF1FgrF1F[pF1FbpF1F1**F\\xF1F bqF1FipF1FhqF1F7**FgxF1FiuF1FeoF1F\\pF1F1**F\\vF1FgrF1F[pF1F\\pF1F7**F jvF1FiuF1FbpF1FeoF1F7**\"$!=F1FiuF1FeoF1F]rF1F1**FhsF1FgrF1F3F1F_pF1F1 **FdwF1FbqF1FipF1F\\uF1F1**FhvF1FgrF1F[pF1F]rF1F7**FhvF1FiuF1FeoF1F[tF 1F7**FjtF1FbqF1FipF1FbpF1F7**FhsF1FiuF1FeoF1FhqF1F1**FfwF1FdoF1FjrF1F \\pF1F1**FfwF1FdoF1FjrF1FhqF1F7**FjtF1FdoF1FjrF1F\\uF1F1**FjtF1FdoF1Fb pF1FjrF1F1*(F`xF1FduF1F\\pF1F1**F\\xF1FbqF1F]rF1FipF1F1**F\\vF1F[pF1Fg rF1F[tF1F1**\"#\")F1FiuF1FeoF1)FbpF/F1F1F1*$)FcpF/F1F7#F1\"$c#F;,$*&,d w*()FdoFgnF1FeoF1F\\pF1!$N\"**\"%sIF1FduF1FbqF1F\\pF1F1*(\"$o(F1FduF1F bqF1F7*(FjtF1FjrF1FgrF1F1*,\"$%QF1)%#H0GFgnF1%$omvGF1)%$ommGF,F1FeoF1F 7**FfwF1FgrF1F]rF1FjrF1F1**\"$K%F1FiuF1F[tF1FipF1F1**\"$C$F1)FdoFSF1F[ pF1F`zF1F1**\"$V#F1FizF1FeoF1)FbpFSF1F1*,FfwF1FgrF1F]rF1FjrF1F\\uF1F7* ,F\\vF1F\\\\lF1F3F1F_pF1FhqF1F1*,\"%?;F1FizF1F]rF1FeoF1F\\pF1F7*(\"#FF 1FizF1FeoF1F1*.\"$i\"F1F\\\\lF1F3F1F_pF1FhqF1F[tF1F7*.FhsF1F\\\\lF1F3F 1F_pF1F\\uF1F]rF1F1*.FjvF1F\\\\lF1F3F1F_pF1F\\uF1FbpF1F7*,\"%'H\"F1Fiu F1FbpF1FipF1FhqF1F7*,\"$k)F1F\\\\lF1FbpF1F[pF1F\\pF1F7*,F]]lF1F\\\\lF1 FbpF1F[pF1FhqF1F1*,\"$5)F1FizF1FeoF1FbpF1F\\pF1F1*,F\\vF1F\\\\lF1F3F1F _pF1F\\pF1F7*,\"$)>F1F\\\\lF1F3F1F_pF1FbpF1F1*,Fi[lF1F\\\\lF1F[pF1F]rF 1F\\pF1F1*,\"$/&F1F\\\\lF1F3F1F_pF1F]rF1F7*,F[]lF1FiuF1F]rF1FipF1F\\pF 1F1*,\"%!3\"F1F\\\\lF1F]rF1F[pF1FhqF1F7*,\"%]8F1FizF1FeoF1F[tF1F\\pF1F 1*,Fi[lF1F\\\\lF1F[pF1F[tF1F\\pF1F1*,\"$'[F1F\\\\lF1F3F1F_pF1F[tF1F1*, \"$0%F1FizF1FeoF1F`zF1F\\pF1F7*,F]]lF1FiuF1FbpF1FipF1F\\uF1F1*,F\\xF1F \\\\lF1FbpF1F[pF1F\\uF1F1*,F`^lF1FizF1FeoF1FbpF1FhqF1F1*,FhsF1F\\\\lF1 F3F1F_pF1F\\uF1F1*,F`^lF1FizF1FeoF1F]rF1FhqF1F7*,Fi[lF1FiuF1F]rF1FipF1 F\\uF1F7*,FgxF1F\\\\lF1F]rF1F[pF1F\\uF1F7*,\"$N\"F1FizF1FeoF1F[tF1FhqF 1F1*,Fi[lF1F\\\\lF1F[pF1F[tF1FhqF1F1*.Fi[lF1F\\\\lF1F3F1F_pF1FhqF1F]rF 1F1*.F^vF1F\\\\lF1F3F1F_pF1F\\pF1F]rF1F1*.FhvF1F\\\\lF1F3F1F_pF1F\\pF1 FbpF1F1*.\"$y$F1F\\\\lF1F3F1F_pF1FhqF1FbpF1F7*.F[\\lF1F\\\\lF1F3F1F_pF 1F\\pF1F[tF1F7**FhvF1F\\\\lF1F[pF1FhqF1F7**Fi[lF1FiuF1FipF1F\\pF1F7**F ]]lF1FiuF1FipF1FhqF1F1**F\\xF1F\\\\lF1FbpF1F[pF1F7**Fi[lF1FiuF1FbpF1Fi pF1F1**FhsF1F\\\\lF1F3F1F_pF1F7**F`^lF1FizF1FeoF1F]rF1F7**\"$?(F1F\\\\ lF1F[pF1F]rF1F1**\"$X*F1FizF1FeoF1F[tF1F1**Fi[lF1FiuF1FipF1F\\uF1F7**F ]]lF1FiuF1F]rF1FipF1F7**F]]lF1F\\\\lF1F[pF1F[tF1F7**F`]lF1FizF1FeoF1F` zF1F7**Fi^lF1FizF1FeoF1FhqF1F7**FgxF1F\\\\lF1F[pF1F\\uF1F7**F[\\lF1F\\ \\lF1F[pF1F\\pF1F1*0Fa[lF1FbpF1Fb[lF1Fd[lF1Fe[lF1FeoF1FcxF1F1*.\"%O:F1 Fb[lF1Fd[lF1Fe[lF1FeoF1F\\uF1F1*.Fa[lF1Fb[lF1Fd[lF1Fe[lF1FeoF1FcxF1F7* .Fe`lF1Fb[lF1Fd[lF1Fe[lF1FeoF1F\\pF1F1*.\"%/BF1Fb[lF1Fd[lF1Fe[lF1FeoF1 FhqF1F7*.Fa[lF1FbpF1Fb[lF1Fd[lF1Fe[lF1FeoF1F1*0Fe`lF1FbpF1Fb[lF1Fd[lF1 Fe[lF1FeoF1F\\pF1F7*0Fi`lF1FbpF1Fb[lF1Fd[lF1Fe[lF1FeoF1FhqF1F1*0Fe`lF1 FbpF1Fb[lF1Fd[lF1Fe[lF1FeoF1F\\uF1F7**\"$g*F1FjrF1FgrF1F\\uF1F1**Fa[lF 1FjrF1FgrF1FcxF1F7**F\\[lF1FduF1FbqF1F\\uF1F1**F^[lF1FduF1FbqF1FcxF1F7 **FjtF1FjrF1FgrF1F\\pF1F7**\"&+7&F1&FfoFRF1FhqF1FdoF1F7**FealF1FfalF1F \\uF1FdoF1F1**\"&+c#F1FfalF1FcxF1FdoF1F7**\"%?^F1FfalF1)F\\pFSF1FdoF1F 1*(FjvF1F\\\\lF1F[pF1F7**\"%3YF1FduF1FbqF1FhqF1F7*,\"%?>F1FgrF1FbpF1Fj rF1F\\pF1F1*,\"%_6F1FgrF1FbpF1FjrF1FhqF1F7*,Fa[lF1FgrF1FbpF1FjrF1F\\uF 1F7*,Fa[lF1FgrF1FbpF1FjrF1FcxF1F1*,F\\[lF1FbqF1FbpF1FduF1F\\pF1F7*,F_b lF1FbqF1FbpF1FduF1FhqF1F1*,F\\[lF1FbqF1FbpF1FduF1F\\uF1F7*,F^[lF1FbqF1 FbpF1FduF1FcxF1F1**FialF1FfalF1F\\pF1FdoF1F1**FfwF1FjrF1FgrF1FhqF1F7*. Fg\\lF1F\\\\lF1F3F1F_pF1F\\pF1F`zF1F1*,F]]lF1FiuF1F[tF1FipF1F\\pF1F7*, Fi[lF1FiuF1F[tF1FipF1FhqF1F1*,\"%G)F1FddlF1FdoF1F\\uF1F1**FcdlF1Fddl F1FdoF1FhqF1F7F1*&FdoF1)FcpFgnF1F7#F;Fg`mFSF]oFgn+3F',$*&,(*&FeoF1F\\p F1F7FeoF1*&F\\pF1&%$sqrGFgoF1F1F1F3F7#!#;FS!\"#,$*&,0FcoF7**F;F1FdoF1F eoF1FbpF1F1*&F/F1F[pF1F7*(F/F1F\\pF1&FdcmF6F1F7*(F/F1F[pF1F\\pF1F1**F, F1FdoF1FeoF1F\\pF1F7**F,F1FdoF1F\\pF1FccmF1F1F1F3F7#F/FSF7,$*&,:*&FdoF 1F[pF1F/*&FbrF1FipF1F1**F;F1FbqF1FbpF1FeoF1F1**F_qF1F[pF1FdoF1FbpF1F7* (F,F1)F3F,F1FccmF1F1*(FSF1F3F1&%%spirGF6F1F7*(FbrF1F\\pF1&FdcmF+F1F1*( FbrF1FipF1F\\pF1F7**F;F1FbqF1FeoF1F\\pF1F7**F;F1FbqF1F\\pF1FccmF1F1**F GF1FdoF1F\\pF1F^dmF1F7**FGF1FdoF1F[pF1F\\pF1F1F1F3F7#F7FSFho,$*&,D*(Fg rF1F\\pF1FccmF1F_q**F_qF1FgrF1FeoF1F\\pF1F7**FdwF1FdoF1F\\pF1F`emF1F1* &FjtF1FjrF1F7*(\"$?\"F1F3F1&F^emF+F1F1*(\"#CF1F[emF1F^dmF1F7*(F\\qF1F[ emF1F4F1F1*(FjtF1F\\pF1&FdcmF:F1F7**F0F1FgrF1FeoF1FbpF1F1**FjvF1FbqF1F [pF1F\\pF1F1**FdwF1FdoF1FipF1F\\pF1F7**FjvF1FbqF1FbpF1F[pF1F7**F\\xF1F doF1FbpF1FipF1F1*(FguF1FdoF1FipF1F7*(F;F1FgrF1FeoF1F1**FjvF1FbqF1F\\pF 1F^dmF1F7*(FjtF1FjrF1F\\pF1F1F1F3F7#F1\"#gF1,$*&,RFhz!#v**Fi_lF1FiuF1F \\pF1F`emF1F1**FablF1FgrF1F\\pF1FefmF1F7**FhamF1FduF1FbqF1F\\pF1F7**Fi [mF1F\\\\lF1F\\pF1F^dmF1F7**FhamF1FbqF1F\\pF1&FdcmFKF1F1**F^glF1F3F1&F ^emF:F1FbqF1F7**\"$S'F1F[emF1F?F1FbqF1F7**FailF1F[emF1F`emF1FbqF1F1*(F hamF1FduF1FbqF1F1*(F_alF1FjrF1FgrF1F1*,F^[lF1Fb[lF1Fd[lF1Fe[lF1FeoF1F7 *(\"#IF1FizF1FeoF1F1**\"#XF1FizF1FeoF1FbpF1F1**Fi_lF1FiuF1FipF1F\\pF1F 7**F`yF1F\\\\lF1FbpF1F[pF1F7**Fi_lF1FiuF1FbpF1FipF1F1**Fi[mF1F\\\\lF1F [pF1F\\pF1F1*.F^[lF1Fb[lF1Fd[lF1Fe[lF1FeoF1F\\pF1F1**FablF1FjrF1FgrF1F \\pF1F1*(F`gmF1F\\\\lF1F[pF1F7*.F^[lF1Fb[lF1Fe[lF1Fd[lF1F\\pF1FccmF1F7 **F^dlF1FgrF1FbpF1FjrF1F7**\"#vF1FizF1F\\pF1FccmF1F1F1*&F3F1FbqF1F7#F7 \"%+7F,,$*&,\\o*,Fb[lF1Fe[lF1Fd[lF1F\\pF1F^dmF1!%sI*0F^[lF1FdoF1Fb[lF1 Fe[lF1Fd[lF1F\\pF1FccmF1F1**FcvF1F`glF1F\\pF1FccmF1F7**\"$+$F1FizF1F\\ pF1F^dmF1F1**\"%!o(F1FgrF1F\\pF1FjgmF1F7**FcamF1F\\pF1&FdcmFRF1FbqF1F1 **F_jmF1FduF1FgrF1F\\pF1F1**F_alF1F\\\\lF1F\\pF1F`emF1F7**F^dlF1FiuF1F \\pF1FefmF1F1**FablF1F[emF1FefmF1FbqF1F1**\"&+#>F1F3F1&F^emFKF1FbqF1F7 **FadlF1F[emF1FBF1FbqF1F7*(FcamF1FfalF1FbqF1F1**FcamF1FfalF1FbqF1F\\pF 1F7**F\\\\mF1FbpF1FduF1FgrF1F7*.Fa[lF1Fb[lF1Fd[lF1Fe[lF1FeoF1FdoF1F7*, F\\[lF1F[pF1Fb[lF1Fd[lF1Fe[lF1F7*(FhamF1FduF1FgrF1F1**F_alF1F\\\\lF1Fi pF1F\\pF1F1**F`yF1FizF1F[pF1FbpF1F1**FcvF1F`glF1FeoF1F\\pF1F1**F]jmF1F izF1F[pF1F\\pF1F7**FfhmF1F`glF1FbpF1FeoF1F7**Fi_lF1F\\\\lF1FipF1FbpF1F 7**F^dlF1FiuF1FjrF1F\\pF1F7**F^dlF1FiuF1FbpF1FjrF1F1*(Fi[mF1F\\\\lF1Fi pF1F7*(FfhmF1F`glF1FeoF1F7*0FcblF1FbpF1Fb[lF1Fd[lF1Fe[lF1FeoF1FdoF1F1* 0F^[lF1Fb[lF1Fd[lF1Fe[lF1FeoF1FdoF1F\\pF1F7*.F\\[lF1F[pF1Fb[lF1Fd[lF1F e[lF1F\\pF1F1*(F_fmF1FizF1F[pF1F1F1*&FbqF1F3F1F7#F7FbhlF;,$*&,jo*.FdoF 1FbpF1F[pF1Fb[lF1Fd[lF1Fe[lF1\"&cA$*(\"'?v5F1FfalF1FgrF1F1*.F_jmF1Fb[l F1Fe[lF1FeoF1Fd[lF1FbqF1F7*.F]^mF1Fb[lF1Fe[lF1F[pF1Fd[lF1FdoF1F7**\"'! 3I%F1FddlF1F\\pF1FbqF1F7**Fj\\nF1F\\pF1&FdcmFfnF1FbqF1F1**\"&gP&F1F[em F1FjgmF1FbqF1F1*,\"&;g)F1Fb[lF1Fe[lF1Fd[lF1FipF1F7**\"'+?nF1F3F1&F^emF RF1FbqF1F7**\"&!orF1F[emF1FWF1FbqF1F7*(FbjlF1F\\\\lF1FjrF1F7*(\"%g7F1F `glF1F[pF1F7*(\"%gLF1FizF1FipF1F1*(\"$?%F1)FdoFGF1FeoF1F1*0\"&/:#F1Fdo F1Fb[lF1Fe[lF1Fd[lF1F\\pF1F^dmF1F1*0F_^nF1FdoF1F[pF1Fb[lF1Fd[lF1Fe[lF1 F\\pF1F7*(Fj\\nF1FddlF1FbqF1F1*0F_jmF1FbqF1Fb[lF1Fd[lF1Fe[lF1FeoF1F\\p F1F1*0F_jmF1FbqF1Fb[lF1Fe[lF1Fd[lF1F\\pF1FccmF1F7*.F`]nF1Fb[lF1Fe[lF1F d[lF1F\\pF1F`emF1F7*.F`]nF1Fb[lF1Fe[lF1Fd[lF1FipF1F\\pF1F1**\"&S1)F1Fi uF1FbpF1FduF1F1**\"'gDKF1FgrF1FbpF1FfalF1F7**\"&!)o#F1F\\\\lF1F\\pF1Fe fmF1F7**F[_nF1F\\\\lF1FjrF1F\\pF1F1**\"%S]F1FizF1FbpF1FipF1F1**Fh]nF1F `glF1FbpF1F[pF1F7**\"'S]@F1FgrF1F\\pF1FajmF1F7**Fa_nF1FgrF1FfalF1F\\pF 1F1**\"&g,#F1F\\\\lF1FbpF1FjrF1F7**\"$N(F1F]^nF1F\\pF1FccmF1F1**Ff_nF1 F]^nF1FeoF1F\\pF1F7**\"%?DF1F`glF1F[pF1F\\pF1F1**Fi_nF1F`glF1F\\pF1F^d mF1F7**\"$:$F1F]^nF1FbpF1FeoF1F1**Fg^nF1FiuF1FduF1F\\pF1F7**Fg^nF1FiuF 1F\\pF1FjgmF1F1**\"%+%)F1FizF1FipF1F\\pF1F7**F``nF1FizF1F\\pF1F`emF1F1 F1*&FbqF1F3F1F7#F7\"'+W8F/F]oFS+1F',&F4F1*&#F;FaglF1*&F3F1&F^emFgoF1F1 F7Fho,&F?F,*&#F;FaglF1*&F3F1F]emF1F1F7F1,(FBF;*&#F;FaglF1*&F3F1F`fmF1F 1F7*(#FH\"#GF1F3F1F)F1F1F,,(*&F3F1F\\hmF1#FioFagl*(FdanF1F3F1F9F1F1*&F /F1FWF1F1F;,(*&F3F1FOF1Fdan*&FSF1F\\oF1F1*&#F;FaglF1*&F3F1FhjmF1F1F7F/ ,(*&F3F1FQF1Fdan*&FgnF1&F5FfnF1F1*&#F;FaglF1*&F3F1Fc]nF1F1F7FSF]oFgn+1 F',&Fi`n#F7F,F^dmF7Fho,&F^anFjbn*&F,F1F`emF1F7F1,&FbanFjbn*&F;F1FefmF1 F7F,,&FganFjbn*&F/F1FjgmF1F7F;,&F`bnFjbn*&FSF1FajmF1F7F/,&FgbnFjbn*&Fg nF1F\\]nF1F7FSF]oFgn" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%(dosolveGR6# %&ineqsG6$%$tmpG%%varsG6\"F+C+-%$mapG6%%'seriesG9$%\"tG-F.6$%)simplify G%\"%G-F.6%%(convertGF6%(polynomG-F.6%%'tcoeffGF6F2>8$F3>8%-%&minusG6$ -%'indetsG6#F?%*knownvarsG-%&printG6#FA-%&solveG6$F?FA-%'assignG6#F6F+ F+F+" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<(&%%spirG6#\"\"!&%%sG_3G6#\" \"\"&%$sqgGF*&%$sqrGF&&F/F*&%%sG_4G6#\"\"#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<(&%%spirG6#\"\"!&F%6#\"\"\"&%%sG_3G6#\"\"#&%$sqgGF-&%$ sqrGF-&%%sG_4G6#\"\"$" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<)&%%sG_5G6# \"\"$&%%spirG6#\"\"\"&F)6#\"\"#&%%sG_3GF&&%$sqgGF&&%$sqrGF&&%%sG_4G6# \"\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<)&%%sG_5G6#\"\"%&%$sqgGF&&%% spirG6#\"\"#&F+6#\"\"$&%%sG_3GF&&%$sqrGF&&%%sG_4G6#\"\"&" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#<*&%%sG_5G6#\"\"$&%$sqgG6#\"\"&&%%spirGF&&F-6#\" \"%&%$sqrGF*&%%sG_3GF*&%%sG_4G6#\"\"'F%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 156 "#Substitute in terms of sigma_0\nsigma(t):\n%;\nsimp lify(series(%,t,7)):\n%:\nsig_ser:=simplify(%);\nsqg[0]:=solve(sig0=co eff(%,t,0),sqg[0]);\n\nsimplify(sig_ser);\n" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#,$*&*(%&kappaG\"\"\")+1%\"tG*&*&)%#H0G\"\"#F'%$ommGF'F' %&omegaG!\"\"F',$*&F-F'F0F'#F'\"\"%F/,$*&*()F.\"\")F'%$omvGF')F0\"\"$F 'F'*$)F1F>F'F2#F'\"#5\"\"&,$*&*(F:F'FF>&Fbo 6#F6F6&Fbo6#FCFC&Fbo6#FKFKFVFQF'F'*&F&F')-%\"SG6#F*F6F'F2F>*&*,F>F'F^o F'FjnF')F.F6F'+3F*&%$sqrGFcoFdo&FjpFXF'&FjpFgoF/&FjpFioF>&FjpF[pF6&Fjp F]pFC&FjpF_pFKFVFQF'F'*(FHF'F&F'FapF'F2F'*&*,#FYF6F',&F'F'%#RcGF2F'F0F 'F-F'-%#qgGFdpF'F'*&F&F')FbpF>F'F2F'F'F'*$)%\"kGF/F'F2F/" }}{PARA 12 " " 1 "" {XPPMATH 20 "6#>%(sig_serG++%\"tG,$*&*&&%$sqgG6#\"\"!\"\"\",&! \"\"F.%#RvGF.F.F.,&F1\"\"%\"\"&F.F0!\"$F-,$*&,,*(%\"kGF.%#B0GF.)F1\"\" #F.F3**\"\"(F.%&omegaGF.F*F.F1F.F.*(F:F.F;F.F1F.F.*(F4F.F:F.F;F.F0*(F? F.F@F.F*F.F0F.,(*$F&%$sqgG6#\"\"!,$*&*&%%sig0G\"\"\",&%#RvG\"\"%\"\"&F,F,F,,&!\"\" F,F.F,F2#F2\"\"$" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#++%\"tG%%sig0G\"\" !,$*&,(*(%\"kG\"\"\"%#B0GF,%#RvGF,\"\"$*(F/F,F+F,F-F,!\"\"*(\"\"(F,%&o megaGF,F%F,F1F,,&F.\"\"%\"#:F,F1#F7\"#9F,,$*&,.*()F+\"\"#F,F%F,F.F,\"# ;**\"#GF,)F4F?F,F%F,F.F,F,*,\"$!=F,F4F,F+F,F-F,F.F,F,*(\"#gF,F>F,F%F,F ,*(\"$:$F,FCF,F%F,F1**FEF,F4F,F+F,F-F,F1F,*&F5F,,&F.F?F7F,F,F1#!#:\"$7 \"F?-%\"OG6#F,F/" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 104 "#This \+ is the solution for sigma\n#velocities\n#photon\nvg(t):=simplify(serie s(expandfully(3/4*qg(t)),t,3));" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>-% #vgG6#%\"tG++F',$*&*&%%sig0G\"\"\",&%#RvG\"\"%\"\"&F-F-F-,&!\"\"F-F/F- F3#F3F0\"\"!,$*&,0*(%&omegaGF-F,F-F/F-!\"%*(F1F-F:F-F,F-F3*,F0F-F:F-F, F-%#RcGF-F/F-F-**F1F-F:F-F,F-F>F-F-*(\"\"#F-%\"kGF-%#B0GF-F3**F0F-FBF- FCF-F/F-F-**FAF-FBF-FCF-)F/FAF-F3F-*$)F2FAF-F3#\"\"$\"#;F-,$*&*&F:F-,> **FBF-FCF-F/F-F>F-\"#7*,\"\"'F-FBF-FCF-FFF-F>F-F3*,F0F-F:F-F,F-F>F-FFF -F3*,\"#DF-F:F-F,F-F>F-F/F-F3**F0F-F:F-F,F-FFF-F-**\"#8F-F:F-F,F-F/F-F -**FQF-FBF-FCF-F/F-F3**FSF-FBF-FCF-FFF-F-*(\"#5F-F:F-F,F-F-*(FSF-FBF-F CF-F-**FSF-FBF-FCF-F>F-F3**FVF-F:F-F,F-F>F-F3*,FQF-F:F-F,F-)F>FAF-F/F- F-**\"#:F-F:F-F,F-F\\oF-F-F-F-*&,(F-F-*&FAF-F/F-F3*$FFF-F-F-F2F-F3#!\" $\"#kFA-%\"OG6#F-FJ" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 63 "#neu trino\n\nvr(t):=simplify(series(expandfully(3/4*qr(t)),t,3));" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>-%#vrG6#%\"tG++F',$*&*&%%sig0G\"\"\", &%#RvG\"\"%\"\"&F-F-F-F/!\"\"#F2F0\"\"!,$*&*(%\"kGF-%#B0GF-,&F2F-F/F-F -F-F/F2#!\"$\"\")F-,$*&*&)F8\"\"#F-F,F-F-F/F2#F-F=FB-%\"OG6#F-\"\"$" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 70 "#neutrino anisotropy stres s\nsimplify(series(expandfully(pir(t)),t,3));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#++%\"tG*&*&%#B0G\"\"\",&!\"\"F(%#RvGF(F(F(F+F*\"\"!,$*& *&%\"kGF(%%sig0GF(F(F+F*#!\"#\"\"$F(,$*&*&F0F(,(*(F0F(F'F(F+F(\"#X**\" #GF(%&omegaGF(F1F(F+F(F(*(F:F(F0F(F'F(F*F(F(*&,&F+\"\"%\"#:F(F(F+F(F*# F*\"#9\"\"#-%\"OG6#F(F4" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 117 "#Note that the above results have corrections from the full evolution , though the leading terms in B0 and sig0 are OK" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 " " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "10 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }