#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <R.h>
#include <Rmath.h>

void irls(double *xmatrix, double *yvec, int *p, int *nits, double *beta)
{
  int i,j;
  int niters=*nits;
  double *mu=0,t1=0,c1=0;
  double Xm=*xmatrix;
  double Wm;
  
  int pval=*p;
  for (i=0;i<pval;i++){
	  t1+=yvec[i];
    c1++;
  }
  for(j=0;j<pval;j++){
    mu[j]=c1;
  }
  for (i=0;i<niters;i++) {
	  for(j=0;j<pval;j++){
      beta[j]=mu[j];
    }
  }
}