The Gauss-Seidel program is kept in the file loadflow_gs .m . This calls the ybus .m function .