power-system code snippets in matlab