Рекурсия: вычисление матричного многочлена

#include "stdafx.h"
#include <stdio.h>
#include <conio.h>
#include <math.h>
#include "locale.h"
const int N=100;
 
void func(int n, int m,int A[N][N], int E[N][N], float p[N])
{
    static int k,i,j;
    static float itog[N][N],g[N][N],h[N][N];
    if(k==n) 
    {
        printf("Результат:\n");
        for(i=0;i<m;i++)
        {
            for(j=0;j<m;j++)
            printf("%6.2f", itog[i][j]);
            printf("\n");
        }
    }
    else
    {    
        if(k==0)
        {
            k++;
        }
        else
        {
            if (k==1)
            {
                for(i=0;i<m;i++)
                    for(j=0;j<m;j++)
                    g[i][j]=p[k-1]*A[i][j];
                for(i=0;i<m;i++)
                    for(j=0;j<m;j++)
                    h[i][j]=p[k]*E[i][j];
                for(i=0;i<m;i++)
                    for(j=0;j<m;j++)
                        itog[i][j]=g[i][j]+h[i][j];
                k++;
            }
            else 
            {
                for(i=0;i<m;i++)
                    for(j=0;j<m;j++)
                    g[i][j]=itog[i][j]*A[i][j];
                for(i=0;i<m;i++)
                    for(j=0;j<m;j++)
                    h[i][j]=p[k]*E[i][j];
                k++;
            }
        }
        func(n,m,A,E,p);
    }
}
int main()
{ 
    setlocale(LC_ALL,"Russian");
    unsigned int m, n, i, j;
    int E[N][N], A[N][N];
    float p[N];
 
    printf("Введите размерность матрицы А[m][m](от 1 до %d):",N);
    scanf("%d",&m);
    printf("\nВведите элементы матрицы А[m][m]:\n");
    for (i=0; i<m; i++)
        for (j=0; j<m; j++)
        {
            printf("A[%d][%d]=",i,j);
            scanf("%d",&A[i][j]);
        }
    printf("\nВведите n(от 1 до %d):",N);
    scanf("%d",&n);
    printf("\nВведите элементы p n, p n-1, ..., p 0:\n");
    for (i=0; i<=n; i++)
    {
        printf("p%d=",n-i);
        scanf("%f",&p[n]);
    }
 
    for(i;i<m;i++)
    {
        for(j;j<m;j++)
            if(i == j)
            {
                E[i][j]=1;
            }
            else
            {
                E[i][j]=0;
            }
         }
    printf("\n\n");
    func(n,m,A,E,p);
    _getch();
    return 0;
}