Hartree Fock (HF) theory is one of the basic theories underlying the current understanding of the electronic structure of materials. It is a simple non-relativistic treatment of many electron system that accounts for the antisymmetric (fermion) nature of electronic wavefunction but does not account for correlation effects. The matrix representation of this theory constitutes the Roothan Equations. This simple MATLAB Code uses the Roothan Equations, within Born Oppenheimer approximation and non-relativistic limit, to calculate the energy of H-H (two Hydrogen atoms) system with different bond lengths using Slater Type Orbital with 4 Gaussian basis functions (STO-4G) for each spatial orbital. It can be used as an example for a class on electronic structure theory or as a starting skeleton for developing more complex ab initio codes, that use a similar basis set.