20 4 3 4
221 -339 138
-15   24 -10
-10   15  -6
 -6    8  -3
 -3    3  -1

first line: n (# responses), d (# c values), p (# y values), M (# solutions)

second line: y = yTbot

rest: R = HRbot, d-by-p

Example from CORVIDS paper:

% python3 -i solve.py
nm = 37 , m = 1.85 , vnn = 5820 , v = 0.765789508836 , s = 0.875094
[6, 13, 0, 0, 1, 0, 0]
[9, 5, 6, 0, 0, 0, 0]
[8, 8, 3, 1, 0, 0, 0]
[7, 11, 0, 2, 0, 0, 0]
search_count = 10312 , solution_count = 4
>>> cTlist
[[6, 13, 0, 0], [9, 5, 6, 0], [8, 8, 3, 1], [7, 11, 0, 2]]
>>> yTbot
Matrix([[221, -339, 138]])
>>> HRbot
Matrix([
[-15, 24, -10],
[-10, 15,  -6],
[ -6,  8,  -3],
[ -3,  3,  -1]])
>>> avg(cTlist)
[7.5, 9.25, 2.25, 0.75]
>>> avg(xTlist)
[7.5, 9.25, 2.25, 0.75, 0.25, 0.0, 0.0]
>>> 
