mat X(5, 5, fill::randu);
vec a = { 1,2,4 };
X.print(“X”);
mat z;
for(int i=0;i<a.size();i++)
{
for (int j = 0; j < X.n_cols; j++)
{
if(a.at(i)==j)
{
z=join_horiz(z, X.col(j));
z.print(“Z:”);
}
}
}
mat X(5, 5, fill::randu);
vec a = { 1,2,4 };
X.print(“X”);
mat z;
for(int i=0;i<a.size();i++)
{
for (int j = 0; j < X.n_cols; j++)
{
if(a.at(i)==j)
{
z=join_horiz(z, X.col(j));
z.print(“Z:”);
}
}
}