原島鮮『初等量子力学』10.2 より
連続関数がデルタ関数へ接近していく様子を可視化します。
ツールは SageShell を使います。
その1
sage: l = [plot(sin(a*x)/pi/x,(x,-15,15))+text("$a={}$".format(a),(15*0.7,a/pi*0.8),fontsize=20, bounding_box={'boxstyle':'round', 'fc':'w'})+text("$y=\sin (ax)/\pi x$",(-15*0.7,a/pi*0.8),fontsize=20,color="black") for a in range(1,31)] sage: animate(l) Launched gif viewer for Animation with 30 frames
その2
sage: l = [plot(a/(a^2+x^2)/pi,(x,-15,15))+text("$a={}$".format(a),(15*0.7,1/a/pi*0.8),fontsize=20,bounding_box={'boxstyle':'round', 'fc':'w'})+text("$y=a/\pi (a^2+x^2)$",(-15*0.7,1/a/pi*0.8),fontsize=20,color="black") for a in [1.1^(-i) for i in range(20)]] sage: animate(l) Launched gif viewer for Animation with 20 frames