#include <stdio.h>
#include <stdlib.h>
#include <time.h>

#define N 1024

void keisan(double a[N][N],double b[N][N],double c[N][N]);

int main()
{
        int i,j;
        static double A[N][N],B[N][N],C[N][N];

        srand(time(NULL));

        for(i=0;i<N;i++) {
                for(j=0;j<N;j++){
                               A[i][j]=(double)rand()/(RAND_MAX);
                               B[i][j]=(double)rand()/(RAND_MAX);
                               C[i][j]=0.0; 
                }
        }

        keisan(A,B,C); 

         printf("end\n");
        return 0;
}

void keisan(double a[N][N],double b[N][N],double c[N][N])
{
        int i,j,k;

        for(i=0;i<N;i++) {
                for(j=0;j<N;j++) {
                            c[i][j]=0.0;
                        for(k=0;k<N;k++) {
                                c[i][j]+=a[i][k]*b[k][j];
                        }
                }
        }
}