function outdata=H_interp2(indata)
od=zeros(length(indata),1);
for i=1:13
    a=i/53;
    c1=a*(a-1)/2;
    c0=-(a-1)*(a+1);
    c_1=a*(a+1)/2;
    od(5+i)=c1*indata(4)+c0*indata(1)+c_1*indata(2);
    od(18+i)=c1*indata(1)+c0*indata(2)+c_1*indata(3);
    od(31+i)=c1*indata(2)+c0*indata(3)+c_1*indata(4);
end
for i=1:10
    a=i/53;
    c1=a*(a-1)/2;
    c0=-(a-1)*(a+1);
    c_1=a*(a+1)/2;
    tmp(i)=c1*indata(3)+c0*indata(4)+c_1*indata(1);
end
od(1:5)=tmp(6:10);od(45:49)=tmp(1:5);
outdata(1:24)=od(1:24);
outdata(25:48)=od(26:49);
outdata=conj(outdata);