采用子函数处理一个未知大小的矩阵 例如矩阵的加法 咋弄啊

juzhenchengfa(float *A1,float *A2,float *B1,float *B2,float *C1,float *C2,int a,int b)
/*为矩阵(A1*jA1)(a*b阶)与矩阵(B1*jB1)(b*a阶)的乘法结果为(C1*jC1)(a*a阶)*/
{
    static float d=0.0,t=0.0;
    for(i=1;i<=a;i++)
    {
        for(j=1;j<=b;j++)
            for(k=1;k<+a;k++)
            {
                d=d+A1[i][k]*B1[k][j]-A2[i][k]*B2[k][j];
                t=t+A1[i][k]*B2[k][j]-A2[i][k]*B1[k][j];
            }
            C1[i][j]=d;
            C2[i][j]=t;
    }
}