(* ::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\\WKO4"];
<< FreeLie.m;
<< AwCalculus.m;
Rs[a_,b_]:= Es[\[LeftAngleBracket]a->LS[0], b->LS[LW@a]\[RightAngleBracket], CWS[0]];
\[Alpha]=LS[{x,y}, \[Alpha]s]; \[Beta]=LS[{x,y}, \[Beta]s];\[Gamma]=CWS[{x,y}, \[Gamma]s];
V=Es[\[LeftAngleBracket]x->\[Alpha], y->\[Beta]\[RightAngleBracket], \[Gamma]];
\[Kappa]=CWS[{x}, \[Kappa]s]; Cap=Es[\[LeftAngleBracket]x->LS[0]\[RightAngleBracket],\[Kappa]];
R4Eqn=V**(Rs[x,z] // d\[CapitalDelta][x,x,y]) \[Congruent] Rs[y,z]**Rs[x,z]**V ;UnitarityEqn=(V**(V//dA[x]//dA[y]) \[Congruent] Es[\[LeftAngleBracket]x->LS[0], y->LS[0]\[RightAngleBracket],CWS[0]]);
CapEqn=((V**(Cap//d\[CapitalDelta][x,x,y])//dc[x]//dc[y]) \[Congruent] (Cap(Cap//d\[Sigma][x,y])//dc[x]//dc[y]));
TwistEqn = (V\[Congruent]Rs[x,y]**(V//d\[Sigma][{x,y}->{y,x}])**(El[\[LeftAngleBracket]x->LS[-LW@y/2], y->LS[-LW@x/2]\[RightAngleBracket], CWS[0]] // \[CapitalGamma]));
\[Alpha]s[x]= 0; \[Alpha]s[y]=-1/2;
SeriesSolve[{\[Alpha], \[Beta],\[Gamma], \[Kappa]}, (\[HBar]^-1 R4Eqn) && UnitarityEqn && CapEqn && TwistEqn];
{V, \[Kappa]}



