Skip to content
Snippets Groups Projects
Commit d3718d25 authored by Luke Naylor's avatar Luke Naylor
Browse files

Reorg some Sage code

parent c9f611fb
No related branches found
No related tags found
No related merge requests found
......@@ -211,9 +211,13 @@ This can be rearranged to express a bound on $d$ as follows:
\begin{sagesilent}
var("r_alt",domain="real") # r_alt = r - R/2 temporary substitution
bgmlv1_with_q_reparam = (bgmlv1_with_q.subs(r == r_alt + R/2)/r_alt).expand()
bgmlv1_d_ineq = ((0 > -bgmlv1_with_q_reparam) + 4*d)/4 # Rearrange for d
bgmlv1_d_ineq = bgmlv1_d_ineq.subs(r_alt == r - R/2).expand() # Resubstitute r back in
bgmlv1_d_ineq = (
((0 >= -bgmlv1_with_q_reparam)/4 + d) # Rearrange for d
.subs(r_alt == r - R/2) # Resubstitute r back in
.expand()
)
bgmlv1_d_lowerbound = bgmlv1_d_ineq.rhs() # Keep hold of lower bound for d
\end{sagesilent}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment