Numerical Solution of Symmetric Positive Differential Equations