第一課題

　第二課題で使用するプログラム2D．fの２つのサブルーティンが未完成なので自分で完成させなさい。"subroutine M_calculation"は全磁気モーメントの磁気モーメントを足し合わせた磁化

を求めている。"subroutine E_calculation"は隣接した磁気モーメントと外部磁場による内部エネルギー

を求めている。但し交換相互作用定数J1、磁気モーメントspin(i,j),外部磁場H,plus1(i)はプログラム内で定義されておりplus1(i)=i+1でありi=Lの時は境界値条件よりplus1(L)=1である。