(* ::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.                                                         *)
(************************************************************************)



SetDirectory["C:\\drorbn\\AcademicPensieve\\Projects\\V1\\OneCo-1606"];
<< OneCo.m


bas=Cases[LBasis[{j,k}],\[Epsilon]*_]


Sum[Subscript[f, i][Subscript[b, j],Subscript[b, k]]bas[[i]], {i,Length@bas}]


\[Rho]0[j_,k_] :=\[Epsilon] Subscript[c, j] Subscript[f, 1][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[c, k] Subscript[f, 2][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[u, j] Subscript[f, 3][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[u, k] Subscript[f, 4][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[w, j] Subscript[f, 5][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[w, k] Subscript[f, 6][Subscript[b, j],Subscript[b, k]]+\[Epsilon] 
\!\(\*SubsuperscriptBox[\(c\), \(j\), \(2\)]\) Subscript[f, 7][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[c, j] Subscript[c, k] Subscript[f, 8][Subscript[b, j],Subscript[b, k]]+\[Epsilon] 
\!\(\*SubsuperscriptBox[\(c\), \(k\), \(2\)]\) Subscript[f, 9][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[c, j] Subscript[u, j] Subscript[f, 10][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[c, k] Subscript[u, j] Subscript[f, 11][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[c, j] Subscript[u, k] Subscript[f, 12][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[c, k] Subscript[u, k] Subscript[f, 13][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[c, j] Subscript[w, j] Subscript[f, 14][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[c, k] Subscript[w, j] Subscript[f, 15][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[u, j] Subscript[w, j] Subscript[f, 16][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[u, k] Subscript[w, j] Subscript[f, 17][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[c, j] Subscript[w, k] Subscript[f, 18][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[c, k] Subscript[w, k] Subscript[f, 19][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[u, j] Subscript[w, k] Subscript[f, 20][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[u, k] Subscript[w, k] Subscript[f, 21][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[c, j] Subscript[u, j] Subscript[w, j] Subscript[f, 22][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[c, k] Subscript[u, j] Subscript[w, j] Subscript[f, 23][Subscript[b, j],Subscript[b, k]]+\[Epsilon] 
\!\(\*SubsuperscriptBox[\(u\), \(j\), \(2\)]\) Subscript[w, j] Subscript[f, 24][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[c, j] Subscript[u, k] Subscript[w, j] Subscript[f, 25][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[c, k] Subscript[u, k] Subscript[w, j] Subscript[f, 26][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[u, j] Subscript[u, k] Subscript[w, j] Subscript[f, 27][Subscript[b, j],Subscript[b, k]]+\[Epsilon] 
\!\(\*SubsuperscriptBox[\(u\), \(k\), \(2\)]\) Subscript[w, j] Subscript[f, 28][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[u, j] 
\!\(\*SubsuperscriptBox[\(w\), \(j\), \(2\)]\) Subscript[f, 29][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[u, k] 
\!\(\*SubsuperscriptBox[\(w\), \(j\), \(2\)]\) Subscript[f, 30][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[c, j] Subscript[u, j] Subscript[w, k] Subscript[f, 31][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[c, k] Subscript[u, j] Subscript[w, k] Subscript[f, 32][Subscript[b, j],Subscript[b, k]]+\[Epsilon] 
\!\(\*SubsuperscriptBox[\(u\), \(j\), \(2\)]\) Subscript[w, k] Subscript[f, 33][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[c, j] Subscript[u, k] Subscript[w, k] Subscript[f, 34][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[c, k] Subscript[u, k] Subscript[w, k] Subscript[f, 35][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[u, j] Subscript[u, k] Subscript[w, k] Subscript[f, 36][Subscript[b, j],Subscript[b, k]]+\[Epsilon] 
\!\(\*SubsuperscriptBox[\(u\), \(k\), \(2\)]\) Subscript[w, k] Subscript[f, 37][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[u, j] Subscript[w, j] Subscript[w, k] Subscript[f, 38][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[u, k] Subscript[w, j] Subscript[w, k] Subscript[f, 39][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[u, j] 
\!\(\*SubsuperscriptBox[\(w\), \(k\), \(2\)]\) Subscript[f, 40][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[u, k] 
\!\(\*SubsuperscriptBox[\(w\), \(k\), \(2\)]\) Subscript[f, 41][Subscript[b, j],Subscript[b, k]]+\[Epsilon] 
\!\(\*SubsuperscriptBox[\(u\), \(j\), \(2\)]\) 
\!\(\*SubsuperscriptBox[\(w\), \(j\), \(2\)]\) Subscript[f, 42][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[u, j] Subscript[u, k] 
\!\(\*SubsuperscriptBox[\(w\), \(j\), \(2\)]\) Subscript[f, 43][Subscript[b, j],Subscript[b, k]]+\[Epsilon] 
\!\(\*SubsuperscriptBox[\(u\), \(k\), \(2\)]\) 
\!\(\*SubsuperscriptBox[\(w\), \(j\), \(2\)]\) Subscript[f, 44][Subscript[b, j],Subscript[b, k]]+\[Epsilon] 
\!\(\*SubsuperscriptBox[\(u\), \(j\), \(2\)]\) Subscript[w, j] Subscript[w, k] Subscript[f, 45][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[u, j] Subscript[u, k] Subscript[w, j] Subscript[w, k] Subscript[f, 46][Subscript[b, j],Subscript[b, k]]+\[Epsilon] 
\!\(\*SubsuperscriptBox[\(u\), \(k\), \(2\)]\) Subscript[w, j] Subscript[w, k] Subscript[f, 47][Subscript[b, j],Subscript[b, k]]+\[Epsilon] 
\!\(\*SubsuperscriptBox[\(u\), \(j\), \(2\)]\) 
\!\(\*SubsuperscriptBox[\(w\), \(k\), \(2\)]\) Subscript[f, 48][Subscript[b, j],Subscript[b, k]]+\[Epsilon] Subscript[u, j] Subscript[u, k] 
\!\(\*SubsuperscriptBox[\(w\), \(k\), \(2\)]\) Subscript[f, 49][Subscript[b, j],Subscript[b, k]]+\[Epsilon] 
\!\(\*SubsuperscriptBox[\(u\), \(k\), \(2\)]\) 
\!\(\*SubsuperscriptBox[\(w\), \(k\), \(2\)]\) Subscript[f, 50][Subscript[b, j],Subscript[b, k]]


R[j_,k_][\[ScriptCapitalE]_] :=  \[ScriptCapitalE] // Subscript[R, 0][j,k] // (#+b[\[Rho]0[j,k],#])&


R3[\[ScriptCapitalE]_] := Expand[(\[ScriptCapitalE]//R[1,2] // R[1,3] // R[2,3])-(\[ScriptCapitalE]//R[2,3] // R[1,3] // R[1,2])]



