function [Cxyo,Fc]=Ogive(Cxy,fs)
% Computes Ogive cumulative function

[Cc,Fc,dfc]=specsmooludot(Cxy,fs);
wc0=real(nansum(Cc.*dfc));
Cxyo=cumsum(Cc.*dfc)./wc0;