function xpadpad=myinterp(x);
if ~rem(length(x),2) % if not odd, make it odd
  xpad=[x(:);0];
else%if
  xpad=x(:);
end%if
% now xpad is odd
XPAD=fft(xpad);
l=length(XPAD);
PAD=zeros(l,1);
XPADPAD=[XPAD(1:(l+1)/2);PAD;XPAD((l+1)/2+1:l)];
xpadpad=ifft(XPADPAD);
%%%if isreal(x)

