Pensieve Header: Comparing the perturbative and the exact solutions.

SetDirectory["C:\\drorbn\\AcademicPensieve\\2012-05\\beta5.1"];
<< betaCalculus.m
Clear[ℏ]; Unprotect[C];
$PerturbativeDegree=8;
PExpand[expr_]:=expr /. sd_SeriesData :> MapAt[Expand, sd, 3];

{V, C, sol}=Get["ExactSolution-120528.m"]

V1=PExpand[Series[# /. {(Subscript[c, 1])->ℏ (Subscript[c, 1]), (Subscript[c, 2]) -> ℏ (Subscript[c, 2])}, {ℏ,0,8}]& /@ V]

C1 = PExpand[Series[# /. {(Subscript[c, 1])->ℏ (Subscript[c, 1]), (Subscript[c, 2]) -> ℏ (Subscript[c, 2])}, {ℏ,0,8}]& /@ C]

{V2,C2,sol2} = Get["SolutionToDegree8-120524.m"]

V1==V2

C1==(C2 /. (Subscript[κ, 1])->0)