gibbs-phenomenon code snippets in matlab