Refactor supplementary conditions in Brans-Dicke Jordan frame

1 files changed, 4 insertions, 2 deletions

diff --git a/03-brans-dicke-jor.mac b/03-brans-dicke-jor.mac
index a5f0a74..3ea8aa5 100644
--- a/03-brans-dicke-jor.mac
+++ b/03-brans-dicke-jor.mac
@@ -179,14 +179,16 @@ e_sc_eqs[1] : expand_equation(expand(λ^2*sc_eqs[1]), coeffs)$
 eq_mono : block(
   [eq : coeff(e_sc_eqs[1], r, -2)=0],
   eq : solve(eq, diff(M,u))[1],
-  trigreduce(expand(eq)));
+  eq : subst(diff(λu_coeffs[2],u), eq),
+  trigreduce(expand(trigsimp(eq))));
 
 /* Supplementary Condition 2 */
-e_sc_eqs[2] : expand_equation(sc_eqs[2], coeffs);
+e_sc_eqs[2] : expand_equation(expand(λ^2*exp(2*β)*sc_eqs[2]), coeffs)$
 
 eq_dipo : block(
   [eq : coeff(e_sc_eqs[2], r, -2)=0],
   eq : -3*solve(eq, diff(N,u))[1],
+  eq : subst(diff(λu_coeffs[2],θ), eq),
   trigreduce(expand(eq)));