function mdot = xcf(t,m);
% XCF solves the 2 x 2 system of ODEs for projectivized
% structure constants of 3d hom sp under the XCF
% Call syntax: mdot = xcf(t,m);
% Inputs are t = time, m = [m(1),m(2)]
% Output is  
%    mdot =
%    -[m(1)*(-3*m(2)^2+1   +m(1)^2 -2*m(1) +2*m(1)*m(2) +2*m(2))*(1-m(1))*(1 + m(1) - m(2))
         m(2)*(m(2)^2   +1 -3*m(1)^2 +2*m(1) +2*m(1)*m(2) -2*m(2))*(1-m(2))*(1 + m(2) - m(1))]
%---------------------------------------------------------
mdot = -[m(1)*(-3*m(2)^2+1   +m(1)^2 -2*m(1) +2*m(1)*m(2) +2*m(2))*(1-m(1))*(1 + m(1) - m(2)); ...
         m(2)*(m(2)^2   +1 -3*m(1)^2 +2*m(1) +2*m(1)*m(2) -2*m(2))*(1-m(2))*(1 + m(2) - m(1))];