  Macroscopically populated quantum superpositions pose a question to what
extent macroscopic world obeys quantum mechanical laws. Recently such
superpositions for light, generated by optimal quantum cloner, were
demonstrated. They are of fundamental and technological interest. We present
numerical methods useful for modeling of these states. Their properties are
governed by Gaussian hypergeometric function, which cannot be reduced to
neither elementary nor easily tractable functions. We discuss the method of
efficient computation of this function for half integer parameters and moderate
value of its argument. We show how to dynamically estimate a cutoff for
infinite sums involving this function performed over its parameters. Our
algorithm exceeds double precision and is parallelizable. Depending on the
experimental parameters it chooses one of the several ways of summation to
achieve the best efficiency. Methods presented here can be adjusted for
analysis of similar experimental schemes.
