
function [param] = boost(X,y,niter)

n = size(X,1);
p = ones(n,1)/n; 

param = [];
for i=1:niter,
  stump = find_stump(X,y,p);

  h = eval_stump(X,stump); % stump predictions on the training set

%  compute the number of votes to give to the new stump
%  alpha = 

%  update the weights
%  p = 

  param = [param,struct('stump',stump,'alpha',alpha)];

end;

