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