From: Ian Jackson <ijackson@chiark.greenend.org.uk>
Date: Mon, 13 Nov 2017 22:56:34 +0000 (+0000)
Subject: symbolic.py: wip p_nosing
X-Git-Url: http://www.chiark.greenend.org.uk/ucgi/~ian/git?a=commitdiff_plain;h=c8687eaa179f080824096ac5196803bbb998134f;p=moebius3.git

symbolic.py: wip p_nosing

Signed-off-by: Ian Jackson <ijackson@chiark.greenend.org.uk>
---

diff --git a/symbolic.py b/symbolic.py
index c3d89db..c70cda6 100755
--- a/symbolic.py
+++ b/symbolic.py
@@ -21,6 +21,24 @@ p_start = Matrix([
   mu * s,
 ])
 
-p_rightvars = p_start.subs(theta, s/r).subs(r, 1/la)
+p_rightvars = p_start.subs( theta, s/r ).subs( r, 1/la )
 
 pprint(p_rightvars)
+
+dirn_rightvars = diff(p_rightvars, s)
+
+pprint(dirn_rightvars)
+
+zeta = Wild('zeta')
+
+p_nosing = (p_rightvars
+            .replace( 1-cos(zeta)  ,   2*sin(zeta/2)**2          )
+            .replace( sin(zeta)**2 ,   zeta*sinc(zeta)*sin(zeta) )
+            )
+p_nosing[1] = (p_nosing[1]
+            .replace( sin(zeta) , zeta * sinc(zeta)                )
+               )
+
+print
+pprint(p_nosing,)
+