function fraction(d)
   frac (1,2,d)
   fprintf('\n')
end

function frac(n,d,m)
   if d <= m
      frac(n*2-1, d*2, m)
      fprintf('%d/%d ', n, d)
      frac(n*2+1, d*2, m)
   end
end

