symbolic.py: wip

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

diff --git a/symbolic.py b/symbolic.py
index 9de7ccd136b559b804ad723923a890cba977ac1c..42a24e3734f23a7a81552ae056aa8d8cacb0b9d8 100755 (executable)
--- a/symbolic.py
+++ b/symbolic.py
@@ -26,10 +26,10 @@ p_rightvars = p_start.subs( theta, s/r ).subs( r, 1/la )
 print
 pprint(p_rightvars)
 
-dirn_rightvars = diff(p_rightvars, s)
+p_dirn_rightvars = diff(p_rightvars, s)
 
 print
-pprint(dirn_rightvars)
+pprint(p_dirn_rightvars)
 
 zeta = Wild('zeta')
 
@@ -44,3 +44,16 @@ p_nosing[1] = (p_nosing[1]
 print
 pprint(p_nosing,)
 
+t = symbols('t')
+
+q_owncoords = p_nosing.replace(s,t).replace(la,-la)
+q_dirn_owncoords = p_dirn_rightvars.replace(s,t).replace(la,-la)
+
+q_dirn_owncoords_0 = q_dirn_owncoords.replace(t,0)
+
+print
+pprint(q_owncoords)
+print
+pprint(q_dirn_owncoords)
+print
+pprint(q_dirn_owncoords_0)