#include <D:\devcpp\all.h>

int main()
{
  long i;
  srand((unsigned)time(NULL));
  Matrix A(8,3);
  A.print(0);

  A.bin_ran(50,0.5).print(0);
  cout<<A.col()<<endl;
  cout<<A.row()<<endl;
  A.col(1).print(0);
  A.row(0).print(0);
  cout<<"mean of A="<<endl;
  A.mean().print(3);
  A.T().mean().print(3);

  Matrix B(3,5);
  B.bin_ran(30,0.6).print(0);
  cout<<"A*B="<<endl;
  (A*B).print(3);

  /*Univat X(100);
  X.print(2);
  X.bin_ran(100,0.3).print(0);
  X.mean().print(3);
  X.variance().print(3);
  cout<<X.size()<<endl;
  Univar Y(100);
  Y.pos_ran(1.0).print(0);
  corr(X,Y).print(3);//相關係數
  */
  cin.get();
  return 0;
}