function [C]=FamilySize2(J,MAX)
% clear all;close all;[C]=FamilySize2(10^6,25)
D=2509;
q=[125*231/(D^2) 1/3 34/D ];q=q/2;X=[1 2 1 0];
X=[1 2 1 0];
C=zeros(MAX,1);
for N=1:MAX
    N
    C(N)=Multinomial(N,q,X)*(J/N);
end;
semilogy(C,'r')


   