这个程序是什么功能

int32 i=0,j=0;
float32 x[10][30],a[10][30],b[10][30],c[10][30];
for(j=1;j<=n;j+=2)
for(i=1;i<=m;i++)
{
x[j][i]=(y[j][i+1]+y[j][2*m-i])/2-
(y[j][i]+y[j][2*m+1-i])/2;
a[j][i]=k/x[j][i];
b[j][i]=y[0][i]-a[j][i]*x[j][i];
}
for(i=0;i<=m;i++)
{a[0][i]=y[0][i];
b[0][i]=y[0][i];}
for(j=1;j<=n;j++)
{a[j][0]=y[j][0];
b[j][0]=y[j][0];}
for(j=2;j<=n;j+=2)
{for(i=1;i<=2*m;i++)
a[j][1]+=y[j][i];
a[j][1]=a[j][1]/(2*m);}
for(j=2;j<=n;j+=2)
{c[j][1]=a[j+2][1]-a[j][1];
a[j][1]=(a[j+2][0]-a[j][0])/c[j][1];
b[j][1]=a[j][0]-a[j][1]*c[j][1];}