example.cpp

#include <Pds/Ra>
 
int main(void)
{
    Pds::Matrix A = "1 -2  2\n"
                    "2  1  1\n"
                    "1  1 -1";
    
    Pds::Matrix Y = "1\n"
                    "2\n"
                    "3";
    
    // Given A*X=Y, then
    Pds::Matrix X = A.Inv()*Y;
    X.Print("Data in matrix X:\n");
    
    (A*X).Print("Data A*X:\n");
    
    Pds::Ra::PrintVersion("Version:\t");
    Pds::Ra::PrintPackage("Package:\t");
    
    return 0;
}