潮流计算中导纳矩阵的计算程序

潮流计算中导纳矩阵的计算程序n=input('请输入节点数:n=');nl=input('请输入支路数:nl=');pr=input('请输入误差精度：pr=');b=input('请输入由支路参数形成的矩阵b=');x=input('请输入由节点号及其对地阻抗形成的矩阵：x=');y=zeros(n);