function yy=cspline(x,y,xx)
       yy=spline(x,abs(y),xx).* ...
                      exp(i*spline(x,unwrap(phase(y)),xx));
return;