CPPLapack Code

#include <cpplapack.h>
#include <ctime>

int main(int argc, char** argv)
{
  const size_t size(3);
  CPPL::dgematrix A(size,size);
  CPPL::dcovector y(size);

  srand(time(NULL));
  for(size_t i=0; i<size; i++){
    for(size_t j=0; j<size; j++){
      A(i,j) =double(rand())/double(RAND_MAX);
    }
    y(i) =double(rand())/double(RAND_MAX);
  }
  std::cout << "A=\n" << A << std::endl;
  std::cout << "y=\n" << y << std::endl;
  
  A.dgesv(y);
  std::cout << "x=\n" << y << std::endl;

  return 0;
}