楕円 の速度ベクトル、加速度ベクトルのアニメーションを作る。可視化に必要な関数を下記で定義する。 sage: t = var('t') sage: x = 9*cos(t) sage: y = 4*sin(t) sage: vx,vy = x.derivative(),y.derivative() sage: ax,ay = vx.derivative(),vy.derivativ…

はじめに 前回のつづきです。今回も数値計算っぽいことをやります。小出昭一郎『量子力学』のp.53より の期待値はとなります。前回と同様 とし、これを数値計算で確かめます。 方法 から をサンプリングして、 の平均値をとれば が計算できます。 からサンプ…