求助，MATLAB程序运行时出错，怎么改？

langia

各位前辈们，有谁特别熟悉MATLAB软件？帮我看下这个程序，为什么运行的时候，总是提示我这个错误：
Warning: Divide by zero.
> In lllllllll at 17
Warning: Divide by zero.
. K( Y2 x& \9 D( H> In lllllllll at 18
. D. x2 `2 G* m, ~??? Attempted to access G(0,0); index must be a positive integer or logical.
; _( f& n( w" K' {$ i3 Z' {Error in ==> lllllllll at 29
    G(i,i)=G(i,i)+gij/(k^2);


# |  v% N2 ^! E, V* A5 c& ^程序如下：
& r6 q0 b* B/ H1 b9 Q%牛顿——拉夫逊法进行潮流计算
n=input('请输入节点数：n=')
nl=input('请输入支路数：nl=')
isb=input('请输入平衡母线节点号：isb=')
$ W' {! L5 q  D7 [9 JA=input('请输入由支路参数形成的矩阵： A=')
A=zeros(nl,5);
0 W% B7 S3 }: Z9 I! ~! w( c& u: iG=zeros(n)
B=zeros(n)
9 L* ?2 K, A# f. Jfor a=1:nl
i=A(a,1);
* S2 Q2 u, I, U$ xj=A(a,2);
r=A(a,3);
x=A(a,4);
k=A(a,5);
z=r*r+x*x;
% Q6 T5 }  c; ], @6 Qgij=r/z;
; a3 Z% u/ E6 D5 f9 Jbij=-x/z;
( N4 Z! T1 N) `       if(k<0)
          k=abs(k);
        G(i,i)=G(i,i)+gij;
8 a- Q3 P# W3 m$ Z        B(i,i)=B(i,i)+bij+k;
        G(j,j)=G(j,j)+gij;
! r+ G( L& w' W9 S  a3 W9 ~+ C        B(j,j)=B(j,j)+bij+k;
        G(i,j)=-gij;
2 [" W$ C7 q2 W2 D& S* u& v) @        B(i,j)=-bij;
3 z% f& G  |% d2 i! Q1 Selseif(k==0)
! }: G5 `5 N: k    k=1;
8 X- |8 b+ U0 e7 |    G(i,i)=G(i,i)+gij/(k^2);
    B(i,i)=B(i,i)+bij/(k^2);
5 Z/ b& C3 g7 w* v. |* Y( c4 o    G(j,j)=G(j,j)+gij;
    B(j,j)=B(j,j)+bij;
$ _1 R8 A) ?) o. y    G(i,j)=-gij/k;
    B(i,j)=-bij/k;
elseif(k>0)
% s; |6 B+ S" b7 V    G(i,i)=G(i,i)+gij/(k*k);
' b% L( x+ ~0 C0 E$ _* f0 N% i7 s    B(i,i)=B(i,i)+bij/(k*k);
        G(j,j)=G(j,j)+gij;
* d& x% I& g# k# m        B(j,j)=B(j,j)+bij;
        G(i,j)=-gij/k;
    B(i,j)=-bij/k;
' b* l0 H' l7 l  e7 N- ]end
% @+ T# n$ n8 C8 K% g1 e7 P* nif(i>j)
           G(j,i)=G(i,j);
           B(j,i)=B(i,j);
end
end
for a=1:n
' V: g9 \2 V0 f; D. y0 P     for b=1:n
        if(a<b)
' C+ s2 O$ [) x7 v, H) n            G(b,a)=G(a,b);
            B(b,a)=B(a,b);
* Q/ ~& k( X% {) ^% q        end
6 p1 b7 I7 _$ ^, G  V: r: r) B     end
end
1 @' Z& o% T& |n=input('请输入节点数：n=')
% v" ?" p* k( L3 Bnl=input('请输入支路数：nl=')
/ G" N+ B! [) N6 Q4 _isb=input('请输入平衡母线节点号：isb=')
- Y- h- y$ D4 V; Y6 @g=0%迭代次数
: P  u; z5 F, y: g1 Z- g/ M9 Ce=input('请输入各节点参数形成的矩阵： e=')%%%%%%%%%%原始数据输入
f=input('请输入各节点参数形成的矩阵： f=')
Pis=input('请输入各节点参数形成的矩阵： Pis=')
1 N( _% O, v1 M/ `9 mQis=input('请输入各节点参数形成的矩阵： Qis=')
# V, P0 d3 x$ z1 N# d e=zeros(n,1)
f=zeros(n,1)
Pis=zeros((n-1),1)
( w: r, }) ~6 S0 |+ I Qis=zeros((n-1),1)
max=1
while(max>1.0e-5)
A=input('请输入由支路参数形成的矩阵： A=')
G=zeros(n,n)
B=zeros(n,n)
for a=1:nl
    i=A(a,1)
* L# V, y& Z( T8 M  r6 [' Ij=A(a,2)
* P2 r! [9 l  w3 |  m7 a# ^r=A(a,3)
- q( c: O# Y' K1 T) j7 P2 fx=A(a,4)
k=A(a,5)
z=r*r+x*x
gij=r/z
, M& I7 t# v) P* K" E/ u       bij=-x/z
       if(k<0)
8 _" `$ {4 \  M* M" R          k=abs(k);
        G(i,i)=G(i,i)+gij;
1 T5 G# v7 ?7 c1 Z1 J        B(i,i)=B(i,i)+bij+k;
        G(j,j)=G(j,j)+gij;
( j) ]. o& t' Q0 t, L+ b        B(j,j)=B(j,j)+bij+k;
/ \5 r/ q1 N! s$ G# [) K$ d        G(i,j)=-gij;
8 b# m2 L1 b8 s. o, P8 }- I4 }+ P, F            B(i,j)=-bij;
4 @& n2 m6 e; Z6 q: _elseif(k==0)
0 p3 i5 e. B2 [  A, C    k=1;
; x) n  r& Q6 t" I3 j( IG(i,i)=G(i,i)+gij/(k*k);
* C# b' z4 v( {/ s    B(i,i)=B(i,i)+bij/(k*k);
        G(j,j)=G(j,j)+gij;
! r/ [5 R, G" a/ H3 d2 A8 T        B(j,j)=B(j,j)+bij;
        G(i,j)=-gij/k;
    B(i,j)=-bij/k;
7 F- u' E& A& K) l8 a! O- ]elseif(k>0)
$ u7 o1 U2 ~# |8 p8 D    G(i,i)=G(i,i)+gij/(k*k);
9 ]5 ^1 ]) q: j5 D    B(i,i)=B(i,i)+bij/(k*k);
$ L9 J2 z# h9 `# l8 e1 s5 F$ {        G(j,j)=G(j,j)+gij;
  j$ z8 T$ S) {, i% I  T        B(j,j)=B(j,j)+bij;
4 b( V: h! ~1 s6 K- y/ T0 v. F0 v( i        G(i,j)=-gij/k;
' C9 \+ H# ~/ M6 T    B(i,j)=-bij/k;
end
if(i>j)
, N- V( x0 o8 i& e* h7 L8 U           G(j,i)=G(i,j);
           B(j,i)=B(i,j);
end
' k  J# B. Y9 L! A8 G8 D5 Fend
3 m/ _  T+ b, K' W6 ]6 v for a=1:n
     for b=1:n
3 p8 w5 N" }) u3 ?+ v3 I+ a        if(a<b)
% V% c8 z: W( k* X4 s; f- y            G(b,a)=G(a,b);
            B(b,a)=B(a,b);
        end
     end
end
7 j0 P7 \6 k9 N0 H这是求节点导纳矩阵。请问，为什么我输进去参数之后，老是提示那个错误？？
3 A5 a: h+ m2 c: `/ G, f  k# T7 ~还有，定义矩阵的时候，定义完之后要不要加;号？？到底怎么区分？？急死了。
大家快来帮帮我。。。。。。。。。。。。。。。。。。
% f3 g0 x0 X9 z. ~先谢谢啦。。。。。。。

qiuwei_002

能否提供一下参数？

sadever

A=input('请输入由支路参数形成的矩阵： A=')" i9 a, L# ^- E) y1 j5 e  E/ Y
A=zeros(nl,5);2 m- H+ g% m- c2 Y
/ z$ Z) v1 G9 i; L% H9 e- V  A  r( G这两句是不是有问题，既然赋值了，为什么又给零矩阵。
( v6 R5 p3 v0 S. s) j这两个错误是很简单的错误，一个是除零，一个是矩阵下标为零，你只要给个断点就可以了。

hhusunguoqiang

楼上的说的很对，做除法的时候有0元素出现！