fn cubic_fit( th0: f64, th1: f64, area: f64, mx: f64, ) -> ArrayVec<(CubicBez, f64, f64), 4>
Returns curves matching area and moment, given unit chord.