Program v C

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

/*  Delo z matrikami: vnos, izpis, množenje.  */

static double **matrika(int, int);
static double **vnos(int *, int *);
static void izpis(double **, int, int);
static double **produkt(double **, int, int, double **, int, int);

int main(void)
{
   double **a, **b, **c;
   int ma, na, mb, nb;

   if ((a = vnos(&ma, &na)) == NULL) return 1;
   if ((b = vnos(&mb, &nb)) == NULL) return 1;
   if ((c = produkt(a, ma, na, b, mb, nb)) == NULL)
      printf(
         "\aNapacni podatki ali premalo pomnilnika.\n"
         "Matrik ne morem zmnoziti.\n"
      );
   else {
      printf("Produkt danih matrik je:\n");
      izpis(c, ma, nb);
   }
   free(a);  free(b);  free(c);
   return 0;
}

/* Pripravi prostor za matriko velikosti m x n.  Če ne uspe, vrne NULL. */
static double **matrika(int m, int n)
{
   double **p, *q;
   int i;

   if ((p = malloc(m * sizeof(q) + m * n * sizeof(*q))) != NULL)
      for (i = 0, q = (double *)&p[m]; i < m; q += n) p[i++] = q;
   return p;
}

/* Zgradi matriko: prebere elemente, vrne kazalec nanje,    */
/* v m in n zapiše velikost matrike. Če ne uspe, vrne NULL. */
static double **vnos(int *m, int *n)
{
   double **a;
   int i, j;

   printf("Vnos matrike:\n");
   printf("Stevilo vrstic: ");  scanf("%d", m);
   printf("Stevilo stolpcev: ");  scanf("%d", n);
   if ((a = matrika(*m, *n)) == NULL) {
      printf("\aPremalo pomnilnika.\n");
      return NULL;
   }
   printf("Vnesi vse elemente (%d) matrike.\n", *n * *m);
   printf("Locis jih lahko s poljubnimi 'belimi' znaki.\n");
   for (i = 0; i < *m; ++i)
      for (j = 0; j < *n; ++j)
         scanf("%lf", &a[i][j]);
   return a;
}

/* Izpiše elemente matrike a, ki je velikosti m x n. */
static void izpis(double **a, int m, int n)
{
   int i, j;

   for (i = 0; i < m; ++i)
      for (j = 0; j < n; ++j)
         printf("%g%c", a[i][j], (j < n - 1) ? '\t' : '\n');
}

/* Zmnoži matriki a in b. Če ne uspe, vrne NULL. */
static double **produkt(double **a, int m, int n, double **b, int p, int q)
{
   double **c;
   int i, j, k;

   if (a == NULL || b == NULL) return NULL;           /* Napačni podatki. */
   if (n != p) return NULL;        /* Matriki nista usklajenih velikosti. */
   if ((c = matrika(m, q)) == NULL) return NULL; /* Ni dovolj pomnilnika. */
   for (i = 0; i < m; ++i)
      for (j = 0; j < q; ++j) {
         c[i][j] = 0;
         for (k = 0; k < n; ++k)
            c[i][j] += a[i][k] * b[k][j];
      }
   return c;
}