function hhu_runpf
disp('          ***歡迎使用Power Flow Calculation Packages進行潮流計算，請選擇潮流計方法***    ')
disp('                 ')
disp('                 ')
disp('              －－－－－－－－－－－－－－－－－－－－－－－－－－   ')
disp('                        [1].直 流 法 潮 流 計 算                                                                         '   )
%disp('|參考文獻: H. Wang, C. E. Murillo-Sánchez, R. D. Zimmerman, R. J. Thomas, “On Computational Issues of                  |')
%disp('| Market-Based Optimal Power Flow”, IEEE Transactions on Power Systems, Vol. 22, No. 3,Aug. 2007, pp. 1185-1193.    \   |')
disp('              －－－－－－－－－－－－－－－－－－－－－－－－－－   ')


disp('                        [2].牛 頓 拉 夫 遜 法 潮 流 計 算   ')
%%disp('        參考文獻:  W. F. Tinney and C. E. Hart, “Power Flow Solution by Newton’s Method”, IEEE Transactions on')
%%disp (' Power Apparatus and Systems, Vol. PAS-86, No. 11, Nov. 1967, pp. 1449-1460. ')
disp('              －－－－－－－－－－－－－－－－－－－－－－－－－－   ')


disp('                        [3].快 速 PQ 分 解 法 潮 流 計 算   ')
%%disp('        參考文獻:  B. Stott and O. Alsac, “Fast decoupled load flow”, IEEE Transactions on Power Apparatus and')
%%disp (' Systems, Vol. PAS-93, June 1974, pp. 859-869. ')

disp('              －－－－－－－－－－－－－－－－－－－－－－－－－－   ')
disp('                        [4].Gauss-Seidel 潮 流 計 算   ')
%%disp('        參考文獻: A. F. Glimm and G. W. Stagg, “Automatic calculation of load flows”, AIEE Transactions (Power')
%%disp (' Apparatus and Systems), vol. 76, pp. 817-828, Oct. 1957. ')
disp('              －－－－－－－－－－－－－－－－－－－－－－－－－－   ')
disp('                        [5].退 出 程 序 運 行   ')
%%disp('        參考文獻: A. F. Glimm and G. W. Stagg, “Automatic calculation of load flows”, AIEE Transactions (Power')
%%disp (' Apparatus and Systems), vol. 76, pp. 817-828, Oct. 1957. ')
disp('              －－－－－－－－－－－－－－－－－－－－－－－－－－   ')
disp('                 ')
disp('                 ')

choice = input('     ****輸入您的選擇[1-5之間]*** :        ');
if isempty(choice)
   disp('                 ')
   disp('                 ')
    disp('                      **** 錯誤： 輸入不能為空，請重新選擇***      ');
   hhu_runpf
end
switch choice
    case 1
     hhu_dcpf
    case 2
     hhu_newtonpf
    case 3
      hhu_pqpf
    case 4
      hhu_gspf
    case 5 
      hhu_break;
       
       
   otherwise
      hhu_again;  
 end

return