function out = vysu2(y)

out = sqrt(1-0.25*y.^2);