def IDST(X): """ Inverse DST. Python -> Matlab """ n = X.shape[0] x = np.real(scipy.fftpack.idst(X,type=1,axis=0)) return x/(n+1.0)