(PARI)a33840_1_2(n)={my(md=3*n*(n+1)/2+1,
T1=矩阵(n,n),T2=矩阵,
M=矩阵(md,md,i,j,0),U=向量(md),
有效(i,j)=i>0&&i<=n&j>0&j<=n&i>=j,k=0,neq=1);
\\边缘列表
对于(i=1,n,对于(j=1,i,T1[i,j]=k++;T2[i,j]=k++;T3[i,j]=k+));
\\所有节点的电流流入和流出
\\左下角三角形,电流从电压源流入
M[1,1]=-1;M[1,2]=-1;M[1,md]=1;
\\三角形左下角上的循环
对于(i=2,n+1,对于(j=1,i,
\\排除三角形顶部的节点
如果(j<n+1,neq++;
如果(有效(i-1,j),M[neq,T1[i-1,j]]=1;M[neq,T3[i-1,j]]=1);
如果(有效(i-1,j-1),M[neq,T2[i-1,j-1]]=1;M[neq,T3[i-1,j-1]]=-1);
如果(有效(i,j),M[neq,T1[i,j]]=-1;M[neq,T2[i,j]]=-1);
\\右下角,电流通过电压源
如果(i==n+1&&j==1,M[neq,md]=-1)
)));
\\顶点高于底面的三角形周围的电压之和
对于(i=1,n,对于(j=1,i,
neq++;M[neq,T1[i,j]]=1;M[neq,T2[i,j]]=-1;M[neq,T3[i,j]]=-1));
\\顶点低于底面的三角形周围的电压之和
对于(i=1,n-1,对于(j=1,i,neq++;
M[neq,T3[i,j]]=1;M[neq,T2[i+1,j]]=1;M[neq,T1[i+1,j+1]]=-1));
\\施加在三角形下角的外部电压
neq=neq++;对于(i=1,n,M[neq,T1[i,1]]=1);
\\方程式右侧;驱动电压1伏
U[neq]=1;1/矩阵解(M,U~)[neq]};
对于(n=1,10,打印1(分子(a33840_1_2(n)),“,”)