function y = circularCdfWrappedNormal(th, mu, rho, varargin)
%
%< circularCdfWrappedNormal >
%  Returns cumulative distribution function value of Wrapped Normal distribution.
%
%  circularCdfWrappedNormal(th, mu, rho)
%       th:  variable
%       mu:  mean direction
%       rho: mean resultant length
%  Optional parameter
%       tol: tolerance level
%
%  Unit is radian.
%

warning('circularCdfWrappedNormal has not been implemented...')
y=NaN;
return;