结构分析的有限元法与MATLAB程序设计 里面k=stiffnessmatrix(ie,1)这段程序stiffnessm

结构分析的有限元法与 MATLAB程序设计 里面k=stiffnessmatrix(ie,1)这段程序stiffnessmatrix什么意思怎么用
就是运行不了，不知道甚么原因，这是程序：
main()
{int nj,ne,nz,i,j,je,i1,i2,i3,i4,i8,j3,j4,j8,njf;
int ie[50][3],zj[20][5],z[15][5];
float vm,p[8]={0,0,0,-5000,0,5000,0,0},t=0.1,a[1],b[1][2],c[1][2],g[1][1],aa[1],zk[15][10];
double e=210.e9;
nj=4;ne=2;nz=4;vm=0.33;
for(i=0;i<nj;i++)
for(j=0;j<=1;j++)
{ scanf("%d",&z [j]);
printf("z[j]=%d",z[j]);
printf("\n");} /*坐标数组*/
for(i=0;i<nz;i++)
for(j=0;j<=1;j++)
{ scanf("%d",&zj[j]);
printf("zj[j]=%d",zj[j]);
printf("\n");} /*约束数组*/
for(i=0;i<ne;i++)
for(j=0;j<3;j++)
{ scanf("%d",&ie[j]);
printf("%d",ie[j]);
printf("\n");} /*布尔矩阵*/
for(je=0;je<ne;je++) /*单刚形成整刚*/
{i1=ie[je][0]-1;
i2=ie[je][1]-1;
i3=ie[je][2]-1;
b[je][0]=z[i2][1]-z[i3][1];
b[je][1]=z[i3][1]-z[i1][1];
b[je][2]=z[i1][1]-z[i2][1];
c[je][0]=z[i3][0]-z[i2][0];
c[je][1]=z[i1][0]-z[i3][0];
c[je][2]=z[i2][0]-z[i1][0];
aa[je]=a[je]=0.5*(b[je][0]*c[je][1]-b[je][1]*c[je][0]);
a[je]=e*t/(4.0*(1.0-vm*vm)*a[je]);
{ for(i=0;i<3;i++)
{ for(j=i;j<3;j++)
{ g[0][0]=b[je]*b[je][j]+0.5*(1.0-vm)*c[je]*c[je][j];
g[1][1]=c[je]*b[je][j]+0.5*(1.0-vm)*b[je]*b[je][j];
g[1][0]=vm*b[je][j]*c[je]+0.5*(1.0-vm)*b[je]*c[je][j];
g[0][1]=vm*c[je][j]*b[je]+0.5*(1.0-vm)*c[je]*b[je][j];
for(i3=0;i3<2;i3++)
{ if(i>=j)
j8=i3;
else j8=0;
for(j3=j8;j3<2;j3++)
{i4=2*(ie[je]-1)+i3;
j4=2*(ie[je][j]-1)+j3;
if(i4>j4)
{i8=i4;i4=j4;j4=i8;}
zk[i4][j4]+=g[i3][j3]*a[je];
}}}}}}
for(i=0;i<nz;i++) /*置大数法*/
{ j= zj[1];
zk[j][j]=zk[j][j]*1.0e25;
p[j]=zk[j][j]*zj[1][2];
}
}