|
（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）），“，”）
| 