void headerml(int niter, int nburn, int nuse)
{
   printf("Marginal likelihood analysis \n");
   printf("Total iterations avaliable:%10d\n",niter);
   printf("        Burn-in iterations:%10d\n",nburn);
   printf("       Analysis iterations:%10d\n",nuse);
}