Add a function that fits a graph to a covariance matrix

[mailund]

Given a covariance matrix, rooted in an arbitrary leaf, as in Felsenstein chapter 23, fit the edge lengths to it. This just means fitting F2 statistics, essentially, and shouldn't be much different from what we do now with F4 and F3 statistics.