(* ::Package:: *)

(************************************************************************)
(* This file was generated automatically by the Mathematica front end.  *)
(* It contains Initialization cells from a Notebook file, which         *)
(* typically will have the same name as this file except ending in      *)
(* ".nb" instead of ".m".                                               *)
(*                                                                      *)
(* This file is intended to be loaded into the Mathematica kernel using *)
(* the package loading commands Get or Needs.  Doing so is equivalent   *)
(* to using the Evaluate Initialization Cells menu command in the front *)
(* end.                                                                 *)
(*                                                                      *)
(* DO NOT EDIT THIS FILE.  This entire file is regenerated              *)
(* automatically each time the parent Notebook file is saved in the     *)
(* Mathematica front end.  Any changes you make to this file will be    *)
(* overwritten.                                                         *)
(************************************************************************)



(* ::Input::Initialization:: *)
wdir=SetDirectory["C:\\drorbn\\AcademicPensieve\\Projects\\PPSA"];
NotebookOpen[wdir<>"\\MakeVSnips.nb"];


(* ::Input::Initialization:: *)
HL[\[ScriptCapitalE]_] := Style[\[ScriptCapitalE], Background -> Yellow];


(* ::Input::Initialization:: *)
$T\[HBar]D=3; $T\[Epsilon]D=2; \[Epsilon] /: \[Epsilon]^d_. /; d>$T\[Epsilon]D :=0;
(* $T\[Epsilon]D can't be \[Infinity] at least because of Quesne. Can't be \[LessEqual]1 at least because of the explicit \[Epsilon]^2 in S\[DoubleStruckCapitalD]$g. *)
SetAttributes[{SS, SST}, HoldAll];
SS[\[ScriptCapitalE]_] := Block[{\[HBar],\[Epsilon]}, (* Shielded Series *)
Collect[Normal@Series[\[ScriptCapitalE], {\[HBar],0,$T\[HBar]D}],\[HBar],Together] ];
SST[\[ScriptCapitalE]_] := Block[{\[HBar],\[Epsilon]},Collect[Normal@Series[\[ScriptCapitalE]/.{Subscript[T, i_]-> E^(\[HBar] Subscript[t, i]/2),T-> E^(\[HBar] t/2)}, {\[HBar],0,$T\[HBar]D}],\[HBar],Together] ];
Simp[\[ScriptCapitalE]_, op_]:= Collect[\[ScriptCapitalE], _CU|_QU, op];
Simp[\[ScriptCapitalE]_]:=Simp[\[ScriptCapitalE],Collect[Normal@Series[#, {\[HBar],0,$T\[HBar]D}],\[HBar],Expand]&];
SimpT[\[ScriptCapitalE]_]:= Collect[\[ScriptCapitalE], _CU|_QU, Collect[Normal@Series[#/.{Subscript[T, i_]-> E^(\[HBar] Subscript[t, i]/2),T-> E^(\[HBar] t/2)}, {\[HBar],0,$T\[HBar]D}],\[HBar],Expand]&];


(* ::Input::Initialization:: *)
Subscript[DP, \[Alpha]_->Subscript[D, x_],\[Beta]_->Subscript[D, y_]][P_][\[Lambda]_]:=Total[CoefficientRules[P,{\[Alpha],\[Beta]}]/.({m_,n_}->c_):>c D[\[Lambda],{x,m},{y,n}]]


(* ::Input::Initialization:: *)
Unprotect[NonCommutativeMultiply]; Attributes[NonCommutativeMultiply]={};
(NCM=NonCommutativeMultiply)[x_] := x;
NCM[x_,y_,z__] := (x**y)**z;
0**_=_**0=0;
(x_Plus)**y_:=(#**y)&/@x; x_**(y_Plus):=(x**#)&/@y;
B[x_,x_]=0; B[x_,y_] := x**y-y**x;


(* ::Input::Initialization:: *)
DeclareAlgebra[U_Symbol, opts__Rule] := Module[{gp,sr,cp,CE,pow,
gs=Generators /. {opts}, cs=Centrals /. {opts} },
(Subscript[#, U]=U@#)&/@gs;
gp=Alternatives@@gs; gp=gp|Subscript[gp, _];(* gen's pattern *)
sr=Thread[gs->Range@Length@gs]; (* sorting rule *)
cp=Alternatives@@cs; (* cent's pattern *)
CE[\[ScriptCapitalE]_]:= Collect[\[ScriptCapitalE], _U, (Expand[#] /. \[HBar]^d_/;d>$T\[HBar]D:>0)&];
Subscript[U, i_][\[ScriptCapitalE]_] := \[ScriptCapitalE] /. {t:cp:>Subscript[t, i],u_U :> Replace[u, x_ :> Subscript[x, i],1]};
Subscript[U, i_][NCM[]]:=U[];
B[U@Subscript[(x_), i_],U@Subscript[(y_), i_]]:=B[U@Subscript[x, i],U@Subscript[y, i]]=Subscript[U, i]@B[U@x,U@y];
B[U@Subscript[(x_), i_],U@Subscript[(y_), j_]] /; i=!=j := 0;
B[U@y_,U@x_]:=CE[-B[U@x,U@y]];
x_**U[]:=x; U[]**x_:=x;
(a_.*x_U)**(b_.*y_U) := If[a b === 0, 0, CE[a b (x**y)]];
U[xx___,x_] ** U[y_,yy___] := If[OrderedQ[{x,y}/.sr],U[xx,x,y,yy], U@xx**(U@y**U@x+B[U@x,U@y])**U@yy];
U@{c_.*(l:gp)^n_,r___}/;FreeQ[c, gp]:= CE[c U@Table[l,{n}]**U@{r}];
U@{c_.*l:gp,r___} := CE[c U[l]**U@{r}];
U@{c_,r___}/;FreeQ[c, gp]:= CE[c U@{r}];
U@{}=U[];
U@{l_Plus, r___} := CE[U@{#,r}& /@ l];
U@{l_,r___} := U@{Expand[l],r};
U[\[ScriptCapitalE]_NonCommutativeMultiply]:=U/@\[ScriptCapitalE];
Subscript[\[DoubleStruckCapitalO], U][poly_, specs___] := Module[{sp,null,vs,us},
sp=Replace[{specs}, l_List:>Subscript[l, null], {1}];
vs=Join@@(First /@ sp);
us=Join@@(sp /. Subscript[l_, s_] :> (l/.Subscript[x_, i_]:>Subscript[x, s]));
CE[Total[
CoefficientRules[poly, vs] /. (p_->c_) :> c U@(us^p)
]] /. Subscript[x_, null]:>x
 ];
pow[\[ScriptCapitalE]_,0]=U[]; pow[\[ScriptCapitalE]_,n_]:=pow[\[ScriptCapitalE],n-1]**\[ScriptCapitalE];
Subscript[\[DoubleStruckCapitalS], U][\[ScriptCapitalE]_,ss___Rule] := CE@Total[
CoefficientRules[\[ScriptCapitalE],First/@{ss}]/.(p_->c_):>c NCM@@MapThread[pow,{Last/@{ss},p}]];
Subscript[S, i_][c_.*u_U]:= CE[(c/.Subscript[S, i][U,Centrals])DeleteCases[u,Subscript[_, i]]**Subscript[U, i][NCM@@Reverse@Cases[u,Subscript[x_, i]:>S@U@x]]];
]


(* ::Input::Initialization:: *)
DeclareMorphism[m_,U_->V_, ongs_List, oncs_List:{}] := (
Replace[ongs, (g_->img_) :> (m[U[g]]=img), {1}];
m[U[]]=V[];
m[U[Subscript[g_, i_]]] := Subscript[V, i][m[U@g]];
m[U[vs__]] := NCM@@(m/@U/@{vs});
m[\[ScriptCapitalE]_] := Simp[\[ScriptCapitalE] /. oncs /.  u_U:>m[u]]; )


(* ::Input::Initialization:: *)
Subscript[S, i_][\[ScriptCapitalE]_Plus] := Simp[Subscript[S, i]/@\[ScriptCapitalE]];


(* ::Input::Initialization:: *)
DeclareAlgebra[CU,Generators->{y,a,x},Centrals->{t}];
B[Subscript[a, CU],Subscript[y, CU]]=-\[Gamma] Subscript[y, CU]; B[Subscript[x, CU],Subscript[a, CU]]=-\[Gamma] Subscript[x, CU];
B[Subscript[x, CU],Subscript[y, CU]]=2\[Epsilon] Subscript[a, CU]-t CU[];
(S@CU@y=-Subscript[y, CU]; S@Subscript[a, CU]=-Subscript[a, CU]; S@Subscript[x, CU]=-Subscript[x, CU];)
Subscript[S, i_][CU,Centrals]={Subscript[t, i]->-Subscript[t, i]};


(* ::Input::Initialization:: *)
DeclareAlgebra[QU, Generators -> {y,a,x}, Centrals -> {t,T}];
q=SS[E^(\[Gamma] \[Epsilon] \[HBar])]; (*T=SS[\[ExponentialE]^(\[HBar] t/2)];*)
B[Subscript[a, QU],Subscript[y, QU]]=-\[Gamma] Subscript[y, QU]; B[Subscript[x, QU],Subscript[a, QU]]=-\[Gamma] QU@x;
B[Subscript[x, QU],Subscript[y, QU]]=(q-1)QU@{y,x}+Subscript[\[DoubleStruckCapitalO], QU][SS[(1-T^2 E^(-2\[Epsilon] a \[HBar]))/\[HBar]], {a}];
(S@Subscript[y, QU]=Subscript[\[DoubleStruckCapitalO], QU][SS[-T^-2 E^(\[HBar] \[Epsilon] a) y],{a,y}]; S@Subscript[a, QU]=-Subscript[a, QU]; S@Subscript[x, QU]=Subscript[\[DoubleStruckCapitalO], QU][SS[-E^(\[HBar] \[Epsilon] a)x],{a,x}];)
Subscript[S, i_][QU,Centrals]={Subscript[t, i]->-Subscript[t, i], Subscript[T, i]->\!\(\*SubsuperscriptBox[\(T\), \(i\), \(-1\)]\)};


(* ::Input::Initialization:: *)
DeclareMorphism[C\[Theta],CU->CU, {y->-Subscript[x, CU], a->-Subscript[a, CU], x->-Subscript[y, CU]}, {t->-t, T->T^-1}];DeclareMorphism[Q\[Theta],QU->QU, {y->Subscript[\[DoubleStruckCapitalO], QU][SS[-T^-1 E^(\[HBar] \[Epsilon] a) x],{a,x}], a->-Subscript[a, QU], x->Subscript[\[DoubleStruckCapitalO], QU][SS[-T^-1 E^(\[HBar] \[Epsilon] a) y],{a,y}]}, {t->-t, T->T^-1}]


(* ::Input::Initialization:: *)
A\[DoubleStruckCapitalD]$f=\[Gamma]/\[HBar] E^(\[HBar](t/2-(a+\[Gamma])\[Epsilon])) (Cosh[\[HBar](a \[Epsilon]+(\[Gamma] \[Epsilon])/2-t/2)]-Cosh[\[HBar] Sqrt[((t-\[Gamma] \[Epsilon])/2)^2+\[Epsilon] \[Omega]]])/(Sinh[(\[Gamma] \[Epsilon] \[HBar])/2](a^2 \[Epsilon]+a \[Gamma] \[Epsilon]-a t-\[Omega]));


(* ::Input::Initialization:: *)
A\[DoubleStruckCapitalD]$\[Omega]=\[Gamma] CU[y,x]+\[Epsilon] CU[a,a]-(t-\[Gamma] \[Epsilon])CU[a];


(* ::Input::Initialization:: *)
DeclareMorphism[A\[DoubleStruckCapitalD],QU->CU,{a->Subscript[a, CU],x->CU@x, y->Subscript[\[DoubleStruckCapitalS], CU][SS[A\[DoubleStruckCapitalD]$f],a->Subscript[a, CU],\[Omega]->A\[DoubleStruckCapitalD]$\[Omega]]**Subscript[y, CU]}]


(* ::Input::Initialization:: *)
S\[DoubleStruckCapitalD]$g=Sqrt[(Cosh[\[HBar]/2 Sqrt[t^2+\[Gamma]^2 \[Epsilon]^2+4 \[Epsilon] \[CurlyPi]]]-Cosh[\[HBar]/2 (t-(2 a+\[Gamma]) \[Epsilon])])/(Sinh[(\[Gamma] \[Epsilon] \[HBar])/2](t (2 a+\[Gamma])-2 a (a+\[Gamma]) \[Epsilon]+2 \[CurlyPi]) \[HBar]/(2\[Gamma]))];


(* ::Input::Initialization:: *)
S\[DoubleStruckCapitalD]$f=FullSimplify[E^(\[HBar](t/2-\[Epsilon] a)) (S\[DoubleStruckCapitalD]$g /. {a->-a, t->-t})];


(* ::Input::Initialization:: *)
S\[DoubleStruckCapitalD]$\[CurlyPi]=\[Gamma] CU[y,x]+\[Epsilon] CU[a,a]-(t-\[Gamma] \[Epsilon])CU[a]-t \[Gamma] CU[]/2;


(* ::Input::Initialization:: *)
DeclareMorphism[S\[DoubleStruckCapitalD],QU->CU,{a->Subscript[a, CU],
x->Subscript[\[DoubleStruckCapitalS], CU][SS[S\[DoubleStruckCapitalD]$f],a->Subscript[a, CU],\[CurlyPi]->S\[DoubleStruckCapitalD]$\[CurlyPi]]**Subscript[x, CU],
y->Subscript[\[DoubleStruckCapitalS], CU][SS[S\[DoubleStruckCapitalD]$g],a->Subscript[a, CU],\[CurlyPi]->S\[DoubleStruckCapitalD]$\[CurlyPi]]**Subscript[y, CU]
}]


(* ::Input::Initialization:: *)
Subscript[E, q_,n_][x_]:=E^(\!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k = 1\), \(n\)]\ 
\*FractionBox[\(
\*SuperscriptBox[\((1 - q)\), \(k\)] 
\*SuperscriptBox[\(x\), \(k\)]\), \(k \((1 - 
\*SuperscriptBox[\(q\), \(k\)])\)\)]\)); Subscript[E, q_][x_]:=Subscript[E, q,$T\[Epsilon]D][x]


(* ::Input::Initialization:: *)
QU[Subscript[R, i_,j_]] := Subscript[\[DoubleStruckCapitalO], QU][SS[E^(\[HBar] Subscript[b, 1] Subscript[a, 2]) Subscript[E, q][\[HBar] Subscript[y, 1] Subscript[x, 2]]/.Subscript[b, 1]->\[Gamma]^-1 (\[Epsilon] Subscript[a, 1]-Subscript[t, i])], Subscript[{Subscript[y, 1],Subscript[a, 1]}, i], Subscript[{Subscript[a, 2],Subscript[x, 2]}, j]];
QU[\!\(\*SubsuperscriptBox[\(R\), \(i_, j_\), \(-1\)]\)] := Subscript[S, j]@QU[Subscript[R, i,j]];


(* ::Input::Initialization:: *)
\[Rho]@Subscript[y, CU]=\[Rho]@Subscript[y, QU]=({
 {0, 0},
 {\[Epsilon], 0}
});\[Rho]@Subscript[a, CU]=\[Rho]@Subscript[a, QU]=({
 {\[Gamma], 0},
 {0, 0}
});
\[Rho]@Subscript[x, CU]=({
 {0, \[Gamma]},
 {0, 0}
}); \[Rho]@Subscript[x, QU]=SS@\!\(\*
TagBox[
RowBox[{"(", "", GridBox[{
{"0", 
RowBox[{
RowBox[{"(", 
RowBox[{"1", "-", 
SuperscriptBox["E", 
RowBox[{
RowBox[{"-", "\[Gamma]"}], " ", "\[Epsilon]", " ", "\[HBar]"}]]}], ")"}], "/", 
RowBox[{"(", 
RowBox[{"\[Epsilon]", " ", "\[HBar]"}], ")"}]}]},
{"0", "0"}
},
GridBoxAlignment->{"Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}, "Items" -> {}, "ItemsIndexed" -> {}},
GridBoxSpacings->{"Columns" -> {Offset[0.27999999999999997`], {Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {Offset[0.2], {Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}, "Items" -> {}, "ItemsIndexed" -> {}}], "", ")"}],
Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]\);
\[Rho][E^\[ScriptCapitalE]_]:=MatrixExp[\[Rho][\[ScriptCapitalE]]];
\[Rho][\[ScriptCapitalE]_]:=(\[ScriptCapitalE]/.{t->\[Gamma] \[Epsilon],T->E^(\[HBar] \[Gamma] \[Epsilon]/2)}/.(U:CU|QU)[u___]:>Fold[Dot,({
 {1, 0},
 {0, 1}
}),\[Rho]/@U/@{u}])


(* ::Input::Initialization:: *)
Subscript[SS, \[Epsilon]][\[ScriptCapitalE]_] := Block[{\[Epsilon]},Collect[Normal@Series[\[ScriptCapitalE], {\[Epsilon],0,$T\[Epsilon]D}],\[Epsilon],Together]];  (* Shielded \[Epsilon]-Series *)
C\[CapitalLambda][t1_,y1_,a1_,x1_,\[Xi]1_,\[Eta]1_,\[Delta]_]:=Module[
{eqn,d,b,c,sol,\[Lambda],q,\[Nu], \[Xi], \[Eta]},
eqn=\[Rho][E^(\[Xi] Subscript[x, CU])].\[Rho][E^(\[Eta] Subscript[y, CU])]==\[Rho][E^(d Subscript[y, CU])].\[Rho][E^(c(t CU[] - 2\[Epsilon] Subscript[a, CU]))].\[Rho][E^(b Subscript[x, CU])];
sol=Solve[Thread[Flatten/@eqn], {d,b,c}][[1]] /. C[1]->0;
\[Lambda]=Simplify[E^\!\(TraditionalForm\`\(-\[Eta]\)\ y - \[Xi]\ x + \[Eta]\ \[Xi]\ t\) Subscript[SS, \[Epsilon]][E^\!\(TraditionalForm\`c\ t\  + \ d\ y\  - \ 2  \[Epsilon]\ c\ a\  + \ b\ x\)/.sol]];
q=\!\(TraditionalForm\`
\*SuperscriptBox[\(E\), \(\[Nu] \((\(-t\)\ \[Xi]\ \[Eta] + \[Eta]\ y + \[Xi]\ x + \[Delta]\ y\ x)\)\)]\);
Collect[\[Nu] q^-1 Subscript[DP, \[Xi]->Subscript[D, x],\[Eta]->Subscript[D, y]][\[Lambda]][q] /. \[Nu]->(1+t \[Delta])^-1,\[Epsilon],Simplify] /. {t->t1,y->y1,a->a1,x->x1,\[Xi]->\[Xi]1, \[Eta]->\[Eta]1}
];


(* ::Input::Initialization:: *)
SetAttributes[C\[DoubleStruckCapitalO], Orderless];
CU@C\[DoubleStruckCapitalO][specs___, \[DoubleStruckCapitalE][L_,Q_,P_]] := Subscript[\[DoubleStruckCapitalO], CU][SS[E^(L+Q) P], specs]


(* ::Input::Initialization:: *)
Subscript[SW, Subscript[x, i_],Subscript[a, j_]][C\[DoubleStruckCapitalO][Subscript[{lh___,Subscript[x, i_],Subscript[a, j_],rh___}, s_],more___,\[DoubleStruckCapitalE][L_,Q_,P_]]] := C\[DoubleStruckCapitalO][Subscript[{lh,Subscript[a, j],Subscript[x, i],rh}, s],more,
With[{q=E^(-\[Gamma] \[Alpha]) \[Xi] Subscript[x, i]+\[Alpha] Subscript[a, j]},
\[DoubleStruckCapitalE][L,E^(-\[Gamma] \[Alpha]) \[Xi] Subscript[x, i]+(Q/.Subscript[x, i]->0), E^-q Subscript[DP, Subscript[x, i]->Subscript[D, \[Xi]],Subscript[a, j]->Subscript[D, \[Alpha]]][P][E^q]] /. {\[Alpha]->\!\(
\*SubscriptBox[\(\[PartialD]\), 
SubscriptBox[\(a\), \(j\)]]L\), \[Xi]->\!\(
\*SubscriptBox[\(\[PartialD]\), 
SubscriptBox[\(x\), \(i\)]]Q\)}]]


(* ::Input::Initialization:: *)
Subscript[SW, Subscript[x, i_],Subscript[y, j_]->k_][C\[DoubleStruckCapitalO][Subscript[{lh___,Subscript[x, i_],Subscript[y, j_],rh___}, s_],more___,\[DoubleStruckCapitalE][L_,Q_,P_]]] := C\[DoubleStruckCapitalO][Subscript[{lh,Subscript[y, k],Subscript[a, k],Subscript[x, k],rh}, s],more,
With[{q=\[Nu](\[Xi] Subscript[x, k]+\[Eta] Subscript[y, k]+\[Delta] Subscript[x, k] Subscript[y, k]-Subscript[t, k] \[Xi] \[Eta])},
\[DoubleStruckCapitalE][L,q+(Q/.Subscript[x, i]|Subscript[y, j]->0), E^-q Subscript[DP, Subscript[x, i]->Subscript[D, \[Xi]],Subscript[y, j]->Subscript[D, \[Eta]]][P][C\[CapitalLambda][Subscript[t, k],Subscript[y, k],Subscript[a, k],Subscript[x, k],\[Xi],\[Eta],\[Delta]]E^q]] /. \[Nu]->(1+Subscript[t, k]\[Delta])^-1/. {\[Xi]->(\!\(
\*SubscriptBox[\(\[PartialD]\), 
SubscriptBox[\(x\), \(i\)]]Q\)/.Subscript[y, j]->0), \[Eta]->(\!\(
\*SubscriptBox[\(\[PartialD]\), 
SubscriptBox[\(y\), \(j\)]]Q\)/.Subscript[x, i]->0), \[Delta]->\!\(
\*SubscriptBox[\(\[PartialD]\), \(
\*SubscriptBox[\(x\), \(i\)], 
\*SubscriptBox[\(y\), \(j\)]\)]Q\)}]]


(* ::Input::Initialization:: *)
Q\[CapitalLambda][T_,y1_,a1_,x1_,\[Xi]1_,\[Eta]1_,\[Delta]_]:=Module[
{adx,G,F,f,unowns,bas,eqns,sol,\[Lambda],q,\[Nu], \[Xi], \[Eta],t},
adx[\[ScriptCapitalE]_]:=Simp[Subscript[x, QU] ** \[ScriptCapitalE] - \[ScriptCapitalE] ** Subscript[x, QU]];
G=Simp[NestList[adx, Subscript[y, QU], $T\[Epsilon]D+1].Table[\[Xi]^k/k!, {k,0, $T\[Epsilon]D+1}]];
F=Sum[Subscript[f, l,i,j,k][\[Eta]]\[Epsilon]^l QU@{y^i,a^j,x^k},{l,0,$T\[Epsilon]D},{i,0, l},{j,0, l},{k,0, Min[l,2l-i-j]}];
unowns=Cases[F, Subscript[f, ___][\[Eta]],\[Infinity]];
bas=Union@@Table[\[Epsilon]^l Cases[Coefficient[F,\[Epsilon],l],_QU,\[Infinity]], {l,0,$T\[Epsilon]D}];
eqns=Flatten[{(Coefficient[F-QU[], #]/.\[Eta]->0)==0, Expand[Coefficient[Simp[F**G-Subscript[y, QU]**F-\!\(
\*SubscriptBox[\(\[PartialD]\), \(\[Eta]\)]F\)], #]]==0}& /@ bas];
{sol}=DSolve[eqns, unowns, \[Eta]];
\[Lambda]=Collect[F/.sol/.{E^_->1, QU->Times}, \[Epsilon], Simplify];
q=\!\(TraditionalForm\`
\*SuperscriptBox[\(E\), \(\[Nu] \((\(-t\)\ \[Xi]\ \[Eta] + \[Eta]\ y + \[Xi]\ x + \[Delta]\ y\ x)\)\)]\);
Collect[\[Nu] q^-1 Subscript[DP, \[Xi]->Subscript[D, x],\[Eta]->Subscript[D, y]][\[Lambda]][q] /. \[Nu]->(1+t \[Delta])^-1/.t->(T^2-1)/\[HBar],\[Epsilon],Simplify] /. {y->y1,a->a1,x->x1,\[Xi]->\[Xi]1, \[Eta]->\[Eta]1}
];



