MatheMatica-Solutions-for-Hubbard-model This project include: Estimation biquadratic and bilinear exchange with multi-orbitals Hubbard model by exact diagonalization method CG-method Slate-Koster-method