#include "libscl.h" 
using namespace scl;

int main(int argc, char** argp, char** envp)
{

  INTEGER n=100;
  INTEGER p=4;
  INT_32BIT seed = 100542;

  realmat data(n,p);
  
  for (INTEGER i=1; i<=n; i++) {
    data(i,2) = 1.0;
    data(i,3) = ran(seed);
    data(i,4) = unsk(seed);
    data(i,1) = data(i,2)+data(i,3)+data(i,4)+0.5*unsk(seed); 
  }

  vecwrite("regr.dat",data);

  return 0;
}