(* ::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\\OneCo-1604";
<< Global.m;
Simp=Expand;
ExportButton


\[Rho]0[j_,k_]:=UU[\[Delta]a[Subscript[ff, 1][Subscript[b, j],Subscript[b, k]],\[CCedilla],j]+\[Delta]a[Subscript[ff, 2][Subscript[b, j],Subscript[b, k]],\[CCedilla],k]+\[Delta]a[Subscript[ff, 3][Subscript[b, j],Subscript[b, k]],j,j]+\[Delta]a[Subscript[ff, 4][Subscript[b, j],Subscript[b, k]],j,k]+\[Delta]a[Subscript[ff, 5][Subscript[b, j],Subscript[b, k]],k,j]+\[Delta]a[Subscript[ff, 6][Subscript[b, j],Subscript[b, k]],k,k]+\[Delta]aa[Subscript[ff, 7][Subscript[b, j],Subscript[b, k]],\[CCedilla],j,j,j]+\[Delta]aa[Subscript[ff, 8][Subscript[b, j],Subscript[b, k]],\[CCedilla],j,j,k]+\[Delta]aa[Subscript[ff, 9][Subscript[b, j],Subscript[b, k]],\[CCedilla],k,j,j]+\[Delta]aa[Subscript[ff, 10][Subscript[b, j],Subscript[b, k]],\[CCedilla],k,j,k]+\[Delta]aa[Subscript[ff, 11][Subscript[b, j],Subscript[b, k]],\[CCedilla],j,k,j]+\[Delta]aa[Subscript[ff, 12][Subscript[b, j],Subscript[b, k]],\[CCedilla],j,k,k]+\[Delta]aa[Subscript[ff, 13][Subscript[b, j],Subscript[b, k]],\[CCedilla],k,k,j]+\[Delta]aa[Subscript[ff, 14][Subscript[b, j],Subscript[b, k]],\[CCedilla],k,k,k]+\[Delta]aa[Subscript[ff, 15][Subscript[b, j],Subscript[b, k]],j,j,j,j]+\[Delta]aa[Subscript[ff, 16][Subscript[b, j],Subscript[b, k]],j,j,j,k]+\[Delta]aa[Subscript[ff, 17][Subscript[b, j],Subscript[b, k]],j,k,j,k]+\[Delta]aa[Subscript[ff, 18][Subscript[b, j],Subscript[b, k]],j,j,k,j]+\[Delta]aa[Subscript[ff, 19][Subscript[b, j],Subscript[b, k]],j,j,k,k]+\[Delta]aa[Subscript[ff, 20][Subscript[b, j],Subscript[b, k]],j,k,k,k]+\[Delta]aa[Subscript[ff, 21][Subscript[b, j],Subscript[b, k]],k,j,k,j]+\[Delta]aa[Subscript[ff, 22][Subscript[b, j],Subscript[b, k]],k,j,k,k]+\[Delta]aa[Subscript[ff, 23][Subscript[b, j],Subscript[b, k]],k,k,k,k]]


R[j_,k_] := R[j,k]=Ea[j,k]**TSD[{
j->UU[a[1,j,h\[Infinity]]]+bb[j,k][\[Rho]0[j,k],UU[a[1,j,h\[Infinity]]]],
k->UU[a[1,k,h\[Infinity]]]+bb[j,k][\[Rho]0[j,k],UU[a[1,k,h\[Infinity]]]]
}];


lhs=R[1,2]**R[1,3]**R[2,3];
rhs=R[2,3]**R[1,3]**R[1,2];


Simp=Simplify


Subscript[ff, 21][__]=0;


Subscript[ff, 5][__]=0;


Subscript[ff, 11][__]=0;


Subscript[ff, 22][__]=0;


Subscript[ff, 18][__]=0;


Subscript[ff, 13][__]=0;


Subscript[ff, 8][x_,y_]:=Subscript[gg, 1][x];


Subscript[ff, 7][__]=0;


Subscript[ff, 1][x_,y_]:=Subscript[gg, 2][y];


Subscript[ff, 12][x_,y_]:=Subscript[gg, 3][y];


Subscript[ff, 10][x_,y_] := -x Subscript[ff, 17][x,y]


Subscript[gg, 1][x_]=Subscript[cc, 1]/x;


Subscript[ff, 17][x_,y_] := -((E^-x (2+E^x (-2+x)+x+(-1+E^x) x Subscript[cc, 1]))/(2 x^3));


Subscript[ff, 14][x_,y_] := -x Subscript[ff, 20][x,y]


Subscript[ff, 16][x_,y_]:=Subscript[gg, 4][x];


Subscript[ff, 9][x_,y_]:=Subscript[gg, 5][x];


Subscript[ff, 20][x_,y_]:=(Subscript[gg, 6][y](E^x-1))/(x E^x);


Subscript[ff, 2][x_,y_]:=Subscript[gg, 7][x]-x Subscript[ff, 4][x,y]


Subscript[gg, 3][y_] := (2-2 E^y+y+E^y y-(1+E^y) y Subscript[cc, 1]-2 E^y y^2 Subscript[gg, 5][y]-2 y^2 Subscript[gg, 6][y]+2 E^y y^2 Subscript[gg, 6][y])/(2 y^2);


Subscript[ff, 4][x_,z_]:=Subscript[gg, 8][x]+(E^x-1)/(x E^x) Subscript[gg, 9][z];


Subscript[gg, 9][y_] := 1/(2 (-1+E^y) y) (2-2 E^y+y+E^y y-y Subscript[cc, 1]+E^y y Subscript[cc, 1]+2 y Subscript[gg, 2][y]+2 y^2 Subscript[gg, 6][y]-2 E^y y^2 Subscript[gg, 6][y]+2 E^y y Subscript[gg, 7][y]-2 E^y y^2 Subscript[gg, 8][y]);



