function hhu_dcpf
disp('  ')
disp('  ')

disp('                ***已經選擇直流潮流法計算，請選擇算例***    ')
disp('                                                                   ')
disp('                        [1].IEEE --9 節點  ')
disp('                －－－－－－－－－－－－－－－－－－－－－－－     ')
disp('                        [2].IEEE --14 節點  ')
disp('                －－－－－－－－－－－－－－－－－－－－－－－     ')
disp('                        [3].IEEE --30 節點 ')
disp('                －－－－－－－－－－－－－－－－－－－－－－－     ')
disp('                        [4].IEEE --39 節點  ')
disp('               －－－－－－－－－－－－－－－－－－－－－－－－－  ')
disp('                        [5].IEEE --57 節點  ')
disp('                －－－－－－－－－－－－－－－－－－－－－－－     ')
disp('                        [6].IEEE --118 節點 ')
disp('                －－－－－－－－－－－－－－－－－－－－－－－     ')
disp('                        [7].IEEE --300 節點  ')
disp('               －－－－－－－－－－－－－－－－－－－－－－－－－  ')
disp('                        [8].返回上級菜單 ')
disp('               －－－－－－－－－－－－－－－－－－－－－－－－－  ')
disp('  ')
 

mpopt = mpoption('PF_DC', 1);
choice = input('   ***\輸入您的選擇 （[1--8]之間）***:       ');

if isempty(choice)
   disp('                 ')
   disp('                 ')
    disp('                      **** 錯誤： 輸入不能為空，請重新選擇***      ');
    hhu_dcpf
end


switch choice
    case 1
    runpf('case9', mpopt);
    case 2
       runpf('case14', mpopt);
    case 3
        runpf('case30', mpopt);
    case 4
       runpf('case39', mpopt);
    case 5
       runpf('case57', mpopt);
    case 6
        runpf('case118', mpopt);
    case 7
        runpf('case300', mpopt);
     case 8
         hhu_runpf
    otherwise
      hhu_agains;  
end

 
 
 return