#define yes 1
#define no  0

// 状態変数の添字
Integer    Ndlt1, Nw1, NeqD1, NgV1, NVvm1;

//  パラメータ（システム１，システム２）
Real    Xd1, Xq1, XdD1, Vinf1, M1, D1, Vts1, efs1, UlV1, LlV1;
//Integer    NPms1, NGV1, NTdoD1, NTV1, NTvm1;
Real	Pms1, GV1, TdoD1, TV1, Tvm1;
Real    Re1, Xtrans1, Xline1;
Real    a11, a21, a31, a41;
Real    c01, c11, c21, c31, c41, c51, c61;
// 変数（システム１，システム２）
Real    id1, iq1, ed1, eq1, edinf1, eqinf1, Pm1, Pe1, Q1, ef1, Vt1;
//  パラメータ（システム１とシステム２に共通）
Real    wB;
Real    fault_time;		// 三相短絡事故が生じる時刻
Real    fault_duration;	// 事故継続時間
Real    reclose_interval;	// 事故除去から再閉路までの時間
Real    twoline_rate;		// 一回線と二回線ある内、二回線の比率
// 二回線上での三相短絡事故が生じる位置 (一回線との接続点からの距離)
Real    fault_position;

////  送電線路図
////                                    2aZe
////                                   ______
////              (1-a)Ze  +-----------|    |---------+
////              ______   |           +----+         |
////  o-----+-----|    |---+                          +------+--------o
////        |     +----+   |  S2   ______   ______    |      |
////        |              +-o/ o--|    |-+-|    |----+      |
////       ---                     +----+ | +----+          ---
////       --- jBe                 2abZe  o  2a(1-b)Ze      --- jBe
////        |                              \                 |
////        |                             o  S1              |
////      -----                           |                -----
////      /////                         -----              /////
////                                    /////
////    fault_time : 三相短絡事故が生じる時刻  i.e. S1=ON,  S2=ON 
////    fault_duration : 事故継続時間          i.e. S1=ON,  S2=OFF
////    reclose_interval : 事故除去から再閉路までの時間    i.e. S1=OFF, S2=ON
////    twoline_rate : 一回線と二回線ある内、二回線の比率     i.e. "a"
////    fault_position : 二回線上での三相短絡事故が生じる位置  i.e. "b"
////

//  コマンドラインオプション
////// Integer opt_f opt_b opt_s opt_c opt_p opt_e opt_x opt_i opt_v;
////// Integer opt_i_arg opt_v_arg;

Integer    comp_dim;	//  計算に用いる次元

Real    h,  t_fin;
Matrix  x_ini,  eps,  xbase;


Integer systype;
#define WithVMdelay 0
#define FirstOrderAVR 1

// Integer    draw_mode;		// 相関図の縦軸，横軸を何にとるか

Integer pid;
Integer auto_clear, update_interval;

// 関数の宣言
Matrix ftrs5();
Matrix Dftrs5(); // 変分系
void set_parameters();
void rem_prm_output();
List kbd_event();