%enter the following
out = 'mlout';     %output file (will contain error messages if any)
simfile= 'client1';  %client's simulation file (output)
hfile = 'sima';    %investigator's simulation file


%do not write anything below
LDPAR=250;
fout=fopen(out,'w');
if fout==-1; disp('failure to open file'); disp(out); quit; end;
simout=fopen(simfile,'w');
if simout==-1; fprintf(fout,'failure to open file %s', simfile); quit; end;
in=fopen(hfile);
if in==-1; fprintf(fout, 'failure to open file %s', hfile); quit; end;
niter=fscanf(in,'%d',1);   errcheck(niter,fout,0);
npar=fscanf(in,'%d',1);    errcheck(npar,fout,0);
if (niter <=0) | (npar<=0);
   fprintf(fout, 'illegal first record of simulation input file  %s %5d %5d', hfile, niter, npar); 
   quit; 
end;
if (npar>LDPAR);
   fprintf(fout, 'input simulation file  %s %5d parameters exceed limit of %5d',hfile, npar, LDPAR); 
   quit; 
end;

for iter=1:niter;
   [jter, wtlog, pri, pdat, thetai]=readrec(npar,in,fout);
   [nc,gc,priorc]=client(npar,thetai,pri);
   if iter==1;
      if (nc < 1) | (nc > LDPAR);
         fprintf(fout, 'Client`s function of interest %2d not in range 1 to %3d',nc,LDPAR);
      end;
      fprintf(simout, '%5d %16d\n', niter, nc); 
   end;
   wtlog=wtlog+priorc-pri;
   pri=priorc;
   writerec(simout,iter,wtlog,pri,pdat,gc); 
end; %for iter

fclose(in);
fclose(fout);
fclose(simout);
quit;