Библиотека Armadillo для решения СЛАУ

ddp = solve(A, B);

************** Exception Text **************
System.Runtime.InteropServices.SEHException (0x80004005): External component has thrown an exception.
   at _CxxThrowException(Void* , _s__ThrowInfo* )
   at arma.?A0x3e4acfb2.arma_stop<char const [52]>($ArrayType$$$BY0DE@$$CBD* x) in d:\programs\visual_cpp\vc\include\armadillo\armadillo_bits\debug.hpp:line 107
   at arma.auxlib.solve<double,class arma::Mat<double> >(Mat<double>* out, Mat<double>* A, Base<double\,arma::Mat<double> >* X, Boolean slow) in d:\programs\visual_cpp\vc\include\armadillo\armadillo_bits\auxlib_meat.hpp:line 3378
   at arma.glue_solve.solve_direct<double,class arma::Mat<double> >(Mat<double>* out, Mat<double>* A, Base<double\,arma::Mat<double> >* X, Boolean slow) in d:\programs\visual_cpp\vc\include\armadillo\armadillo_bits\glue_solve_meat.hpp:line 28
   at arma.glue_solve.apply<class arma::Mat<double>,class arma::Mat<double> >(Mat<double>* out, Glue<arma::Mat<double>\,arma::Mat<double>\,arma::glue_solve>* X) in d:\programs\visual_cpp\vc\include\armadillo\armadillo_bits\glue_solve_meat.hpp:line 62
   at arma.Mat<double>.operator=<class arma::Mat<double>,class arma::Mat<double>,class arma::glue_solve>(Mat<double>* , Glue<arma::Mat<double>\,arma::Mat<double>\,arma::glue_solve>* X) in d:\programs\visual_cpp\vc\include\armadillo\armadillo_bits\mat_meat.hpp:line 4810
   at arma.solve<class arma::Mat<double>,class arma::Mat<double> >(Mat<double>* out, Base<double\,arma::Mat<double> >* A, Base<double\,arma::Mat<double> >* B, Boolean slow, Double* junk) in d:\programs\visual_cpp\vc\include\armadillo\armadillo_bits\fn_solve.hpp:line 124
   at F(fixed<16>* , fixed<16>* p) in d:\studying\c++\visual_studio\dynamics\double_rod_with_slider\opengl_on_a_windows_form\ode.h:line 97

#include <armadillo>
 
int main(){
 
    double A_[5][5] = { {1, 5, 5, 6, 1},
                        {6, 2, 6, 1, 12},
                        {1, 6, 1, 1, 13},
                        {6, 1, 1, 1, 1},
                        {5, 12,12,1, 5} };
    double B_[5] = {5, 1, 5, 12, 13};
 
    arma::mat::fixed<5,5> A;
    arma::vec::fixed<5> B;
 
    for(int i = 0 ; i < 5; i ++){
        for(int j = 0; j < 5; j++)
            A(i,j) = A_[i][j];
        B(i) = B_[i];
    }
    std::cout << solve(A,B);
    return 0;
}