Вы можете использовать решение линейных наименьших квадратов в закрытой форме.
C = inv (transpose (X) * X) * transpose (X) * y
В приведенном выше примере make первая строка X все единицы, а следующие строки x1, x2, ...
C будут содержать соответствующие константы. Первая запись в C - c.