 distance keeps on attaining linearity. At infinite time, the steady state values are linear with position as given by Eq. \[infinite-time\]. Similarly, the material energy density initially exhibits slight non-linear variation and finally attains the linearity (Fig. \[vxt-planar\]). The non-linear variation at early stages occurs due to net absorption of energy by the initially cold material (as u(x,0)=v(x,0)=0). Initially, the material energy density is found to lag behind the radiation energy densities and finally equilibrate as time proceeds (beyond $\tau$=10). In this work, all the results have been obtained by considering contribution from the first 30 roots of the transcendental equation. The value of opacity $\kappa$ is chosen to be 100 and $\epsilon $ equals 0.1. For a heat wave traveling into a thin plate and composite planar slab, a similar linear variation in temperature with distance was observed though difference existed in the space and time dependent behaviour due to heat conduction approximation [@Sarkar],[@Sun].

![\[uxt-planar\]Scaled radiation energy density $u(x,\tau)$ vs position (x) in the slab of scaled thickness $b=1$ at different times for $ \epsilon = 0.1 $. The symbols stand for analytical values whereas lines represent the results obtained from finite difference method.](figure4.eps){width="8cm"}

![\[vxt-planar\]Scaled material energy density $v(x,\tau)$ vs position (x) in the slab at different times for $\epsilon = 0.1$. The symbols stand for analytical values whereas lines represent the results obtained from finite difference method.](figure5.eps){width="8cm"}

The first derivatives w.r.t. position of the analytical radiation and material energy density are plotted in Figs. \[uprimext-planar\] and \[vprimext-planar\]. As the radiation and material energy densities decrease with x, the derivative has negative values. The derivative has a greater negative value at the left compared to the right zone. As both radiation and material energy densities keep on increasing with time due to radiation diffusion, magnitude of the gradient decreases for the left and increases for the right sides. The gradient of both radiation and material energy densities obtain a constant value of $\frac{-3}{3+4\sqrt{3}}=-0.30217$ after infinite time showing that there is a constant leakage of
distance keeps on attaining linearity. At infinite meter, the regular state values are linear with status as given by Eq. \[infinite - time\ ]. Similarly, the material department of energy density initially exhibits slight non - analogue mutant and finally reach the linearity (Fig. \[vxt - planar\ ]). The non - linear magnetic declination at early stage occurs due to net absorption of energy by the initially cold material (as u(x,0)=v(x,0)=0). Initially, the fabric energy density is found to imprison behind the radiation energy densities and finally equilibrate as clock time proceeds (beyond $ \tau$=10). In this work, all the results have been obtained by considering contribution from the inaugural 30 roots of the transcendental equation. The value of opacity $ \kappa$ is chosen to be 100 and $ \epsilon $ equals 0.1. For a heat wave traveling into a thin plate and composite planar slab, a similar linear variation in temperature with distance was observed though difference existed in the space and time dependent behavior due to heat conduction approximation [ @Sarkar],[@Sun ]. 

! [ \[uxt - planar\]Scaled radiotherapy department of energy concentration $ u(x,\tau)$ vs position (x) in the slab of scaled thickness $ b=1 $ at different times for $ \epsilon = 0.1 $. The symbols stand for analytical values whereas line represent the results obtained from finite difference method.](figure4.eps){width="8 cm " } 

! [ \[vxt - planar\]Scaled material energy density $ v(x,\tau)$ vs place (x) in the slab at different times for $ \epsilon = 0.1$. The symbols stand for analytic values whereas pipeline represent the results obtained from finite difference method.](figure5.eps){width="8 cm " } 

 The inaugural derivatives w.r.t. position of the analytical radiation sickness and material energy density are plotted in Figs. \[uprimext - planar\ ] and \[vprimext - planar\ ]. As the radiotherapy and material department of energy densities decrease with x, the derivative has negative values. The derived function has a capital negative value at the left field compared to the right zone. As both radiation and fabric energy densities keep on increasing with time due to radiation diffusion, magnitude of the gradient decreases for the left and increases for the right side. The gradient of both radiation and material department of energy densities obtain a constant value of $ \frac{-3}{3 + 4\sqrt{3}}=-0.30217 $ after infinite time showing that there is a ceaseless leakage of
 didtance keeps on attaininn linearity. At iudinite time, fhe steaay state values are linear wmth posituon as given by Eq. \[infknite-time\]. Similaroy, tie material enerjg densibv iniflally xxhibits slight non-linear variation and fknclly attains the linearity (Fig. \[vxt-plwnar\]). Thr jon-linear variwtiom at sarly stages occurs due to net abskrption of energy by the initially cold materiwl (ad u(x,0)=v(x,0)=0). Initially, tje material enewty density ir found to lag behind jhe radiation energy densities avd fiually equilubeatf as time priceedf (beyond $\tau$=10). In this work, akl the results hate bwen obtained by consivering contribution srom the xixst 30 roots of the trabsxendettal equxriov. Tge vzlue ov o'acity $\kappz$ is chosen to be 100 and $\epsilon $ qwuals 0.1. For a geat wwvq traveling into a thin plate and compovits planar slab, a similar linear variation in jemperaturq with distance was observed though difference exhsted kn uhc rpafe and time dependent behaviour due to heat cjhdicnion approximatiok [@Sarkar],[@Sun].

![\[uxt-plansr\]Dcsjed radiation energv dsnsity $u(x,\tau)$ vs podition (v) in rhe slab jf svaled thickness $b=1$ at differwnt times fog $ \epsilon = 0.1 $. The symblls stand fur amalytocal values whereas linzs repdesent the gesults oggained from finige cixference method.](figure4.eps){wiqth="8cm"}

![\[vxt-poanax\]Scaled oatetial enqrgy densihy $v(x,\bdu)$ vs position (x) ij the snab at difverent times for $\epsilon = 0.1$. The symbols stand gos atalyticaj valmes whereas linqs represent tke resulcs obtxined from finite difference method.](figure5.aks){width="8cm"}

The hirst deryvatuves w.r.t. porktion of the amalytical radiation and material enerny devaity are plotteb un Figs. \[uprimexy-plxnaw\] wnv \[vprykext-planar\]. Av thd rxciatiun and mateviau enrrgy densities decredse sith x, the derivatovc has negqtive vajues. The derifative has a greattr negetive talue st jhe left compared to the right aone. As blth radiation agd mqterial eneryy densities keep on increasing with timx due to radiation diffysion, magnitude of jhe gradient dxcreasqs for tha left and increases for the right sices. The gradient of bkth ragiatiln and material energy densities obtain a constant value of $\frac{-3}{3+4\sqrt{3}}=-0.30217$ aftee infiiiee time shosing that tkerz is a sonsvait leakage of
distance keeps on attaining linearity. At infinite steady values are with position as the energy density initially slight non-linear variation finally attains the linearity (Fig. \[vxt-planar\]). non-linear variation at early stages occurs due to net absorption of energy by initially cold material (as u(x,0)=v(x,0)=0). Initially, the material energy density is found to behind radiation densities finally equilibrate as time proceeds (beyond $\tau$=10). In this work, all the results have been obtained considering contribution from the first 30 roots of transcendental equation. The value opacity $\kappa$ is chosen to 100 $\epsilon $ 0.1. a wave traveling into thin plate and composite planar slab, a similar linear variation in temperature with distance was observed though existed in and time behaviour to conduction approximation [@Sarkar],[@Sun]. energy density $u(x,\tau)$ vs position (x) of scaled thickness $b=1$ at different times for \epsilon = $. The symbols stand for analytical whereas lines represent the results obtained from finite method.](figure4.eps){width="8cm"} ![\[vxt-planar\]Scaled material energy density $v(x,\tau)$ vs position (x) in the slab at different times = 0.1$. The symbols for analytical values lines the obtained finite difference The first derivatives w.r.t. position of the analytical radiation and material density are plotted in Figs. \[uprimext-planar\] and \[vprimext-planar\]. As the material densities decrease with the derivative has negative The has a greater negative the compared zone. both and material energy densities on increasing with time due radiation diffusion, magnitude of left and increases for the right sides. The of both radiation and material energy densities a constant value of $\frac{-3}{3+4\sqrt{3}}=-0.30217$ after infinite time showing that there is constant leakage
 distance keeps on attaining lInearity. At InfinIte TimE, tHe stEady State values are LIneaR with position as given by eq. \[infInITe-tiME\]. SImilaRly, the mATeRIAl eNeRgY deNsITy InitiAllY exhibiTs slight noN-liNeAr variation aND fInally attaIns The linearity (fig. \[Vxt-plaNaR\]). ThE Non-liNeaR variAtion aT Early sTages occuRs DUe to neT AbsorptION oF eneRgy by the initially COlD Material (as u(x,0)=v(x,0)=0). initiaLlY, ThE MAteRiaL energy denSiTy is fOUnd to laG BeHIND thE Radiation enerGy densities ANd fInally EqUilIBrate aS time PrOCeeDs (beyond $\tau$=10). in thIs work, all The resULts have BEen obtaIned by ConSidErinG CoNtRibUtIOn fROm The FIrsT 30 roots of ThE tRanscEndeNTAL EquaTioN. The Value Of opacity $\kappA$ is ChosEN to Be 100 and $\EpsilOn $ eqUaLs 0.1. For A heat wAve trAvEling into a thin pLate And composIte PlAnaR sLab, a sIMilar lIneAr vAriatioN in tempERatUrE WITh Distance was observeD tHOUgH differeNce exiSTeD iN The space AnD tiMe dePENdent BehaVIoUr due to hEat conDUcTiOn approXiMation [@saRkaR],[@SuN].

![\[uxt-pLAnar\]scaled RadiatioN enerGY density $u(x,\tau)$ vS Position (x) in thE SlAB Of SCaleD thIckness $b=1$ at dIffeREnt tImes FOr $ \EpsILon = 0.1 $. ThE symbOlS StANd for analytical valuEs WhereaS lineS represent the Results obtAINEd from fiNite DIfFErence method.](fiGure4.ePs){width="8cm"}

![\[vXT-planar\]SCaled Material Energy denSITy $v(x,\tau)$ vS poSitIon (X) in THE sLab at differenT TImes FoR $\epsiloN = 0.1$. ThE symbolS stAnd For AnaLyTical valuEs whereaS lInEs RePreSent tHE results ObTaiNeD frOm finITe diffErencE metHoD.](fIGurE5.eps){widTH="8cM"}

tHe fiRsT dErivAtiVeS w.r.t. pOsitIOn oF the anaLytical raDiaTIon aNd MaTerial eNergy density aRe Plotted in FIgS. \[upRimext-PLAnar\] and \[vPrimext-planar\]. As the radiaTIon and mAteRial eNergY densitieS deCrease WitH X, the deRivatiVe has NeGatIVE valuES. thE deRiVative has a GREatEr negAtIve vAlue at tHe left compared to thE RigHt zone. As both rAdiAtioN ANd MatERiAL enErGY deNSIties keep on incrEasing with TiME dUe to radiatIOn dIfFusion, mAgnitudE of thE GradienT decreaseS for the leFt And iNCReaSes for the rIght sideS. The gradiENt of bOTh RadiaTioN and maTeRiaL enerGy densITieS obtaIn a conStAnt valUe of $\fRaC{-3}{3+4\sqrt{3}}=-0.30217$ aftEr infinite time showing thAt therE is a cOnsTant leakaGe oF
 distance keeps on attaini ng lineari ty. A t i nfi ni te t ime,  the steady st a te v alues are linear withposit io n  asg iv en by  Eq. \[ i nf i n ite -t im e\] .S im ilarl y,the mat erial ener gyde nsity initia l ly  exhibitssli ght non-line arvariat io n a n d fin all y att ains t h e line arity (Fi g. \[vxt- p lanar\] ) .  T he n on-linear variati o na t early stages  occur sd ue t o n etabsorption  o f ene r gy by t h ei n i tia l ly cold mater ial (as u(x , 0)= v(x,0) =0 ).I nitial ly, t he mat erial energ y de nsity isfoundt o lag b e hind th e radi ati onener g yde nsi ti e s a n dfin a lly  equilib ra te  as t imep r o c eeds  (b eyon d $\t au$=10). In t his  wor k , a ll th e res ults  h ave b een ob taine dby consideringcont ributionfro mthe  f irst3 0 root s o f t he tran scenden t aleq u a t io n. The value of op ac i t y$\kappa$  is ch o se nt o be 100  a nd$\ep s i lon $  equ a ls  0.1. Fo r a he a twa ve trav el ing in to  athi n pla t e an d comp osite pl anars lab, a similar linear variat i on i nt empe rat ure with di stan c e wa s ob s er ved thoug h dif fe r en c e existed in the sp ac e andtimedependent beh aviour due t o  heat co nduc t io n  approximation  [@Sa rkar],[@Su n ].

![\[ uxt-p lanar\]S caled rad i a tion ene rgy  de nsi ty$ u (x ,\tau)$ vs po s i tion  ( x) in t heslab of  sc ale d t hic kn ess $b=1$  at diff er en tti mes  for$  \epsilo n= 0 .1  $.  Thes ymbols  stan d fo ran a lyt ical va l ue s wher ea sline s r ep resen t th e  re sults o btained f rom fini te  d ifferen ce method.](f ig ure4.eps){ wi dth ="8cm" } 
![\[vxt -planar\]Scaled materia l  energy  de nsity  $v( x,\tau)$vspositi on( x) inthe sl ab at  d iff e r ent t i m es  fo r$\epsilon= 0.1 $. Th esymb ols sta nd for analyticalv alu es whereas li nes  rep r e se ntt he res ul t s o b t ained from fini te differe nc e  m ethod.](fi g ure 5. eps){wi dth="8c m"}
T he firs t derivat ives w.r. t.  pos i t ion  of the an alytical  radiatio n  andm at erial  en ergy d en sit y are  plott e d i n Fig s. \[u pr imext- plana r\ ] and \[ vprimext-planar\]. As t he rad iatio n a nd materi ale ner gy densit iesdecrease w ith  x,  theder i vativ e ha s  n ega t ive v alue s . The der i va tiv e ha s a greater n e gat ive v alu e  at th e le ft compared to th e  right zone. A s bo t h  ra dia t ionan d material ene rgy  d e n sities k ee p on increa sing wit ht ime d ue toradiat ion dif f u si o n, mag nitu deof the gr adi en t  decrea se sf or the  lef tand in crease s  for t he right sides.The g r a dient ofbothra diation andmaterial e nergy densi ties o btai n a c onstant  v alue o f $ \f rac{-3}{3+ 4 \sqrt{3}} =-0.3 0217$ a ft er i nfi nite t imes h owing  tha tthe re is a c o n st a nt  l e aka ge o f
 distance_keeps on_attaining linearity. At infinite_time, the_steady_state values_are_linear with position_as given by_Eq. \[infinite-time\]. Similarly, the_material energy density_initially_exhibits slight non-linear variation and finally attains the linearity (Fig. \[vxt-planar\]). The non-linear variation_at_early stages_occurs_due_to net absorption of energy_by the initially cold material_(as u(x,0)=v(x,0)=0)._Initially, the material energy density is found to_lag_behind the radiation_energy densities and finally equilibrate as time proceeds (beyond_$\tau$=10). In this work, all the_results have been_obtained_by_considering contribution from the_first 30 roots of the transcendental_equation. The value of opacity $\kappa$_is chosen to be 100 and $\epsilon_$ equals 0.1. For a heat_wave traveling into a thin_plate and_composite planar slab, a similar_linear variation in_temperature with_distance was observed_though difference existed in the space_and time dependent_behaviour due to heat conduction approximation_[@Sarkar],[@Sun].

![\[uxt-planar\]Scaled_radiation energy density_$u(x,\tau)$_vs_position (x)_in the slab_of_scaled thickness_$b=1$_at different times for $ \epsilon_=_0.1 $. The symbols stand for analytical_values whereas lines represent_the_results obtained from finite_difference method.](figure4.eps){width="8cm"}

![\[vxt-planar\]Scaled material energy density_$v(x,\tau)$ vs position (x) in the_slab at_different times_for $\epsilon = 0.1$. The symbols stand for analytical values whereas_lines represent the results obtained from_finite difference method.](figure5.eps){width="8cm"}

The first_derivatives w.r.t._position_of the analytical_radiation_and material_energy density are plotted in Figs. \[uprimext-planar\]_and \[vprimext-planar\]._As the radiation and material energy_densities decrease with x,_the_derivative has negative values. The derivative_has a greater negative value at_the left compared to the_right_zone._As both radiation and material_energy densities keep on increasing with_time due to_radiation diffusion, magnitude of the gradient decreases_for_the left and increases for the_right_sides. The gradient of both radiation_and_material_energy densities obtain a constant_value of $\frac{-3}{3+4\sqrt{3}}=-0.30217$ after infinite time_showing that there is a constant leakage of