# Gibbs Phenomenon for Fourier Integral Represention S:=(x,v)->(1/Pi)*(Si(v+v*x)+Si(v-v*x)); # make an animation of Gibbs overshoot plots[animate]( plot, [S(x,v), x=-1.5..1.5],v=20..40);