Matpower计算潮流问题，程序如下

Huntsman

程序：
; j5 x2 H- c% Mfunction [baseMVA,bus,gen,branch]=casefz

baseMVA = 100;

bus =[
    1.0000    3         0         0         0         0    1.0000    1.0000         0  116.0000    1.0000    1.1000    0.9000
) _% @4 L. `* K, }    2.0000    4         0         0         0         0    1.0000    1.0000         0  115.5000    1.0000    1.1000    0.9000
    3.0000    4         0         0         0    1.0000    1.0000         0         0  110.0000    1.0000    1.1000    0.9000
( I8 \3 e- y1 T7 U2 P8 r    4.0000    1   20.0000   15.0000         0         0    1.0000    1.0000         0   11.0000    1.0000    1.1000    0.9000
. T6 M1 }4 m* u" }    5.0000    1   10.0000    8.0000         0         0    1.0000    1.0000         0    6.3000    1.0000    1.1000    0.9000
  s& i$ t* C( b! J# E# u9 i    6.0000    1    8.0000    6.0000         0         0    1.0000    1.0000         0   11.0000    1.0000    1.1000    0.9000];
% N/ A) S& p) s5 |$ Wgen =[
& P0 s( W+ R+ n     1     0     0    30   -30     1   100     1   250     0
     5    24    18    30   -30     1   100     1   300     0];
+ b" S7 l% O7 q% I
branch =[
    1    2   16.8800   17.1600    0.0004  250.0000  250.0000  250.0000         0         0    1.0000
& |% \) J- m$ ~    2    4    2.2400   43.7800         0  250.0000  250.0000  250.0000   10.5000         0    1.0000
$ g* ?% f( o$ D$ e2 Y    2    3    21.1000   21.4500    0.0000  150.0000  150.0000  150.0000         0         0    1.0000
) F8 D0 J7 k3 s% R* Z    3    5    4.3200   76.8000         0  300.0000  300.0000  300.0000   19.2000         0    1.0000
    3    6    29.7400  175.1500         0  150.0000  150.0000  150.0000   10.0000         0    1.0000];

return;
: F: t1 i+ ^$ a! ^. d
3 f" B. o/ L2 U7 v. m" k- j! p
报错：
  {: x% F/ x5 k2 C9 f
1 v3 a4 t8 c, p/ d

报错