以下は、用C / C + +プログラムウィンドウおよびソースコードを計算し、 zbusは、与えられた電気ybusとのネットワークを任意のサイズです。は、入力システムが書かれ" rx.txt "は、次の列の形式; 'からバス' 、 'をバス' 、 '研究' 、および、 [ x ]をします。 'よりバス'と'をバス'しなければならない正の整数を連続でゼロからスタートします。 バスのゼロが想定されるスイングの緩みまたはバスです。 複雑なヘッダファイルを使用して表の操作を支援する。 zbusは、単純に計算さybus反転しています。

//-----------------------------
機能＃ pragma hdrstop
//----------------------------
機能＃ pragma argsused
＃含める < stdio.h >
＃含める< conio.h >
場合＃ include " matrix.h "
＃含める< complex.h >
を使用してのstd ： ：複雑な;
を使用して名前空間の数学;
定義typedef複雑な複雑な;
定義typedef表表;
無効メイン（無効）
（
/ /変数を宣言する
int CNTの、私は、無し、 matsize = 1 ;
ダブルtmp3 、 tmp4 ;
複雑なtmp6 ;
表 y 、 ybus 、 zbus ;
ファイル* 、 ;
y. 、 setSize （ matsize 、 matsize ） ; / /一時的な表のサイズを設定する
/ /データを入手して支店
で= fopen （ " rx.txt " 、 "研究" ） ; / /ファイルをオープン
もし（ ==ヌル）
（ printf （ " \ nrx.txtが見つかりませんでした" ） ;
getch （ ） ;
）
CNTの=無し=私= 0 ;
を（ ; ; ）
（もし（読み込むfscanf （には、 " ％ d個" 、 ＆私） == EOFの）
（休憩; ）
もし（私> matsize ）
（ matsize =私; ）
読み込むfscanf （には、 " ％ d個" 、 ＆無し） ;
もし（無し> matsize ）
（ matsize =無し; ）
y. 、 setSize （ matsize +1 、 matsize +1 ） ;
読み込むfscanf （には、 " ％のLF " 、 ＆ tmp3 ） ;
読み込むfscanf （には、 " ％のLF " 、 ＆ tmp4 ） ;
tmp6 =団地（ tmp3 、 tmp4 ） ;
Yに変更（私は、無し） = 1.0/tmp6 ;
CNTの+ + ;
）
fclose （ ） ;
/ /最後に就学許可証の支店のデータ
/ /作成ybus
ybus.setsize （ matsize +1 、 matsize +1 ） ;
/ /対角線のエントリ
を（私= 0 ;私< = matsize ;私+ + ）
（を（無し= 0 ;無し< = matsize ; J + +の）
（ ybus （私、私） = ybus （私、私） + Yに変更（私は、無し） ; ）
）
を（無し= 0 ;無し< = matsize ; J + +の）
（は（私= 0 ;私< = matsize ;私+ + ）
（ ybus （ジャック、無し） = ybus （ジャック、無し） + Yに変更（私は、無し） ; ）
）
/ /終了
/ /オフ対角線のエントリ
CNTの= 0 ;
を（無し= CNTの;無し< = matsize ; J + +の）
（は（私= 0 ;私< = matsize ;私+ + ）
（もし（無し！ =私）
（ ybus （私は、無し） = ybus （幼、私） =- 1.0 * Yに変更（私は、無し） ; ）
）
CNTの+ + ;
）
/ /終了
/ /最後の作成yバス
/ /バス0ではない表に含まれてybus
/ /ため、これは、リファレンスバス
を（無し= 1 ;無し< = matsize ; J + +の）
（は（私= 1 ;私< = matsize ;私+ + ）
（ ybus （私- 1 、のJ - 1 ） = ybus （私は、無し） ;
）
）
ybus.setsize （ matsize 、 matsize ） ;
zbus.setsize （ matsize 、 matsize ） ;
printf （ " \ n \ nybus表\ n " ） ;
裁判所< < ybus < < endl ;
zbus = ！ ybus ; / /入手して、逆ybus
printf （ " \ n \ nzbus表\ n " ） ;
裁判所< < zbus < < endl ;
getch （ ） ;
）