计算C++中矩阵的伪逆

#include<iostream>
#include <Eigen/Dense>

using namespace std;
using namespace Eigen;

int main()
{
    MatrixXd A(3, 2);
    A << 1, 2,
         3, 4,
         5, 6;

    MatrixXd pinvA = A.completeOrthogonalDecomposition().pseudoInverse();

    cout << "The pseudo-inverse of A is:\n" << pinvA<< endl;

    return 0;
}