  We present a Monte Carlo study of the helicity moduli of an anisotropic
classical three-dimensional (3D) XY-model of YBCO in superconducting state. It
is found that both the ab-plane and the c-axis helicity moduli, which are
proportional to the inverse square of the corresponding magnetic field
penetration depth, vary linearly with temperature at low temperatures. The
result for the c-axis helicity modulus is in disagreement with the experiments
on high quality samples of YBCO. Thus we conclude that purely classical phase
fluctuations of the superconducting order parameter cannot account for the
observed c-axis electrodynamics of YBCO.
