SSL Minutes 2003-10-16 (Thursday)
Note-taker: Hal
-- Minutiae: --
Karl will take notes next time.
Martin will provide refreshments on Tuesday.
We will continue to meet at 2:30 on Tuesdays.
-- Old Business: --
Getting graph.tcl (http://math.wisc.edu/~propp/software.html)
to work. Hal suggests having XWin-32 installed in B107 by the
administrator.
Jim will talk to Yvonne Nagel about this and the computational
cluster (see his email for follow-up).
-- Reports: --
* Martin is working on (or thinking about) a program that will
find a generating function from a matrix. He wants to apply
it to 3-D snakes. Will try and finish it by next week.
* Hal and Karl have been looking at the Apollonian Gasket
(because it is so cool.) We were trying to discover the
Soddy relationship on our own, but were unsuccessful. So we
googled to
http://mathworld.wolfram.com/SoddyCircles.html
here we found the relationship given by this maple function
(Hal got garbage when he tried this the first time because
he left out a * in the formula.)
To see some neat math, plug the following into maple:
w := (x,y,z) -> simplify((x*y*z)/( (x*y) + (y*z) + (x*z)
+ 2 * ((x*y*z*(x+y+z))^(1/2)) ));
f := (x,y,z) -> (w(x,y,z),y,z);
g := (x,y,z) -> (x,w(x,y,z),z);
h := (x,y,z) -> (x,y,w(x,y,z));
f(1,1,1);
g(f(1,1,1));
f(g(f(1,1,1)));
h(g(f(1,1,1)));
f(f(h(g(f(1,1,1)))));
Note: hal is going to put some of his notes at:
http://ups.physics.wisc.edu/~hal/SSL/2003/apollonian_gasket_maple.html
While still on the Apollonian Gasket, we talked about how
the equation has a plus-or-minus part. These correspond to
the inside and outside circles. Stephen asked if there was
a nice automorphism of R^2 that leaves the original triangle
the same and swaps out the inner and outer circle. This
leaves us the question: what does it mean to turn a circle
inside out? (negative curvature)
* Abby has been looking at the 491 notes.
(http://math.wisc.edu/~propp/491/) She also has still been
looking at 3-d snakes.
* Paul noticed something neat about the 3-d snakes. Let G be
the graph of a flat cubicular snake. Let E' be the set of
edges that connect the two layers. Let phi be a subset of
E'. Let M_phi be the number of matching on one layer of G
with the vertices removed that correspond to the edges in
phi. Then, it can be seen that:
M(G) = Sum(for each phi in E') of (M_phi)^2
Later on, this was generalized to any graph G. Let {A,B}
be a partition of the vertices of G. Define E' to be those
edges that join vertices in A to vertices in B. Choose A and
B such that no two elements in E' are incident to the same
vertex.
Then define M(A,phi) to be the number of matchings on A
after all of the vertices incident to phi have been removed.
Then
M(G) = Sum(for each phi in E') of ( M(A,phi) * M(B,phi) ).
Jim also talked about another way to define levels on a graph
such that each level contains only vertices from one of the
partite sets. Then you can always know how many edges in a
perfect matching will go from one level to the next.
* Emilie looked at the tree papers: Kuo, Itsara, and Somos.
* Jon also read through the articles.
* Stephen has been studying up on cluster algebras.
* Sam has been reading Kuo and Itsara.
* Jim has been playing with snakes. We saw in after.mws the
function -> simplify((y^2+z^2)/x).
Jim defined a sequence with this:
s_0 = x
s_1 = y
s_3 = z
s_n = (s_{n-1}^2 + s_{n-2}^2)/s_{n-3}
and then looked at the coefficients. He didn't find anything
good looking. So he tried:
s_0 = x
s_1 = y
s_3 = z
s_n = (a * s_{n-1}^2 + b * s_{n-2}^2)/s_{n-3}
where a and b are formal variables, just like x,y,and z.
Then he gets an interesting sequence for the biggest
coefficient in each polynomial.:
1, 2, 6, 60, 3052, 2519660, ...
Why?
Here are Jim's notes from the agenda:
-------------------------------------------------------
Maximal coefficients of the Scott-2 coefficients
(round, therefore combinatorial):
s := proc(n) option remember;
if n=0 then x
elif n=1 then 1
elif n=2 then z else
simplify( (a*s(n-1)^2+b*s(n-2)^2)/s(n-3) ); fi; end;
seq(max(coeffs(numer(s(n)))),n=3..8);
1, 2, 6, 60, 3052, 2519660
NOT IN THE HANDBOOK!
seq(ifactor(max(coeffs(numer(s(n))))),n=3..8);
1, (2), (2)*(3), (2)^4*(5), (2)^6*(3)*(11)*(13),
(2)^7*(7)^2*(13)^2*(17)*(19)*(23)
subs(z=0,numer(s(8)));
seq([coeffs(subs(z=0,numer(s(n))))],n=4..8);
seq([coeffs(subs(y=0,numer(s(n))))],n=4..8);
seq([coeffs(subs(x=0,numer(s(n))))],n=4..8);
seq(degree(subs(x=1,y=1,numer(s(n)))),n=3..8);
-------------------------------------------------------
Then we broke into groups and discussed the papers until we were out
of time.