#fix dt/reset のXmax決める用 temp = 0 timing = 0 for i in range(M-1): if abs(x_detec1[i+1] - x_detec1[i])>temp: temp = abs(x_detec1[i+1] - x_detec[i]) timing = i print(temp) print(timing)

2) 10K以上の弾性波伝播で，detectorが波を受けるまでの間にゆっくりと座標を変化させている（右にずれている）。弾性波を流さずに長いタイムステップnveで放置したらどうなるか見て，この現象についても考える。10Knve_doNothing(nve_relax)など →前も書い…