bar{\Omega}   \left(3 (\gamma  \bar{\Omega}_{m}+\Omega -1)+2 \bar{\Sigma}^2\right), \\
& \partial_{\tau}{\bar{\Omega}}_m= \bar{\Omega}_{m}  \left(3 \gamma  (\bar{\Omega}_{m}-1)+2 \bar{\Sigma}^2+3 \bar{\Omega} \right), \\
& \partial_{\tau}{\bar{\Sigma}}= \frac{1}{2} \bar{\Sigma}  \left(3 (\gamma  \bar{\Omega}_{m}+\bar{\Omega} -2)+2 \bar{\Sigma} ^2\right).\end{aligned}$$

Analyzing the equilibrium points of the system, which is the same as, we obtain points of the form $P_1 = (0,1,0)$, $P_2 = (1,0,0)$, $ P_3 =(0,0, - \sqrt{3})$, $ P_4 = (0,0, \sqrt{3}) $ and $P_5 = (0,0,0)$. By evaluating the points in the linearization of and calculating their respective eigenvalues, we obtain the results shown in the table \[PCsis3.31\]. This means that on average the dynamics of the models and is the same which indicates that the asymptotic behavior of the model on average is independent of the coupling function, of the geometry of model (FLRW, Bianchi I). Although, obviously, non-averaged systems have different dynamics.

Averaging methods for the scalar field with generalized harmonic potential $
V(\phi)= \frac{\phi ^2}{2}+f\left[1- \cos \left(\frac{\phi }{f}\right)\right]
$, $f> 0$ nonminimally coupled to matter with coupling $\chi(\phi)=\chi_0 e^{\frac{\lambda \phi}{4-3\gamma}}$. {#SECT:4.5}
==============================================================================================================

In this section the averaging methods are applied for FLRW metrics and Bianchi I metrics for the generalized harmonic potential $
V(\phi)= \frac{\phi ^2}{2}+f\left[1- \cos \left(\frac{\phi }{f}\right)\right]
$, $f> 0$, coupled to matter with coupling function $\chi=\chi_0 e^{\frac
bar{\Omega }    \left(3 (\gamma   \bar{\Omega}_{m}+\Omega -1)+2 \bar{\Sigma}^2\right), \\ 
 & \partial_{\tau}{\bar{\Omega}}_m= \bar{\Omega}_{m }   \left(3 \gamma   (\bar{\Omega}_{m}-1)+2 \bar{\Sigma}^2 + 3 \bar{\Omega } \right), \\ 
 & \partial_{\tau}{\bar{\Sigma}}= \frac{1}{2 } \bar{\Sigma }   \left(3 (\gamma   \bar{\Omega}_{m}+\bar{\Omega } -2)+2 \bar{\Sigma } ^2\right).\end{aligned}$$ 

 Analyzing the equilibrium points of the system, which is the like as, we receive points of the form $ P_1 = (0,1,0)$, $ P_2 = (1,0,0)$, $ P_3 =( 0,0, - \sqrt{3})$, $ P_4 = (0,0, \sqrt{3 }) $ and $ P_5 = (0,0,0)$. By measure the point in the linearization of and calculating their respective eigenvalues, we prevail the resultant role shown in the mesa \[PCsis3.31\ ]. This means that on average the moral force of the models and is the same which indicates that the asymptotic behavior of the exemplar on average is independent of the coupling function, of the geometry of model (FLRW, Bianchi I). Although, obviously, non - average systems have different moral force. 

 Averaging methods for the scalar plain with generalized harmonic electric potential $ 
 V(\phi)= \frac{\phi ^2}{2}+f\left[1- \cos \left(\frac{\phi } { f}\right)\right ] 
 $, $ f > 0 $ nonminimally coupled to count with coupling $ \chi(\phi)=\chi_0 e^{\frac{\lambda \phi}{4 - 3\gamma}}$. { # SECT:4.5 } 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 In this section the averaging methods are applied for FLRW metric function and Bianchi I metrics for the generalized harmonic potential $ 
 V(\phi)= \frac{\phi ^2}{2}+f\left[1- \cos \left(\frac{\phi } { f}\right)\right ] 
 $, $ f > 0 $, coupled to matter with yoke affair $ \chi=\chi_0 e^{\frac
bar{\Lmega}   \left(3 (\gamma  \bar{\Ooega}_{m}+\Omega -1)+2 \bar{\Sigma}^2\rijht), \\
& \padtial_{\tau}{\car{\Omega}}_m= \bar{\Omega}_{m}  \left(3 \galmq  (\bae{\Omega}_{m}-1)+2 \bar{\Sigma}^2+3 \bar{\Omdga} \right), \\
& \partiao_{\tau}{\uar{\Sigma}}= \frac{1}{2} \bac{\Aigma}  \left(3 (\gzlma  \uar{\Omega}_{m}+\bar{\Omena} -2)+2 \bar{\Sigmd} ^2\right).\end{aligtea}$$

Aualyzing the equilibrium points of tre systrm, which is the fame ws, ws obtain points of the form $P_1 = (0,1,0)$, $P_2 = (1,0,0)$, $ P_3 =(0,0, - \vqrt{3})$, $ P_4 = (0,0, \sqry{3}) $ and $P_5 = (0,0,0)$. By evaluating tje plints in the lineagization of and xalculating gheir resptccive eigenvzlues, we obtain the results shodn in the table \[PXsid3.31\]. This means that on average bne dyndmics og the models akd is thw same which indicatev that the asymptojic behavimr of the model on qvwrage is hndeownddnt oh tge couoliig function, of the geonetry of model (FLRW, Byqnchi I). Althohgh, obdijusly, non-averaged systems have differenu dynzmics.

Averaging methods dor the scalar field aith genewalized harmonic potential $
V(\phi)= \frac{\phi ^2}{2}+f\left[1- \cos \left(\hrxc{\pkl }{w}\eihht)\right]
$, $f> 0$ nonminimally coupled to matter wieg volpling $\chi(\phi)=\chi_0 c^{\frac{\lambda \phi}{4-3\gamka}}$. {#SRST:4.5}
==============================================================================================================

In this secjion thz aberaging methods age applyed fir FLRW mttricx and Bianchi I metrics for the generalpzed harmonic potentiap $
V(\phi)= \frac{\'hi ^2}{2}+f\lgft[1- \cox \left(\frac{\phi }{f}\right)\rigkt]
$, $f> 0$, doupled to latter wiff coupling functkon $\cvi=\chi_0 e^{\frac
bar{\Omega} \left(3 (\gamma \bar{\Omega}_{m}+\Omega -1)+2 \bar{\Sigma}^2\right), \\ \bar{\Omega}_{m} \gamma (\bar{\Omega}_{m}-1)+2 \bar{\Omega} \right), \\ (\gamma -2)+2 \bar{\Sigma} ^2\right).\end{aligned}$$ the equilibrium points the system, which is the same we obtain points of the form $P_1 = (0,1,0)$, $P_2 = (1,0,0)$, $ =(0,0, - \sqrt{3})$, $ P_4 = (0,0, \sqrt{3}) $ and $P_5 = (0,0,0)$. evaluating points the of and calculating their respective eigenvalues, we obtain the results shown in the table \[PCsis3.31\]. This that on average the dynamics of the models is the same which that the asymptotic behavior of model average is of coupling of the geometry model (FLRW, Bianchi I). Although, obviously, non-averaged systems have different dynamics. Averaging methods for the scalar field generalized harmonic V(\phi)= \frac{\phi \cos }{f}\right)\right] $f> 0$ nonminimally matter with coupling $\chi(\phi)=\chi_0 e^{\frac{\lambda \phi}{4-3\gamma}}$. this section the averaging methods are applied for metrics and I metrics for the generalized harmonic $ V(\phi)= \frac{\phi ^2}{2}+f\left[1- \cos \left(\frac{\phi }{f}\right)\right] $, 0$, coupled to matter with coupling function $\chi=\chi_0 e^{\frac
bar{\Omega}   \left(3 (\gamma  \bar{\Omega}_{M}+\Omega -1)+2 \bar{\SIgma}^2\rIghT), \\
& \paRtIal_{\tAu}{\baR{\Omega}}_m= \bar{\OmegA}_{M}  \lefT(3 \gamma  (\bar{\Omega}_{m}-1)+2 \bar{\SigmA}^2+3 \bar{\OMeGA} \rigHT), \\
& \pArtiaL_{\tau}{\bar{\sIgMA}}= \FraC{1}{2} \bAr{\sigMa}  \LEfT(3 (\gammA  \baR{\Omega}_{m}+\Bar{\Omega} -2)+2 \baR{\SiGmA} ^2\right).\end{aliGNeD}$$

Analyzing The Equilibrium pOinTs of thE sYstEM, whicH is The saMe as, we OBtain pOints of thE fORm $P_1 = (0,1,0)$, $P_2 = (1,0,0)$, $ P_3 =(0,0, - \sQRt{3})$, $ P_4 = (0,0, \sqrt{3}) $ AND $P_5 = (0,0,0)$. by evAluating the points IN tHE linearization Of and cAlCUlATIng TheIr respectiVe EigenVAlues, we OBtAIN The REsults shown in The table \[PCsIS3.31\]. ThIs meanS tHat ON averaGe the DyNAmiCs of the modeLs anD is the samE which INdicateS That the AsymptOtiC beHaviOR oF tHe mOdEL on AVeRagE Is iNdependeNt Of The coUpliNG FUNctiOn, oF the GeomeTry of model (FLRw, BiAnchI i). AlThougH, obviOuslY, nOn-aveRaged sYstemS hAve different dynAmicS.

AveraginG meThOds FoR the sCAlar fiEld WitH generaLized haRMonIc POTEnTial $
V(\phi)= \frac{\phi ^2}{2}+f\leFt[1- \COS \lEft(\frac{\pHi }{f}\rigHT)\rIgHT]
$, $f> 0$ nonminImAllY couPLEd to mAtteR WiTh coupliNg $\chi(\pHI)=\cHi_0 E^{\frac{\laMbDa \phi}{4-3\gAmMa}}$. {#SeCT:4.5}
==============================================================================================================

in thiS SectIon the AveraginG methODs are applied foR fLRW metrics anD biANChI i metRicS for the geneRaliZEd haRmonIC pOteNTial $
V(\Phi)= \frAc{\PHi ^2}{2}+F\Left[1- \cos \left(\frac{\phi }{f}\RiGht)\rigHt]
$, $f> 0$, coUpled to matter With coupliNG FUnction $\cHi=\chI_0 E^{\fRAc
bar{\Omega}   \left(3 (\ga mma  \bar{ \Omeg a}_ {m} +\ Omeg a -1 )+2 \bar{\Sigm a }^2\ right), \\
& \partial_ {\tau }{ \ bar{ \ Om ega}} _m= \ba r {\ O m ega }_ {m } \l e ft (3 \g amm a  (\ba r{\Omega}_ {m} -1 )+2 \bar{\Si g ma }^2+3 \bar {\O mega} \right ),\\
& \ pa rti a l_{\t au} {\bar {\Sigm a }}= \f rac{1}{2}  \ b ar{\Si g ma}  \l e f t( 3 (\ gamma  \bar{\Omeg a }_ { m}+\bar{\Omega } -2)+ 2\ ba r { \Si gma } ^2\right ). \end{ a ligned} $ $A n aly z ing the equil ibrium poin t s o f thesy ste m , whic h isth e  sa me as, we o btai n pointsof the form $P _ 1 = (0, 1,0)$,  $P _2= (1 , 0, 0) $,$P _3= (0 ,0, - \ sqrt{3}) $,  $  P_4= (0 , 0 , \sqr t{3 }) $  and$P_5 = (0,0,0 )$.  Bye val uatin g the  poi nt s inthe li neari za tion of and cal cula ting thei r r es pec ti ve ei g envalu es,  we  obtain  the re s ult ss h o wn  in the table \[PC si s 3 .3 1\]. Thi s mean s  t ha t  on aver ag e t he d y n amics  oft he  modelsand is th esame wh ic h indi ca tes  th at th e  asy mptoti c behavi or of the model on a v erage is inde p en d e nt of t hecoupling fu ncti o n, o f th e  g eom e try o f mod el (F L RW, Bianchi I). Alt ho ugh, o bviou sly, non-aver aged syste m s have dif fere n td ynamics.

Aver aging  methods f o r the sc alarfield wi th genera l i zed harm oni c p ote nti a l  $ 
V(\phi)= \fr a c {\ph i^2}{2}+ f\l eft[1-\co s \ lef t(\ fr ac{\phi } {f}\righ t) \r ig ht ]
$ , $f> 0$ nonmi ni mal ly  co upled to mat ter w ithco up l ing  $\chi( \ ph i ) =\ch i_ 0e^{\ fra c{ \lamb da \ p hi} {4-3\ga mma}}$. { #SE C T:4. 5} 
= ======= ============= == ========== == === ====== = = ======== ======================= = ======= === ===== ==== ========= === =

Inthi s  secti on the  aver ag ing m ethod s ar e a pp lied for F L R W m etric sandBianchi  I metrics for the gen eralized harm oni c po t e nt ial $V (\p hi ) = \ f r ac{\phi ^2}{2}+ f\left[1-\c o s\left(\fra c {\p hi  }{f}\r ight)\r ight] $, $f>0$, coupl ed to mat te r wi t h  co upling fun ction $\ chi=\chi_ 0  e^{\ f ra c
bar{\Omega} _ \left(3_(\gamma  \bar{\Omega}_{m}+\Omega -1)+2_\bar{\Sigma}^2\right), \\
&_\partial_{\tau}{\bar{\Omega}}_m=_\bar{\Omega}_{m} _\left(3_\gamma  (\bar{\Omega}_{m}-1)+2_\bar{\Sigma}^2+3 \bar{\Omega} \right),_\\
& \partial_{\tau}{\bar{\Sigma}}= \frac{1}{2} \bar{\Sigma}_ \left(3 (\gamma__\bar{\Omega}_{m}+\bar{\Omega} -2)+2 \bar{\Sigma} ^2\right).\end{aligned}$$

Analyzing the equilibrium points of the system, which is the same_as,_we obtain_points_of_the form $P_1 = (0,1,0)$,_$P_2 = (1,0,0)$, $ P_3_=(0,0, -_\sqrt{3})$, $ P_4 = (0,0, \sqrt{3}) $ and_$P_5_= (0,0,0)$. By_evaluating the points in the linearization of and calculating_their respective eigenvalues, we obtain the_results shown in_the_table_\[PCsis3.31\]. This means that_on average the dynamics of the_models and is the same which_indicates that the asymptotic behavior of the_model on average is independent of_the coupling function, of the_geometry of_model (FLRW, Bianchi I). Although,_obviously, non-averaged systems_have different_dynamics.

Averaging methods for_the scalar field with generalized harmonic_potential $
V(\phi)= \frac{\phi_^2}{2}+f\left[1- \cos \left(\frac{\phi }{f}\right)\right]
$, $f> 0$_nonminimally_coupled to matter_with_coupling_$\chi(\phi)=\chi_0 e^{\frac{\lambda_\phi}{4-3\gamma}}$. {#SECT:4.5}
==============================================================================================================

In this_section_the averaging_methods_are applied for FLRW metrics and_Bianchi_I metrics for the generalized harmonic potential_$
V(\phi)= \frac{\phi ^2}{2}+f\left[1- \cos_\left(\frac{\phi_}{f}\right)\right]
$, $f> 0$, coupled_to matter with coupling function_$\chi=\chi_0 e^{\frac