function f=nfrankx(m,k,h,n);

% Programme auxiliaire dans l'implémentation 
% du programme nfrankl. 


j=2^m*3^k*5^h;
   C=zeros(j+1,n*j+1);
for k=0:n
    C(1,1+k)=(n-k)/n;
end;
y=size(C);
for r=2:(y(1)-1)
  d=1+n*(r-2);
  for kc=0:n
      i=n-kc;
      C(r,d+n-i)=(n-i)/n;
      C(r,d+2*n-(n-i))=(n-i)/n;
  end;
end;

for i=0:n
    C(y(1),n*j+1-i)=(n-i)/n;
end;
f=C;