program ymnoz_matr_on_matr;

{$APPTYPE CONSOLE}

uses
  SysUtils,
  windows;

const
  n=3;
  m=3;
  r=3;
type
  TElem=integer;
  TMatrixA=array [1..n,1..n] of TElem;
  TMAtrixB=array[1..m,1..m] of TElem;
  TMatrixC=array[1..r,1..r] of TElem;
var
  A: TMatrixA;
  B: TMAtrixB;
  C: TMAtrixC;

procedure input_matrixs (var A: TMatrixA;n: TElem;var B: TMAtrixB;m: TElem);
var i,j,k,l: integer;
begin
  writeln('введите матрицу А');
  for i:=1 to n do
  begin
    for j:=1 to n do
    read(A[i,j]);
    readln
  end;
  writeln('введите матрицу B');
  for k:=1 to m do
  begin
    for l:=1 to m do
    read(B[k,l]);
    readln
  end;
end;

procedure Ym_matr_on_matr(A: TMatrixA;B: TMatrixB;var C: TMatrixC;r: TElem);// умножение матриц
var p,f: integer;
//---------
 function Scalar_Vector_A_on_Vector_B (Row,Col: integer): TElem;//умножение строки Row матрицы А и столбца Col матрицы B
 var i,j: integer;
 begin
 result:=0;
 i:=0;
 for j:=1 to n do
 inc(i);
 result:=A[Row,j]*B[i,Col]+result;
 end;

begin
  for p:=1 to r do
  begin
    for f:=1 to r do
    C[p,f]:= Ym_matr_on_matr(i,j);
end;

procedure print_matrixC (C: TMatrixC;r: TElem);
var
  p,f: integer;
begin
  writeln('матрица С');
  for p:=1 to r do
  begin
    for f:=1 to r do
    write(C[p,f],' ');
    writeln
  end;
end;

begin
SetConsoleCP(1251);
SetConsoleOutputCP(1251);
input_matrixs(A,n,B,m);
Ym_matr_on_matr(A,B,C,r);
print_matrixC(C,r);
readln
end.

Цитата(Insomnia @ 6.06.2011 20:57)

Uses Crt;
const Nmax=10; //ранг матрицы не больше 10
type matrice=array[1..Nmax,1..Nmax] of integer;
var a,c,b:matrice; i,j:byte; N:1..Nmax;
procedure Multiply(var a,b:matrice);
var i,j,k:byte;
begin
for i:=1 to N do
   for j:=1 to N do
    for k:=1 to N do
   C[i,j]:=C[i,j]+a[i,k]*b[k,j]; //умножаем
end;
procedure Print(var a:matrice); //выводим на экран
var i,j:byte;
begin
for i:=1 to N do  
 begin
   for j:=1 to N do
    write(a[i,j],' ');
    writeln;
  end;
end;
Begin ClrScr;
readln(N);
writeln('A');
for i:=1 to N do
    for j:=1 to N do
       readln(a[i,j]);
writeln('B');
for i:=1 to N do
    for j:=1 to N do
       readln(b[i,j]);
Multiply(a,b);
Print(c);
readln;
end. 

const
  n=3;
  m=3;
  r=3;
type
  TElem=integer;
  TMatrixA=array [1..n,1..n] of TElem;
  TMAtrixB=array[1..m,1..m] of TElem;
  TMatrixC=array[1..r,1..r] of TElem;
var
  A: TMatrixA;
  B: TMAtrixB;
  C: TMAtrixC;

const
   n = 3;
type
   TElem = integer;
   TMatrix = array[1 .. n, 1 .. n] of TElem;
var
   A, B, C : TMatrix;

procedure Mult(A, B : TMatrix; var C : TMatrix);

   function ScalarVecA_VecB (RowA, ColB : integer) : TElem;
   var i : integer;
   begin
      result := 0;
      for i := 1 to n do
         result := result + A[RowA, i] * B[i, ColB];
   end;

var
   i, j : integer;

begin
   for i := 1 to n do
      for j := 1 to n do
         C[i, j] := ScalarVecA_VecB(i, j);
end;