if(exist("t")==0 || t<0) t = tmin # tが定義されていない場合やt<0の場合はtにtminを代入 filename=sprintf("wavefunction_%07d.dat", t) # データファイルの名前を設定 timelabel=sprintf("t = %5d*10^{-4} (hbar/MeV)", 100*t) # 図中に表示するラベル(時刻情報)を設定(なくてもよい) unset label # 一つ前の時刻のときのラベルを削除 set label timelabel at 10, 0.4 # 座標(10,0.4)にラベルを配置 # データのプロット。実部、虚部と絶対値をそれぞれプロット plot \ filename using 1:2 t "Re psi(x,t)" w l lt 1 lw 1,\ filename using 1:3 t "Im psi(x,t)" w l lt 2 lw 1,\ filename using 1:(($2)**2+($3)**2)**0.5 t "|psi(x,t)|" w l lt 3 lw 3 # tをdtだけすすめる t = t + dt # tがtmaxを超えない場合はこのファイルの冒頭に戻って動作を繰り返す if ( t