 against linear homogeneous perturbations in the background of $f(\mathcal{G},T)$ gravity. For EU, $a(t)=a_{0}=$ constant and consequently, the field equations (\[10\]) and (\[11\]) reduce to $$\begin{aligned}
\label{13}
\frac{3}{a_{0}^2}&=&\kappa^2\rho_{0}+\frac{1}{2}
f(\mathcal{G}_{0},T_{0})+(\rho_{0}+P_{0})f_{T}(\mathcal{G}_{0},T_{0}),
\\\label{14}-\frac{1}{a_{0}^2}&=&\kappa^2P_{0}-\frac{1}{2}f(\mathcal{G}_{0},T_{0}),\end{aligned}$$ where $\mathcal{G}_{0}=\mathcal{G}(a_{0})=0,~T_{0}=\rho_{0}-3P_{0},~\rho_{0}$ and $P_{0}$ are the unperturbed energy density and pressure, respectively. To explore the stability regions, we consider linear form of equation of state as $P(t)=\omega\rho(t)$ and define linear perturbations in the scale factor and energy density as follows $$\label{15}
a(t)=a_{0}+a_{0}\delta
a(t),\quad\rho(t)=\rho_{0}+\rho_{0}\delta\rho(t),$$ where $\delta a(t)$ and $\delta \rho(t)$ represent the perturbed scale factor and energy density, respectively. Applying the Taylor series expansion in two variables upto first order with the assumption that $f(\mathcal{G},T)$ is an analytic function, we have $$\label{16}
f(\mathcal{G},T)=f(\mathcal{G}_{0},T_{0})+
f_{\mathcal{G}}(\mathcal{G}_{0},T_{0})\delta\mathcal{G}
+f_{T}(\mathcal{G}_{0},T_{0})\delta T,$$ where $\delta\mathcal{G}$ and $\delta T$ have the following expressions $$\label{17}
\delta\mathcal{G}=\frac{24}{a_{0}^{2}}\delta\ddot{a},\quad \delta
T=T_{0}\delta\rho,$$ where $\delta\ddot{a}=\frac{d^2}{dt^2}(\delta a)$.
against linear homogeneous perturbations in the background of $ f(\mathcal{G},T)$ gravity. For EU, $ a(t)=a_{0}=$ changeless and consequently, the sphere equations (\[10\ ]) and (\[11\ ]) reduce to $ $ \begin{aligned } 
 \label{13 } 
 \frac{3}{a_{0}^2}&=&\kappa^2\rho_{0}+\frac{1}{2 } 
 f(\mathcal{G}_{0},T_{0})+(\rho_{0}+P_{0})f_{T}(\mathcal{G}_{0},T_{0 }), 
 \\\label{14}-\frac{1}{a_{0}^2}&=&\kappa^2P_{0}-\frac{1}{2}f(\mathcal{G}_{0},T_{0}),\end{aligned}$$ where $ \mathcal{G}_{0}=\mathcal{G}(a_{0})=0,~T_{0}=\rho_{0}-3P_{0},~\rho_{0}$ and $ P_{0}$ are the unperturbed department of energy density and pressure, respectively. To explore the constancy region, we consider analogue human body of equation of department of state as $ P(t)=\omega\rho(t)$ and define linear perturbation in the scale agent and energy density as succeed $ $ \label{15 } 
 a(t)=a_{0}+a_{0}\delta 
 a(t),\quad\rho(t)=\rho_{0}+\rho_{0}\delta\rho(t),$$ where $ \delta a(t)$ and $ \delta \rho(t)$ represent the perturbed plate factor and energy concentration, respectively. Applying the Taylor serial expansion in two variables upto first order with the assumption that $ f(\mathcal{G},T)$ is an analytic function, we have $ $ \label{16 } 
 f(\mathcal{G},T)=f(\mathcal{G}_{0},T_{0})+ 
 f_{\mathcal{G}}(\mathcal{G}_{0},T_{0})\delta\mathcal{G } 
 + f_{T}(\mathcal{G}_{0},T_{0})\delta T,$$ where $ \delta\mathcal{G}$ and $ \delta T$ own the following expressions $ $ \label{17 } 
 \delta\mathcal{G}=\frac{24}{a_{0}^{2}}\delta\ddot{a},\quad \delta 
 T = T_{0}\delta\rho,$$ where $ \delta\ddot{a}=\frac{d^2}{dt^2}(\delta a)$.
 agwinst linear homogeneous perturbations nb the uackgrohnd of $f(\oathcal{G},T)$ gravity. For EU, $a(t)=a_{0}=$ cinstabt and consequently, thd field eeuations (\[10\]) anv (\[11\]) reduce to $$\begmh{aligned}
\label{13}
\rvac{3}{a_{0}^2}&=&\kcp'a^2\rho_{0}+\frac{1}{2}
f(\mathcsl{G}_{0},T_{0})+(\rho_{0}+P_{0})f_{T}(\kathcal{G}_{0},T_{0}),
\\\label{14}-\xrxc{1}{c_{0}^2}&=&\kappa^2P_{0}-\frac{1}{2}f(\mathcal{G}_{0},T_{0}),\end{aligned}$$ wherq $\mathcsl{H}_{0}=\mathcal{G}(a_{0})=0,~T_{0}=\rho_{0}-3K_{0},~\rho_{0}$ sgd $P_{0}$ are the unperturbed energy densitg and pgessure, respectivrly. To explore the stabilihy rfgions, we consider linear forn of wquation of rtate as $P(u)=\okega\rho(t)$ ahd define linear perturbations kn thz scale facjod wtd energy dxnsity as follows $$\label{15}
a(t)=a_{0}+d_{0}\delta
a(y),\quad\rho(t)=\rho_{0}+\rhp_{0}\denta\eho(t),$$ where $\delta a(t)$ aid $\delta \rho(t)$ represgnt the pestbrbed scale factor ane wnergi denvity, eesoecuivxly. Applylng the Taylod series expansion in two varisbjvx upto first order wyth the assumption that $f(\mathcal{G},T)$ is at ahalytic function, we havw $$\label{16}
f(\mathcal{G},T)=f(\matjcal{G}_{0},T_{0})+
f_{\maehcal{G}}(\mathcal{G}_{0},T_{0})\delta\mathcal{G}
+f_{T}(\mathcal{G}_{0},T_{0})\delta T,$$ whare $\dxlga\mcbmcal{E}$ ajd $\delta T$ have the following expressions $$\labqm{17}
\dtlts\mathcal{G}=\frac{24}{a_{0}^{2}}\belta\ddot{a},\quad \drlha
Y=J_{0}\delta\rho,$$ wherg $\delta\bsof{a}=\frac{d^2}{dt^2}(\delta a)$.
against linear homogeneous perturbations in the background gravity. EU, $a(t)=a_{0}=$ and consequently, the reduce $$\begin{aligned} \label{13} \frac{3}{a_{0}^2}&=&\kappa^2\rho_{0}+\frac{1}{2} \\\label{14}-\frac{1}{a_{0}^2}&=&\kappa^2P_{0}-\frac{1}{2}f(\mathcal{G}_{0},T_{0}),\end{aligned}$$ where $\mathcal{G}_{0}=\mathcal{G}(a_{0})=0,~T_{0}=\rho_{0}-3P_{0},~\rho_{0}$ $P_{0}$ are the unperturbed energy density pressure, respectively. To explore the stability regions, we consider linear form of equation state as $P(t)=\omega\rho(t)$ and define linear perturbations in the scale factor and energy as $$\label{15} a(t),\quad\rho(t)=\rho_{0}+\rho_{0}\delta\rho(t),$$ $\delta a(t)$ and $\delta \rho(t)$ represent the perturbed scale factor and energy density, respectively. Applying the series expansion in two variables upto first order the assumption that $f(\mathcal{G},T)$ an analytic function, we have f(\mathcal{G},T)=f(\mathcal{G}_{0},T_{0})+ +f_{T}(\mathcal{G}_{0},T_{0})\delta T,$$ $\delta\mathcal{G}$ $\delta have the following $$\label{17} \delta\mathcal{G}=\frac{24}{a_{0}^{2}}\delta\ddot{a},\quad \delta T=T_{0}\delta\rho,$$ where $\delta\ddot{a}=\frac{d^2}{dt^2}(\delta a)$.
 against linear homogeneous pErturbatioNs in tHe bAckGrOund Of $f(\mAthcal{G},T)$ gravitY. for Eu, $a(t)=a_{0}=$ constant and consequEntly, ThE FielD EqUatioNs (\[10\]) and (\[11\]) reDUcE TO $$\beGiN{aLigNeD}
\LaBel{13}
\frAc{3}{a_{0}^2}&=&\Kappa^2\rhO_{0}+\frac{1}{2}
f(\mathCal{g}_{0},T_{0})+(\Rho_{0}+P_{0})f_{T}(\mathcaL{g}_{0},T_{0}),
\\\Label{14}-\frac{1}{a_{0}^2}&=&\KapPa^2P_{0}-\frac{1}{2}f(\mathCal{g}_{0},T_{0}),\end{aLiGneD}$$ Where $\MatHcal{G}_{0}=\MathcaL{g}(a_{0})=0,~T_{0}=\rho_{0}-3p_{0},~\rho_{0}$ and $P_{0}$ aRe THe unpeRTurbed eNERgY denSity and pressure, reSPeCTively. To explorE the stAbILiTY RegIonS, we consideR lInear FOrm of eqUAtION Of sTAte as $P(t)=\omega\rHo(t)$ and definE LinEar perTuRbaTIons in The scAlE FacTor and energY denSity as folLows $$\laBEl{15}
a(t)=a_{0}+a_{0}\dELta
a(t),\quAd\rho(t)=\Rho_{0}+\Rho_{0}\DeltA\RhO(t),$$ WheRe $\DEltA A(t)$ And $\DEltA \rho(t)$ repReSeNt the PertURBED scaLe fActoR and eNergy density, rEspEctiVEly. applyIng thE TayLoR seriEs expaNsion In Two variables uptO firSt order wiTh tHe AssUmPtion THat $f(\maThcAl{G},t)$ is an anAlytic fUNctIoN, WE HaVe $$\label{16}
f(\mathcal{G},T)=f(\MaTHCaL{G}_{0},T_{0})+
f_{\mathCal{G}}(\maTHcAl{g}_{0},t_{0})\delta\maThCal{g}
+f_{T}(\mATHcal{G}_{0},t_{0})\delTA T,$$ Where $\delTa\mathCAl{g}$ aNd $\delta t$ hAve the FoLloWinG exprESsioNs $$\labeL{17}
\delta\maThcal{g}=\Frac{24}{a_{0}^{2}}\delta\ddot{A},\Quad \delta
T=T_{0}\deLTa\RHO,$$ wHEre $\dEltA\ddot{a}=\frac{d^2}{Dt^2}(\deLTa a)$.
 against linear homogeneou s perturba tions  in  th eback grou nd of $f(\math c al{G },T)$ gravity. For EU,  $a(t )= a _{0} = $const ant and co n s equ en tl y,th e  f ieldequ ations(\[10\]) a nd(\ [11\]) reduc e  t o $$\begin {al igned}
\labe l{1 3}
\fr ac {3} { a_{0} ^2} &=&\k appa^2 \ rho_{0 }+\frac{1 }{ 2 }
f(\m a thcal{G } _ {0 },T_ {0})+(\rho_{0}+P_ { 0} ) f_{T}(\mathcal {G}_{0 }, T _{ 0 } ),\\\ label{14}- \f rac{1 } {a_{0}^ 2 }& = & \ kap p a^2P_{0}-\fra c{1}{2}f(\m a thc al{G}_ {0 },T _ {0}),\ end{a li g ned }$$ where $ \mat hcal{G}_{ 0}=\ma t hcal{G} ( a_{0})= 0,~T_{ 0}= \rh o_{0 } -3 P_ {0} ,~ \ rho _ {0 }$a nd$P_{0}$ar ethe u nper t u r b ed e ner gy d ensit y and pressur e,resp e cti vely.  To e xplo re  thestabil ity r eg ions, we consid er l inear for m o fequ at ion o f  state  as  $P (t)=\om ega\rho ( t)$  a n d de fine linear pertur ba t i on s in the  scale fa ct o r and en er gydens i t y asfoll o ws  $$\labe l{15}a (t )= a_{0}+a _{ 0}\del ta 
a( t), \quad \ rho( t)=\rh o_{0}+\r ho_{0 } \delta\rho(t), $ $ where $\del t aa ( t) $  and  $\ delta \rho( t)$r epre sent th e p e rturb ed sc al e  f a ctor and energy den si ty, re spect ively. Applyi ng the Tay l o r  seriesexpa n si o n in two varia blesupto first order wi th th e assump tion that $ f(\mathc al{ G}, T)$  is a nanalytic func t i on,we  have $ $\l abel{16 }
f (\m ath cal {G },T)=f(\m athcal{G }_ {0 }, T_ {0} )+
f_ { \mathcal {G }}( \m ath cal{G } _{0},T _{0}) \del ta \m a thc al{G}
+ f _{ T } (\ma th ca l{G} _{0 }, T_{0} )\de l taT,$$ wh ere $\del ta\ m athc al {G }$ and$\delta T$ ha ve  the follo wi ngexpres s i ons $$\l abel{17}
\delta\mathcal { G}=\fra c{2 4}{a_ {0}^ {2}}\delt a\d dot{a} ,\q u ad \de lta
T= T_{0} \d elt a \ rho,$ $ wh ere  $ \delta\ddo t { a}= \frac {d ^2}{ dt^2}(\ delta a)$.
 against_linear homogeneous_perturbations in the background_of $f(\mathcal{G},T)$_gravity._For EU,_$a(t)=a_{0}=$_constant and consequently,_the field equations_(\[10\]) and (\[11\]) reduce_to $$\begin{aligned}
\label{13}
\frac{3}{a_{0}^2}&=&\kappa^2\rho_{0}+\frac{1}{2}
f(\mathcal{G}_{0},T_{0})+(\rho_{0}+P_{0})f_{T}(\mathcal{G}_{0},T_{0}),
\\\label{14}-\frac{1}{a_{0}^2}&=&\kappa^2P_{0}-\frac{1}{2}f(\mathcal{G}_{0},T_{0}),\end{aligned}$$ where_$\mathcal{G}_{0}=\mathcal{G}(a_{0})=0,~T_{0}=\rho_{0}-3P_{0},~\rho_{0}$_and $P_{0}$ are the unperturbed energy density and pressure, respectively. To explore the stability_regions,_we consider_linear_form_of equation of state as_$P(t)=\omega\rho(t)$ and define linear perturbations_in the_scale factor and energy density as follows $$\label{15}
a(t)=a_{0}+a_{0}\delta
a(t),\quad\rho(t)=\rho_{0}+\rho_{0}\delta\rho(t),$$_where_$\delta a(t)$ and_$\delta \rho(t)$ represent the perturbed scale factor and energy_density, respectively. Applying the Taylor series_expansion in two_variables_upto_first order with the_assumption that $f(\mathcal{G},T)$ is an analytic_function, we have $$\label{16}
f(\mathcal{G},T)=f(\mathcal{G}_{0},T_{0})+
f_{\mathcal{G}}(\mathcal{G}_{0},T_{0})\delta\mathcal{G}
+f_{T}(\mathcal{G}_{0},T_{0})\delta T,$$ where_$\delta\mathcal{G}$ and $\delta T$ have the following_expressions $$\label{17}
\delta\mathcal{G}=\frac{24}{a_{0}^{2}}\delta\ddot{a},\quad \delta
T=T_{0}\delta\rho,$$ where $\delta\ddot{a}=\frac{d^2}{dt^2}(\delta a)$.