vscodeを用いたプログラム作成依頼

分子動力学でレナードジョーンズ法12-6を用いて、分子の動きを求めるプログラムを作成しています。 分子数、初期位置(x,y)、初速度(vx,vy)、時間刻みをファイルから読み込むことで、位置、速度、力を計算するプログラムをつくりました。 次に、これにカットオフ距離と周期境界条件を周辺8領域に適応したプログラムを作りたいのですが、作り方が分かりません。 どなたか作っていただけないでしょうか？ 今できているプログラムを添付しております。 このプログラムは、カットオフ距離rcを考慮した力の式に書き換えていますが、rcよりもrが大きい場合、ポテンシャルを0と見なす、というような部分までは出来ていません。 よろしくお願い致します。