C++ Eigen 线性回归（最小二乘法、牛顿法）

代码：

#include "stdafx.h"
#include <iostream>
#include <Dense>
using namespace std;
using namespace Eigen;

MatrixXd newton(MatrixXd, MatrixXd, int, float, float);
MatrixXd least_squre(MatrixXd, MatrixXd);
int get_min_m(MatrixXd, MatrixXd, float, float, MatrixXd, MatrixXd, MatrixXd);
MatrixXd get_error(Matrix<double, -1, -1>, Matrix<double, -1, -1>, MatrixXd);
MatrixXd first_derivative(MatrixXd, MatrixXd, MatrixXd);
MatrixXd second_derivative(MatrixXd);
int main()
{
   
	MatrixXd feature(5,3);
	feature <<
		1, 1, 1,
		1, 2, 3,
		1, 4, 1,
		1, 5, 6,
		1, 3,