(* ::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:: *)
MyCollect[\[ScriptCapitalE]_,vs_List]:=MyCollect[\[ScriptCapitalE],vs,Identity];
MyCollect[\[ScriptCapitalE]_,vs_List,simp_]:=Total[CoefficientRules[\[ScriptCapitalE],vs]/.((ps_->c_):>simp[c]Times@@(vs^ps))];
MyCollect[\[ScriptCapitalE]s_List,vs_List] := MyCollect[#,vs]&/@\[ScriptCapitalE]s;
MyCollect[\[ScriptCapitalE]s_List,vs_List,simp_] := MyCollect[#,vs,simp]&/@\[ScriptCapitalE]s;
MyCollect[sd_SeriesData,vs_List] := MapAt[MyCollect[#,vs]&, sd, 3];
MyCollect[sd_SeriesData,vs_List,simp_] := MapAt[MyCollect[#,vs,simp]&, sd, 3];


(* ::Input::Initialization:: *)
r=\[Omega]^-1+((-2T \[Omega] d\[Omega])/(T-1) x y+2T \[Omega] d\[Omega] a +(\[Omega] T d\[Omega]-(T-1)^2/T p1)) \[Epsilon]/\[Omega]^3+(2 T \[Omega]^2 (2 d\[Omega]^2 T-d\[Omega] \[Omega]-dd\[Omega] T \[Omega])a^2+((2 (-1+T) p1 ((1+T) \[Omega]-3 (-1+T) T d\[Omega]))/T+2 \[Omega] ((-1+T)^2 dp1+2 T^2 (d\[Omega])^2-T \[Omega](d\[Omega]+T dd\[Omega])))\[Omega] a+(T \[Omega]^2 (4 d\[Omega]^2 (-1+T) T-d\[Omega] (-3+T) \[Omega]-2 dd\[Omega] (-1+T) T \[Omega]))/((-1+T)^3)  x^2 y^2-(2 \[Omega](-3 d\[Omega] p1 (-1+T) T+dp1 (-1+T) T \[Omega]+p1 (1+T) \[Omega]))/T  x y+(4 T \[Omega]^2 (2 d\[Omega]^2 (1-T) T-d\[Omega] \[Omega]-dd\[Omega] (1-T) T \[Omega]))/((-1+T)^2)  a x y-(p2/2+\[Omega] (3 d\[Omega] p1 (-1+T)^2-((p1 (-1+T^2)+T (dp1 (-1+T)^2+2 d\[Omega]^2 T^2)) \[Omega])/T+T (d\[Omega]+dd\[Omega] T) \[Omega]^2))) \[Epsilon]^2/\[Omega]^5



