 With (\[tauAB\]), we see that both $\tau_{j+1}^{A}$ and $\tau_{j+1}^{B}$ are $\mathcal{O}\left(  \Delta_{j}^{2}\right)  $.

In a similar manner, we obtain for the symmetrized matrix (\[T2\]) $\tau_{j+1}^{\psi}=\mathcal{O}\left(  \Delta_{j}^{3}\right)  $ and $\tau
_{j+1}^{\psi^{\prime}}=\mathcal{O}\left(  \Delta_{j}^{3}\right)  $. More precisely, the calculation yields for a constant $m^{\ast}$ $$\begin{aligned}
\tau_{j+1}^{\psi}  &  =\frac{1}{24}\left(  k_{j}^{2}\right)  ^{\prime}\psi
_{j}\Delta_{j}^{3}+\mathcal{O}\left(  \Delta_{j}^{4}\right) ,\nonumber\\
\tau_{j+1}^{\psi^{\prime}}  &  =-\frac{1}{12}\left(  k_{j}^{2}\right)
^{\prime\prime}\psi_{j}\Delta_{j}^{3}-\frac{1}{24}\left(  k_{j}^{2}\right)
^{\prime}\psi_{j}^{\prime}\Delta_{j}^{3}+\mathcal{O}\left(  \Delta_{j}%
^{4}\right)\end{aligned}$$ (and a somewhat more complicated expression for a position dependent effective mass). This means that $\tau_{j+1}^{A}$ and $\tau_{j+1}^{B}$ are now $\mathcal{O}\left(  \Delta_{j}^{3}\right)  $.

Piecewise linear potential approximation
----------------------------------------

For computing the LDEs (\[tau\]) and (\[taus\]) of the Airy function approach, we proceed in a manner similar as above. Equations (\[airy1\]) and (\[airy2\]) give the relation between values $\psi_{j}$, $\psi_{j}^{\prime}$ at $z_{j}$ and the numerical result $\psi_{j+1}$, $\psi_{j+1}^{\prime}$ obtained at $z_{j+1}$ from the Airy function approach: $$\begin{aligned}
\psi_{j+1}  &  =\frac{1}{D_{j}}\left[  \mathrm{
With (\[tauAB\ ]), we see that both $ \tau_{j+1}^{A}$ and $ \tau_{j+1}^{B}$ are $ \mathcal{O}\left (  \Delta_{j}^{2}\right)   $. 

 In a similar manner, we obtain for the symmetrized matrix (\[T2\ ]) $ \tau_{j+1}^{\psi}=\mathcal{O}\left (  \Delta_{j}^{3}\right)   $ and $ \tau 
 _ { j+1}^{\psi^{\prime}}=\mathcal{O}\left (  \Delta_{j}^{3}\right)   $. More precisely, the calculation render for a changeless $ m^{\ast}$ $ $ \begin{aligned } 
 \tau_{j+1}^{\psi }   &   = \frac{1}{24}\left (  k_{j}^{2}\right)   ^{\prime}\psi 
 _ { j}\Delta_{j}^{3}+\mathcal{O}\left (  \Delta_{j}^{4}\right), \nonumber\\ 
 \tau_{j+1}^{\psi^{\prime } }   &   = -\frac{1}{12}\left (  k_{j}^{2}\right) 
 ^{\prime\prime}\psi_{j}\Delta_{j}^{3}-\frac{1}{24}\left (  k_{j}^{2}\right) 
 ^{\prime}\psi_{j}^{\prime}\Delta_{j}^{3}+\mathcal{O}\left (  \Delta_{j}% 
 ^{4}\right)\end{aligned}$$ (and a somewhat more complicated expression for a position subject effective mass). This means that $ \tau_{j+1}^{A}$ and $ \tau_{j+1}^{B}$ are immediately $ \mathcal{O}\left (  \Delta_{j}^{3}\right)   $. 

 Piecewise analogue potential approximation 
 ---------------------------------------- 

 For computing the LDEs (\[tau\ ]) and (\[taus\ ]) of the Airy function approach, we continue in a manner similar as above. Equations (\[airy1\ ]) and (\[airy2\ ]) pass the relation between values $ \psi_{j}$, $ \psi_{j}^{\prime}$ at $ z_{j}$ and the numerical result $ \psi_{j+1}$, $ \psi_{j+1}^{\prime}$ obtained at $ z_{j+1}$ from the Airy function access: $ $ \begin{aligned } 
 \psi_{j+1 }   &   = \frac{1}{D_{j}}\left [   \mathrm {
 Wihh (\[tauAB\]), we see that botm $\tau_{j+1}^{A}$ and $\tau_{j+1}^{Y}$ are $\methcal{O}\meft(  \Deuta_{j}^{2}\right)  $.

In a similar mannxr, ww obtqin for the symmetrizea matrix (\[N2\]) $\tau_{j+1}^{\psi}=\nathral{O}\left(  \Delta_{j}^{3}\cjght)  $ and $\tah
_{m+1}^{\psi^{\pxine}}=\mathcal{O}\left(  \Delta_{j}^{3}\richt)  $. More prewirepy, the calculation yields for a conftant $m^{\ssh}$ $$\begin{aligned}
\jau_{j+1}^{\pxy}  &  =\frac{1}{24}\left(  k_{j}^{2}\right)  ^{\prime}\psi
_{j}\Delta_{n}^{3}+\mathcan{O}\left(  \Delta_{k}^{4}\right) ,\nonumber\\
\tau_{j+1}^{\psi^{\primf}}  &  =-\frac{1}{12}\left(  k_{j}^{2}\right)
^{\orime\prime}\pwi_{j}\Dqota_{j}^{3}-\frac{1}{24}\left(  k_{j}^{2}\right)
^{\pgnme}\psi_{j}^{\primg}\Delta_{j}^{3}+\mathcal{O}\left(  \Delta_{j}%
^{4}\right)\dnd{alngned}$$ (and a simeavat more conplicwted expresslpn for a posiyion dependent efhectuve mass). This means tiat $\tau_{j+1}^{A}$ and $\tau_{j+1}^{B}$ wre now $\mdtkcal{O}\left(  \Delta_{j}^{3}\righr)  $.

Piecgwise lindqr oottntmal approdimetion
----------------------------------------

For cojputing the LDEs (\[tau\]) and (\[taus\]) pf nne Airy funcfion a[pwoach, we proceed in a manner similar as abkve. Equations (\[airy1\]) and (\[qiry2\]) give the relatioj between values $\psi_{j}$, $\psi_{j}^{\prime}$ at $z_{j}$ and the numerical revult $\'sk_{j+1}$, $\ksl_{m+1}^{\priow}$ lbtained at $z_{j+1}$ from the Airy function approacr: $$\btgik{aligned}
\psi_{j+1}  &  =\nrac{1}{D_{j}}\left[  \mathrm{
With (\[tauAB\]), we see that both $\tau_{j+1}^{A}$ are \Delta_{j}^{2}\right) $. a similar manner, matrix $\tau_{j+1}^{\psi}=\mathcal{O}\left( \Delta_{j}^{3}\right) $ $\tau _{j+1}^{\psi^{\prime}}=\mathcal{O}\left( \Delta_{j}^{3}\right) More precisely, the calculation yields for constant $m^{\ast}$ $$\begin{aligned} \tau_{j+1}^{\psi} & =\frac{1}{24}\left( k_{j}^{2}\right) ^{\prime}\psi _{j}\Delta_{j}^{3}+\mathcal{O}\left( \Delta_{j}^{4}\right) ,\nonumber\\ \tau_{j+1}^{\psi^{\prime}} & k_{j}^{2}\right) ^{\prime\prime}\psi_{j}\Delta_{j}^{3}-\frac{1}{24}\left( k_{j}^{2}\right) ^{\prime}\psi_{j}^{\prime}\Delta_{j}^{3}+\mathcal{O}\left( \Delta_{j}% ^{4}\right)\end{aligned}$$ (and a somewhat more complicated expression for position effective This that $\tau_{j+1}^{A}$ and $\tau_{j+1}^{B}$ are now $\mathcal{O}\left( \Delta_{j}^{3}\right) $. Piecewise linear potential approximation ---------------------------------------- For computing LDEs (\[tau\]) and (\[taus\]) of the Airy function we proceed in a similar as above. Equations (\[airy1\]) (\[airy2\]) the relation values $\psi_{j}^{\prime}$ $z_{j}$ and the result $\psi_{j+1}$, $\psi_{j+1}^{\prime}$ obtained at $z_{j+1}$ from the Airy function approach: $$\begin{aligned} \psi_{j+1} & =\frac{1}{D_{j}}\left[ \mathrm{
 With (\[tauAB\]), we see that both $\tau_{J+1}^{A}$ and $\tau_{j+1}^{B}$ Are $\maThcAl{O}\LeFt(  \DeLta_{j}^{2}\Right)  $.

In a similaR MannEr, we obtain for the symmetRized MaTRix (\[T2\]) $\TAu_{J+1}^{\psi}=\mAthcal{O}\LEfT(  \dEltA_{j}^{3}\RiGht)  $ AnD $\TaU
_{j+1}^{\psi^{\PriMe}}=\mathcAl{O}\left(  \DelTa_{j}^{3}\RiGht)  $. More preciSElY, the calculAtiOn yields for a ConStant $m^{\AsT}$ $$\beGIn{aliGneD}
\tau_{j+1}^{\Psi}  &  =\fraC{1}{24}\Left(  k_{j}^{2}\Right)  ^{\primE}\pSI
_{j}\DeltA_{J}^{3}+\mathcaL{o}\LeFt(  \DeLta_{j}^{4}\right) ,\nonumber\\
\TAu_{J+1}^{\Psi^{\prime}}  &  =-\frac{1}{12}\leFt(  k_{j}^{2}\riGhT)
^{\PrIME\prIme}\Psi_{j}\Delta_{j}^{3}-\FrAc{1}{24}\lefT(  K_{j}^{2}\right)
^{\PRiME}\PSi_{j}^{\PRime}\Delta_{j}^{3}+\matHcal{O}\left(  \DeLTa_{j}%
^{4}\Right)\eNd{AliGNed}$$ (and A someWhAT moRe complicatEd exPression fOr a posITion depENdent efFectivE maSs). THis mEAnS tHat $\TaU_{J+1}^{A}$ aND $\tAu_{j+1}^{b}$ Are Now $\mathcAl{o}\lEft(  \DeLta_{j}^{3}\RIGHT)  $.

PieCewIse lInear Potential apprOxiMatiON
----------------------------------------

FoR compUting The LdES (\[tau\]) aNd (\[taus\]) Of the aiRy function approAch, wE proceed iN a mAnNer SiMilar AS above. equAtiOns (\[airy1\]) And (\[airy2\]) GIve ThE RELaTion between values $\pSi_{J}$, $\PSi_{J}^{\prime}$ at $Z_{j}$ and tHE nUmERical resUlT $\psI_{j+1}$, $\psI_{J+1}^{\Prime}$ ObtaINeD at $z_{j+1}$ froM the AiRY fUnCtion apPrOach: $$\beGiN{alIgnEd}
\psi_{J+1}  &  =\Frac{1}{d_{j}}\left[  \Mathrm{
 With (\[tauAB\]), we seethat both$\tau _{j +1} ^{ A}$and$\tau_{j+1}^{B } $ ar e $\mathcal{O}\left( \Delt a_ { j}^{ 2 }\ right )  $.
I na sim il ar  ma nn e r,  we o bta in forthe symmet riz ed  matrix (\[T 2 \] ) $\tau_{j +1} ^{\psi}=\mat hca l{O}\l ef t(\Delt a_{ j}^{3 }\righ t )  $ a nd $\tau_{ j +1}^{\ p si^{\pr i m e} }=\m athcal{O}\left( \ De l ta_{j}^{3}\rig ht)  $ .M or e pre cis ely, the c al culat i on yiel d sf o r  ac onstant $m^{\ ast}$ $$\be g in{ aligne d} 
\t a u_{j+1 }^{\p si }   &   =\frac{1} {24} \left(  k _{j}^{ 2 }\right )   ^{\pr ime}\p si_{j }\De l ta _{ j}^ {3 } +\m a th cal { O}\ left(  \ De lt a_{j} ^{4} \ r i g ht),\n onum ber\\ 
\tau_{j+1}^{ \ps i^{\ p rim e}} &  =- \fra c{ 1}{12 }\left (  k_ {j }^{2}\right)
^{ \pri me\prime} \ps i_ {j} \D elta_ { j}^{3} -\f rac {1}{24} \left(k_{ j} ^ { 2 }\ right)
^{\prime}\p si _ { j} ^{\prime }\Delt a _{ j} ^ {3}+\mat hc al{ O}\l e f t(  \ Delt a _{ j}%
^{4} \right ) \e nd {aligne d} $$ (an da s ome whatm orecompli cated ex press i on for a posit i on dependente ff e c ti v e ma ss) . This mean s th a t $\ tau_ { j+ 1}^ { A}$ a nd $\ ta u _{ j +1}^{B}$ are now $\ ma thcal{ O}\le ft(  \Delta_{ j}^{3}\rig h t )   $.

Pi ecew i se linear potenti al ap proximatio n 
------- ----- -------- --------- - - -------- -
For  co mpu t i ng  the LDEs (\[ t a u\])  a nd (\[t aus \]) ofthe  Ai ryfun ct ion appro ach, wepr oc ee dina man n er simil ar  as  a bov e. Eq u ations  (\[a iry1 \] )a nd(\[airy 2 \] ) give  t he  rel ati on  betw eenv alu es $\ps i_{j}$, $ \ps i _{j} ^{ \p rime}$at $z_{j}$ an dthe numeri ca l r esult$ \ psi_{j+1 }$, $\psi_{j+1}^{\prime } $ obtai ned  at $ z_{j +1}$ from  th e Airy  fu n ctionapproa ch: $ $\ beg i n {alig n e d} 
\p si _{j+1}  &= \fr ac{1} {D _{j} }\left[   \mathrm{
 With_(\[tauAB\]), we_see that both $\tau_{j+1}^{A}$_and $\tau_{j+1}^{B}$_are_$\mathcal{O}\left( _\Delta_{j}^{2}\right)_ $.

In a_similar manner, we_obtain for the symmetrized_matrix (\[T2\]) $\tau_{j+1}^{\psi}=\mathcal{O}\left(__\Delta_{j}^{3}\right)  $ and $\tau
_{j+1}^{\psi^{\prime}}=\mathcal{O}\left(  \Delta_{j}^{3}\right)  $. More precisely, the calculation yields_for_a constant_$m^{\ast}$_$$\begin{aligned}
\tau_{j+1}^{\psi}_ &  =\frac{1}{24}\left( _k_{j}^{2}\right)  ^{\prime}\psi
_{j}\Delta_{j}^{3}+\mathcal{O}\left(  \Delta_{j}^{4}\right)_,\nonumber\\
\tau_{j+1}^{\psi^{\prime}} _&  =-\frac{1}{12}\left(  k_{j}^{2}\right)
^{\prime\prime}\psi_{j}\Delta_{j}^{3}-\frac{1}{24}\left(  k_{j}^{2}\right)
^{\prime}\psi_{j}^{\prime}\Delta_{j}^{3}+\mathcal{O}\left( _\Delta_{j}%
^{4}\right)\end{aligned}$$_(and a somewhat_more complicated expression for a position dependent effective mass)._This means that $\tau_{j+1}^{A}$ and $\tau_{j+1}^{B}$_are now $\mathcal{O}\left(__\Delta_{j}^{3}\right)_ $.

Piecewise linear potential_approximation
----------------------------------------

For computing the LDEs (\[tau\]) and_(\[taus\]) of the Airy function approach,_we proceed in a manner similar as_above. Equations (\[airy1\]) and (\[airy2\]) give_the relation between values $\psi_{j}$,_$\psi_{j}^{\prime}$ at_$z_{j}$ and the numerical result_$\psi_{j+1}$, $\psi_{j+1}^{\prime}$ obtained_at $z_{j+1}$_from the Airy_function approach: $$\begin{aligned}
\psi_{j+1}  & _=\frac{1}{D_{j}}\left[  \mathrm{