状態密度の計算

各原子軌道(2s,2p_x,2p_y,2p_x)ごとに状態密度を求め，

その値を積分のプログラムを用いて幅を持たせてグラフにする．

計算方法

原子軌道iの状態密度 (Σ(S_j,i x C_k,j))/N

合計は，ss,sp,ppそれぞれ分けてやる．

(S:重なり行列 C:固有ベクトル)

(N:原子数 k:エネルギー固有値に対応する番号 i,j:原子軌道の番号)

計算プログラム