function hybridhw;

% INPUT
% s      .. number of seats
% v      .. vector with the votes for each party
% 
% OUTPUT
% e      .. array with the number of seats for various methods

% inputs
s=input('give the number of seats  : ');
if s-fix(s)>0,error('---------->Give integer number of seats');end;
len=input('give the number of parties: ');
fprintf('     Put the parties in order of priority for ties.\n');
v=zeros(len,1);
for j1=1:len,
    fprintf('give votes for party - %1.0f  : ',j1);
    v(j1)=input('');                                   
end;
v=v.*(1+(len*1e-10:-1e-10:1e-10)');                                   % disturbance for handling ties
suma=sum(v);

p=s*v/suma;   % fair share

              % Hamilton
resids=sort(p-fix(p),'descend');tresh=resids(s-sum(fix(p)))-100*eps;
eh=fix(p)+(p-fix(p)>tresh);

              % Hybrid Hamilton/Webster)
resids=sort(p./(ceil(p)-1/2),'ascend');tresh=resids(sum(ceil(p))-s)+100*eps;
em=ceil(p)-(p./(ceil(p)-1/2)<tresh);

e= [apport(p,s,1:s+1) ...                                             % Jefferson | D'Hond 
    apport(p,s,1/2:s+1/2) ...                                         % Webster | Sainte-Lague
    eh ...                                                            % Hamilton
    em ...                                                            % Hybrid
    apport(p,s,sqrt((100*eps:s+100*eps).*(1+100*eps:s+1+100*eps))) ...% Huntington-Hill
    apport(p,s,(100*eps:s+100*eps).*(1:s+1)./(.5:s+.5)) ...           % Dean
    apport(p,s,100*eps:s+100*eps)];                                   % Adams    

fprintf('\n');
fprintf('   votes     %%   JEF  WEB  HAM  H-W  H/H  DEA  ADA \n');
for j1=1:len,
    fprintf('%8.0f %5.2f %4.0g %4.0g %4.0g %4.0g %4.0g %4.0g %4.0g\n',v(j1),100*v(j1)/suma,e(j1,:))
end;
return

%------------------------------------------------------------------
function e=apport(p,s,diairetis);
aver=kron(p,1./diairetis)+1e-11*randn(length(p),s+1);
divid=sort(reshape(aver,(s+1)*length(p),1),'descend');
metro=divid(s)-1e-12;
e=sum(aver'>=metro)';
return